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The  objective  of  this  work  is  to  characterize  the  acidity  H3PW12O40  in  different 
systems.  This  acid  can  catalyze  reactions  in  homogeneous  and  heterogeneous  systems. 
Acidity  and  solubility  can  be  controlled  by  substituting  the  proton  for  different  cations. 

The  calorimetric  titration  of  H3PW12O40  (H3PW)  with  a  series  of  bases  in  CH3CN 
gives  equilibrium  constants  and  enthalpies  for  the  three  deprotonation  steps.  These 
results  show  that  solutions  of  H3PW  in  any  solvent  whose  basicity  is  equal  to  or  weaker 
than  CH3CN  will  furnish  only  one  strong  proton  per  molecule  of  acid  for  catalytic 
reactions.  The  best  procedure  to  prepare  the  anhydrous  H3PW  is  established. 

The  ECW  model  gives  a  Ca*/Ea*  ratio  of  2.88,  which  indicates  significant 
covalency  in  the  interactions  of  this  acceptor.  The  enthalpies  of  reactions  of  pyridine 
with  a  series  of  Bronsted  acids  in  acetonitrile  solvent  leads  to  the  acidity  order: 
H3PW12O40  >  CF3SO3H  >  P-CH3C6H4SO3H  =  H2SO4  >  CF3COOH  >  m-ClC6H4COOH. 

Thermodynamic  constants  for  H3PW  and  CsxHs.xPW  derivatives  reacting  with 
pyridine  in  cyclohexane  solution  are  provided.    The  Cal-ad  results  indicate  that  the 
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protons  of  H3PW  in  the  solid  behave  like  those  in  solution,  with  all  having  different 
strengths.  The  enthalpy  for  reaction  of  the  most  acidic  proton  of  H3PW  with  pyridine  is 
equivalent  for  the  solid  and  in  solution  of  acetonitrile. 

Incomplete  neutralization  of  the  total  available  protons  of  H3PW,  CSH2PW,  and 
CS2HPW  with  pyridine  in  the  calorimetric  titrations  and  spectroscopic  methods  suggests 
formation  of  layered  products.  These  methods  also  support  that  the  distribution  of  the 
protons  on  the  cesium  salts  is  random  and  homogeneous  in  the  solids. 

Measurements  of  enthalpy  of  neutralization  for  the  strongest  protons  in  each 
compound  lead  to  the  following  order  of  acidity:  H3PW  >  CsaHPW  =  CS2.5H0  5PW  > 
CSH2PW  >  CS3PW. 

Thermodynamic  constants  for  H3PW  supported  on  silica,  alumina,  and  carbon 
reacting  with  pyridine  in  cyclohexane  solution  are  reported.  Loading  H3PW  up  to  25 
wt.%  on  these  carriers  does  not  increase  the  acidity  of  the  free  acid.  However,  Cal-ad 
results  show  that  H3PW  is  well  dispersed  on  these  supports. 

The  acidity  order  obtained  is:  H3PW  >  Si02-H3PW  supported  >  AI2O3-H3PW 
supported  >  Carbon-H3PW  supported. 
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CHAPTER  1 

HETEROPOLY  COMPOUNDS:  CLASSES,  STRUCTURES,  AND  RELEVANCE 

Definition  and  Historical  Aspects 
Heteropolyanions  are  polymeric  oxoanions  represented  by  the  general  formula 
[XxMmOy]''"  (x  <  m)  formed  by  the  condensation  of  more  than  two  different  oxoanions 
[1].  Free  acid  or  acid  forms  of  these  species  are  called  heteropoly  acids.  For  instance, 
reaction  of  tungstate  with  phosphate  salts  in  acidic  medium  produce  the  heteropoly  acid 
H3PW,204o: 

12W04^-  +  HP04^-  +  26H^    ^    H3PW,204o  +  I2H2O 
Of  the  elements  comprising  the  structure,  M  is  usually  molybdenum,  tungsten,  vanadium, 
niobium,  tantalum,  or  mixtures  of  these  elements  in  their  highest  oxidation  states  (d^,  d'), 
while  X  may  comprise  a  wide  variety  of  group  I  to  VII  elements  (usually  P,  Si,  As,  Sb, 
etc.)  [2]. 

The  first  heteropoly  compound,  now  known  as  ammonium  1 2-molybdophosphate, 
was  described  by  Berzelius  in  1826  as  the  yellow  precipitate  that  is  produced  when 
ammonium  molybdate  is  added  to  phosphoric  acid.  This  compound  was  introduced  into 
analytical  chemistry  by  Svanberg  and  Struve  in  1848.  The  PMoi204o^'  ion  was  used  in 
quantitative  determination  of  phosphate.  When  the  tungstosilicic  acid  and  its  salts  were 
discovered  by  Marignac  in  1862,  the  analytical  compositions  of  such  heteropoly  acids 
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were  precisely  determined.  By  1910,  about  sixty  heteropoly  compounds  had  been 
described  [2]. 

The  initial  attempts  to  understand  the  composition  of  heteropolyanions  were  based 
on  Werner's  coordination  theory.  An  advance  occurred  with  the  Miolati -Rosenheim 
theory,  according  to  which,  heteropolyanions  are  based  on  six-coordinate  heteroatoms 
with  MO4  '  or  M2O7  "  anions  as  ligands  or  bridging  groups.  In  1929,  Pauling  [3] 
proposed  that  complexes  with  12:1  ratio  (M  to  X)  had  an  arrangement  of  twelve  MOe 
octahedra,  each  one  sharing  three  comers  with  an  adjoining  octahedron,  surrounding  a 
central  XO4  tetrahedron.  Although  the  detailed  structure  was  not  completely  correct,  the 
resulting  formulas  predicted  the  correct  basicity  and  brought  attention  to  the  close  relation 
of  structural  principles  between  heteropolyanions  and  infinite  lattices.  In  1933  Keggin 
[4a,  4b]  determined  the  structure  of  H3PWi204o-5H20  by  X-ray  diffraction  showing  that 
the  anion  was  indeed  based  on  WOe  octahedral  units  linked  by  shared  edges  and  as  well 
as  comers.  Other  isomorphous  complexes  were  reported,  but  the  next  polyanion  stmcture 
appeared  in  1948  when  Evans  [5]  described  the  stmcture  of  TeMo6024^'.  A  different 
stmcture  for  the  ion  (H6P2M018O62)  was  reported  by  Lindqvist  [6a]  in  1950,  and  later  by 
Evans  [6b].  Many  other  heteropoly  compounds  were  isolated  and  their  stmctures  were 
determined.  Taking  into  account  those  species  with  different  m/x  ratios,  with  more  than 
one  type  of  heteroatom,  and  with  mixtures  of  polyatoms,  the  resulting  field  of 
heteropolyanions  is  vast,  and  several  hundreds  of  complexes  are  known  [2]. 
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Main  Structures  of  Heteropolyanions 
A  variety  of  polyanion  structures  exist  with  different  formulas.  A  few  examples 
include  Keggin  (a,  P,  y,  and  e-types;  XM12O40),  Silverton  (XM12O42),  Dawson 
(X2M18O62),  Anderson  (XM6O24),  and  Waugh  (XM9O32).  The  X  element  represents  the 
central  atom  or  heteroatom  and  M  is  called  peripheral,  poly  or  addenda  atoms. 
Heteropoly  compounds  may  be  classified  according  to  the  ratio  of  the  number  of  central 
atoms  to  the  addenda  atoms.  Compounds  with  the  same  number  of  atoms  in  the  anion 
usually  are  isomorphous  and  have  similar  chemical  properties  [7-10].  Typical 
heteropolyanions  [1 1]  are  listed  in  Table  1-1. 


Table  1-1.  Main  Classes  of  Heteropolyanions. 


Ratio  X:M 

Structure 

Polyatom  (M) 

Heteroatom  (X) 

Formula 

1:12 

Keggin 

Mo,  W 

P(V),  As(V),  Si(IV),  Ge(IV) 

[X""M,204o]^'-"'- 

1:12 

Silverton 

Mo 

Ce(IV),  Th(IV) 

[X'*^Mi2042]'- 

2:18 

Dawson 

Mo,  W 

P(V),  As(V), 

[X2'"M,8062]^- 

1:9 

Waugh 

Mo 

Mn(IV),  Ni(V) 

1:6 

Anderson 

Mo,  W 

Te(VI),  I(VII) 

[X""M6024]^"-">- 

The  Keggin  structure  (XM12O40),  illustrated  in  Figure  1-1,  is  based  on  tetrahedral 
heteroatoms  [2-4].  The  structure  has  an  overall  Ta  symmetry  and  is  formed  by  a  central 
XO4  tetrahedron  surround  by  twelve  MOe  octahedra  arranged  in  four  groups  of  three 
edge-shared  octahedra,  M3O13.  These  groups  (called  M3  triplets)  are  linked  by  shared 
comers  to  each  other  and  have  a  common  oxygen  vertex  to  the  central  heteroatom.  The 
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(b) 


Figure  1-1.  The  Keggin  Structure  (a-XMi204o)  Based  on  Tetrahedral  Heteroatoms 
Represented  by  (a)  Bond,  Polyhedral,  and  Space-filling  Drawings;  (b)  One  M3O13  Group 
is  Removed  to  Show  More  Clearly  the  Inner  Cavity  that  Accommodates  the  Heteroatom. 
The  Black  Dots  Show  the  Oxygen  Atoms  bonded  to  the  Heteroatom  [1,2]. 
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existence  of  isomers  for  the  Keggin  anion  has  been  established  [12].  In  the  P-type,  one  of 
the  three  edge  shared  W3O13  triplets  of  the  a  structure  is  rotated  by  60°  reducing  the 
symmetry  of  the  anion  from  Td  to  Csy.  The  other  isomers  involving  the  60*^  rotation  of 
two,  three,  and  all  four  W3O13  are  called  the  y,  5,  and  e-type,  respectively.  A  much  more 
complex  analysis  is  predicted  for  mixed  Keggin  anions  of  formula  XMxM  12-XO40",  which 
presents  numerous  possibilities  of  what  is  known  as  positional  isomers  [2]. 

The  Silverton  structure  PCM12O42)  is  based  on  icosahedral  coordinated 
heteroatoms  and  faced-shared  pairs  of  MOe.  This  structure  is  observed  for  molybdenum 
complexes,  but  no  analogous  tungstates  have  been  reported.  The  structure  [2]  is 
illustrated  in  Figure  1-2.    The  first  anion  described  by  the  Silverton  structure  was 

Q 

CeMoi2042  '.  The  icosahedron  of  CeOi2  is  surround  by  six  face-shared  M02O9  units, 
linked  together  by  comer-sharing.  The  structure  is  known  for  molybdates  with  Ce(ni  and 
VI),  Th(IV),  U(IV  and  V),  and  Np(V)  heteroatoms. 

The  Dawson  structure  (X2M18O62)  can  be  understood  in  terms  of  an  intermediate 
lacunary  Keggin  anion.  The  stability  of  Keggin  heteropolyanions  of  W  and  Mo  in 
solution  is  dependent  on  the  pH.  As  the  pH  increases,  X"^Mi204o**'"^'  undergoes 
hydrolysis  forming  the  lacunary  X""^M  11039^'^""^"  (pH  =  2)  which  degrades  to  give 
X"'^M903/^""^"  (pH  =  8).  Two  units  of  that  lacunary  anion  fuse  to  form  the  Dawson 
structure  (Figure  1-3).  There  are  three  isomers,  depending  on  the  number  of  rotated 
M3O13  groups  [2]. 


Figure  1-2.  The  Silverton  Structure  of  CeMoi2042*',  Which  is  Based  on  1 2-coordinated 
Heteroatoms  [2]. 


Figure  1-3.  The  Dawson  Structure  of  a-P2Wi8062^"  [2]. 
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The  Anderson  structure  (XM6O24)  is  based  upon  octahedral  heteroatoms.  The 
polyanion  arrangement  is  made  up  of  seven  edge-shared  octahedra  with  a  planar  structure 
averaging  Dad  symmetry.  Although  this  structure  was  originally  proposed  for  the 
heptamolybdate  anion  (Mo7024^"),  Lindqvist  later  demonstrated  that  the  heptamolybdate 
had  a  bent  structure  with  an  averaged  C2v  symmetry  [2].  The  first  authentic  Anderson 
structure  [13]  was  observed  for  TeMo6024^"  (Figure  1-4).  Heteroatoms  in  oxidation  states 
+2,  +3,  +4,  +6,  and  +7  have  been  observed  in  this  heteropolyanion.  With  the  lower 
oxidation  states  (+2,  +3)  the  anion  carries  six  non-acidic  protons  which  are  attached  to  the 
oxygens  of  the  central  X0(,  octahedron. 

The  Waugh  structure  (XM9O32)  may  be  constructed  from  hypothetical  XMe 
Anderson  species.  For  instance  [14],  Mn'^Mo9032^"  is  composed  by  removing  three 
alternate  MoOe  octahedra  and  placing  three  MoOe  octahedra  above  and  below  the 
MnMo3  unit  (Figure  1-5).  The  final  structure  has  a  D3  symmetry. 

Structural  Principles  for  Formation  of  Polyanions 
The  cluster  structures  discussed  before  belong  to  a  more  general  class  of 
compounds  called  metal-oxygen  anions,  which  consist  of  so-called  heteropoly  and 
isopolyanions  (polyanions  formed  by  condensation  of  only  one  kind  of  oxoanion).  These 
clusters  contain  highly  symmetrical  core  assemblies  of  MOx  imits  (M  =  V,  Mo,  W)  and 
frequently  adopt  quasi-spherical  structures  based  upon  Platonic  solids  of  topological 
interest  [15].  The  structures  of  heteropoly  and  isopolyanions  appear  to  be  governed  by 
the  electrostatic  and  radius-ratio  principles  observed  for  extended  ionic  lattices. 
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Figure  1-5.  The  Waugh  Structure  of  MnMo9032*'  [2]. 
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Thus,  it  is  easier  to  describe  these  structures  in  terms  of  assemblages  of  metal-centered 
MOx  polyhedra  that  are  linked  by  shared  comers  and  edges,  and  occasionally,  faces.  The 
ionic  radii  for  polyanion  addenda  atoms  [16]  are  given  in  Table  1-2.  In  all  cases,  the 
metal  ion  is  displaced  from  the  center  towards  a  vertex  or  edge  of  its  own  polyhedron. 


Table  1-2.  Effective  Ionic  Radii  (A)  for  oxides  (O"^  =  1.40  A). 


Ion 

Td  radius 

Oh  radius 

Observed  coordination 
number  in  polyanions 

y+5 

0.49 

0.68 

4,  5,  6,  7 

Nb"5 

0.62 

0.78 

6 

Ta"^ 

0.62 

0.78 

6 

Mo"^ 

0.55 

0.73 

4,  6,7 

0.56 

0.74 

6 

Polyoxometalates  structures  are  governed  by  two  general  principles.  First,  each 
metal  atom  occupies  an  MOx  coordination  polyhedron.  The  most  common  is  an 
octahedron  or  square  pyramid  on  which  the  metal  atoms  are  displaced  towards  those 
polyhedral  vertices  that  form  the  surface  of  the  structure.  This  is  a  result  of  M-0  n 
bonding.  The  second  principle,  the  Lipscomb  principle,  establishes  that  structures  with 
MOe  octahedra  do  not  contain  more  than  two  unshared  oxygen  atoms.  This  has  been 
rationalized  in  terms  of  the  strong  trans  influence  of  the  terminal  M-0  bonds,  which 
facilitates  dissociation  of  neutral  MO3  from  the  cluster.   However,  a  few  examples  of 
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complexes  have  been  characterized  with  anti-Lipscomb  structures  [17,  18]  suggesting  this 
rule  is  not  absolute. 

The  existence  of  heteropoly  and  isopoly  compounds  is  considered  as  a  special 
class  of  compounds  of  d"  metals  of  groups  V  and  VI  [19].  The  two  decisive  factors  are: 
(1)  cis  arrangement  of  the  multiple-bonded  ligands  in  octahedral  complexes  of  d°  metals 
in  contrast  to  their  trans  arrangement  in  complexes  of  the  d^  metals;  (2)  the  pronounced 
trans  influence  of  the  metal-ligand  multiple  bond,  leading  to  considerable  weakening  and 
lengthening  of  the  bond  with  its  trans  partner.  When  applied  to  oligomeric  or  polymeric 
oxygen  compounds,  the  first  factor  results  in  localization  of  the  n  interaction  in  -M-O-M- 
O-  chains  of  d°  metals  in  M-0  bonds  through  one,  which  leads  to  a  cis  arrangement  of  the 
M-0  multiple  bonds  in  each  of  the  M0(,  octahedra.  In  none  of  the  linked  MOe  octahedra 
of  the  anion  can  terminal  oxygen  atoms  be  found  in  a  trans  position  to  one  another.  Only 
linkage  of  octahedra  along  common  edges  is  considered,  since  bonding  by  a  common  face 
results  in  metal  atoms  being  too  close  to  one  another,  and  connection  through  a  common 
vertex  either  keeps  the  trans  arrangement  of  the  terminal  oxygen  atom  or  leads  to 
constructions  with  too  low  coordination  numbers  for  the  internal  oxygen  atoms. 
Therefore,  the  anion  that  satisfies  this  main  requirement  and  has  the  least  number  of 
metal  atoms  is  the  trimer  M3O13.  Because  the  mean  charge  per  oxygen  atom  for  this 
species  is  too  high  (-0.6  e),  this  unit  is  metastable,  and  must  form  larger  clusters 
promptly.  The  second  and  third  factors  are  applicable  to  the  distribution  in  the  trans 
positions  relative  to  M-0  multiple  bonds  (in  particular,  the  trans  position  relative  to  a 
terminal  oxygen)  of  the  oxygen  atoms  having  the  highest  coordination  number  in  relation 
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to  the  metal.  Taken  together,  all  this  leads  to  the  appearance  of  oligomeric  species 
constructed  out  of  octahedra  linked  by  common  edges  in  such  a  way  that  multiple  bonded 
oxygen  atoms  are  peripheral  and  in  the  trans  positions  to  them  with  the  internal  oxygen 
atoms  having  high  coordination  numbers. 

Structures  of  the  iso  and  heteropolyanion  type  disappear  [19]  on  going  to  the 
metals  of  Groups  VII  and  VIII,  mainly  because  the  d'  and  d^  states,  characteristic  of  these 
metals,  take  the  place  of  d°  state  for  Groups  V  and  VI.  In  addition,  the  tendency  towards 
localization  of  the  n  interaction  is  reversed,  being  delocalized  for  Groups  VII  and  VIII. 
Those  structures  are  also  not  encountered  for  the  metals  situated  to  the  left  of  Groups  V 
and  VI  because  the  trans  influence  disappears.  With  the  increasing  polarity  of  M-0 
bonds  on  moving  across  the  Periodic  Table,  the  n  interaction  is  generally  significant. 
There  is  also  an  increase  in  the  tendency  for  coordination  numbers  higher  than  6.  An 
extra  principle  must  be  considered  for  understanding  the  appearance  of  iso  and 
heteropolyanions.  It  relates  solely  to  the  ability  of  metal  of  Group  V  and  VI  to  change 
their  coordination  number  from  4  to  7  in  relation  to  oxygen.  Whatever  the  specific 
mechanism  of  polymerization  of  these  polyanions  in  solution,  they  can  proceed  only  with 
a  change  in  the  coordination  number  of  the  metal  in  the  copolymerizing  species  and  with 
the  participation  of  protons  forming  hydroxyl  and  or  water  molecules. 

Synthesis  of  Heteropoly  Acids 
Heteropolyanions  are  most  commonly  prepared  by  acidification  of  aqueous 
solutions  of  simple  oxoanions  and  the  necessary  heteroatoms.  When  the  central  atom  is 
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not  a  transition  element,  a  soluble  molybdate  or  tungstate  may  be  dissolved  with  a  soluble 
salt  containing  the  heteroatom  in  the  appropriate  oxidation  state.  Then,  the  mixture  is 
acidified  to  an  adequate  pH  range.  When  the  central  atom  is  a  transition  metal,  a  simple 
salt  of  that  element  is  mixed  hot  with  a  soluble  molybdate  or  tungstate  in  the  proper  pH 
range.  If  the  central  atom  must  be  raised  to  an  unusual  oxidation  state,  peroxide, 
persulfate  or  bromine  water  are  frequently  employed  or  electrolytic  oxidation  may  also  be 
used  [9].  Detailed  syntheses  of  heteropoly  acids  are  described  in  the  literature  [2,  9]. 

Isolation  of  the  heteropoly  acid,  after  acidification  of  an  aqueous  solution  of  the 
salt  of  the  hetero  element  and  the  alkali  metal  molybdate  or  tungstate,  is  obtained  through 
extraction  with  ether  or  ion  exchange  [2].  The  solution  of  heteropolyanion  is  shaken  with 
excess  diethyl  ether  and  three  phases  are  formed,  consisting  of  an  ether  layer,  an  aqueous 
layer,  and  a  heavy  oily  etherate.  The  lowest  layer  is  drawn  off,  shaken  with  excess  ether 
to  remove  entrained  aqueous  solution,  and  separated  again.  The  etherate  is  decomposed 
by  the  addition  of  water,  the  ether  is  removed,  and  the  aqueous  solution  of  the  heteropoly 
acid  is  allowed  to  evaporate  until  crystallization  occurs.  The  exact  nature  of  the  etherates 
[2]  is  not  clear  at  present,  although  Et20H*  has  been  detected  by  IR  spectroscopy.  Acids 
such  as  12-tungstophosphoric,  12-tungstosilicic,  and  1 2-molybdophosphoric  are 
commercially  available  as  crystalline  hydrates 

A  new  method  using  electrodialysis  has  been  reported  for  preparation  of 
H3PW12O40  with  much  better  yields  and  low  by-product  formation  [20a,  20b].  The 
heteropoly  acid  is  formed  in  the  anode  compartment  of  a  dialyzer  from  Na2W04  and 
H3PO4  in  an  aqueous  solution.  Acidification  is  provided  by  electrolysis  of  water.  Na* 
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ions  are  transferred  by  the  electric  potential  from  the  anode  to  cathode  through  a  cation- 
exchange  membrane.  As  a  result,  pure  H3PW12O40  at  the  anode  and  NaOH  at  the  cathode 
are  formed.  The  heteropoly  acid  is  isolated  by  crystallization  in  almost  100%  yield. 

Applications 

Applications  of  polyoxometalates  are  mainly  based  on  their  unique  properties, 
such  as  size,  mass,  electron-proton  transfer,  thermal  stability,  lability  of  lattice  oxygen, 
and  high  Bronsted  acidities.  They  have  uses  in  many  fields  including  analytical 
chemistry,  biochemistry,  solid  state  physics,  and  catalysis. 

Analysis  of  P,  As,  Si,  and  Ge  has  been  made  by  formation  of  XM012O40""  anions, 
which  can  be  either  precipitated  (gravimetry)  or  reduced  (colorimetry)  [21].  Variations 
contain  the  use  of  molybdate-vanadate  or  tungstate-vanadate  mixtures  [22].  A  recent 
study  uses  differential-pulse  anodic  voltametry  [23]  to  determine  simultaneously  P,  Si, 
Ge,  and  As. 

Salts  of  heteropolyanions  have  been  shown  to  be  biologically  active.  Such 
activity  includes  selective  inhibition  of  enzyme  functions,  in  vitro  and  in  vivo  antitumoral 
activity,  antiviral  activity  (e.g.,  against  rabies),  and  antiretroviral  activity  (e.g.,  against 
HIV  infections  connected  with  AIDS).  The  heteropolytungstate  NaSb307(SbW7024)'^" 
(called  HPA-23)  has  been  used  [15]  for  treatment  of  AIDS  patients  in  France.  A  few 
reasons  can  be  pointed  out  [24]  for  biological  activity  of  heteropoly  acids,  such  as  ionic 
size  and  shape,  electron  transfer  and  reservoir  properties,  stability  or  non-lability  at 
micromolar  concentrations  and  physiological  pH. 
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Many  properties  of  heteropoly  compounds  in  solution  are  also  observed  in  the 
crystalline  salts  and  free  acids,  since  the  large  anions  tend  to  have  very  weak  attractions  to 
countercations  and  solvate  molecules.  Particularly,  the  electron-  and  proton-transfer 
characteristics  are  often  conserved  in  the  solid  state.  Crystalline  acids  such  as 
H3PWi204o-21H20  are  protonic  conductors  and  proton-conducting  films  [25a-c]  have 
been  described.  Electrochromic  and  ion-selective  electrodes  have  also  been  reported 
[26a,  26b]. 

The  widest  use  of  heteropolyanions  lies  in  oxidation  and  acid  catalysis.  They  are 
also  excellent  model  systems  [27]  for  investigating  fundamental  problems  of  catalysis. 
Recent  reviews  of  catalytic  chemical  applications  of  heteropoly  compounds  in  solution 
and  solid  state  have  been  published,  emphasizing  homogeneous  oxidations  [28],  synthesis 
of  fine  chemicals  [29a,  29b],  and  acid  catalysis  [30,  31].  Table  1-3  summarizes  some 
reactions  [31]  catalyzed  by  heteropoly  compounds  in  homogeneous  and  heterogeneous 
solution. 

Solubility  of  heteropolyanions  in  water  and  many  polar  organic  solvents,  and  their 
insolubility  in  non-polar  solvents  makes  these  compounds  usefiil  for:  liquid-phase 
homogeneous  reactions  (aqueous  and  organic  media);  liquid-phase  heterogeneous 
reactions  (non-polar  solvents,  e.g.,  hydrocarbons);  vapor-phase  heterogeneous  reactions. 

The  mechanism  of  homogeneous  acid  catalysis  [29a,  29b]  by  heteropoly  acids 
involves  their  functioning  as  proton  donors.  In  the  simplest  case,  this  catalysis  can  be 
represented  by: 

S  +  HA  SH^  +  A-  (I);  SH^  ^  )    P  +  HA  (H), 
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where  S  is  the  substrate,  HA  the  acid  catalyst,  and  P  the  product.  As  stronger  acids, 
heteropoly  acids  are  able  to  activate  the  substrate  more  effectively  than  the  usual 
inorganic  acids.  However,  heteropoly  acids  play  other  important  roles  in  promoting 
reactions.  They  can  stabilize  organic  intermediates  or  cause  a  strong  salt  effect  on  the 
reaction  rate. 


Table  1-3.  Organic  Reactions  Catalyzed  by  Heteropoly  Acids. 


Reaction 


Substrate 


Isomerization 

Hydration 

Dehydration 

Ether  formation 

Esterification 

Rearrangement 

Friedel-Crafts 

Prins  reaction 

Polymerization 


1  -butene,  m-xylene,  p-xylene 
propene,  isobutene,  cyclohexene 
2-propanol,  cyclohexanol 
methanol  +  isobutene  (t-BuOH) 
carboxylic  acid  +  alcohol  (alkene) 
cumene  hydroperoxide 
alkylation,  acylation 
styrene  +  formaldehyde 
tetrahydrofuran 


Methanol  to  gasoline  conversion  (MTG)  methanol 


The  hydration  of  isobutene  is  a  known  example  [29a]  to  illustrate  those  features. 
A  two  path  mechanism  (Figure  1-6)  is  accepted  for  this  reaction  with  formation  of  the 
intermediate  carbenium  ion,  which  operates  with  mineral  acids  (path  I).  Formation  of  a 
complex  between  the  heteropoly  acid  and  the  olefin  (path  II)  is  responsible  for  promoting 
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the  reaction  and  increasing  the  order  with  respect  to  the  heteropoly  acid.  The  salt  effect, 
established  empirically,  determines  the  formal  order  of  isobutene  hydration  reaction  with 
respect  to  heteropoly  acids.  In  diluted  solutions  (0.02-0.1  M)  the  order  is  1.2-1.5,  while 
in  concentrated  solutions  (up  to  0.5  M)  is  2.0-2.2. 


(CHglj  C=CH2 


==(CH3)jC/©\ 


CH, 


(CH3)3C+ 


(CHajjCOH  +  H* 


(CH3);C^®\  CH2...SiW,20l"c 


(CH3)3C+...SiW,204o 


(CHgjaCOH  +  HSiW,20^o 


Path  I 


Path  II 


Figure  1-6.  The  Two  Path  Mechanism  for  Hydration  of  Isobutene  catalyzed  by 
H4SiWi204o  [29a]. 

Most  research  in  homogeneous  catalysis  by  heteropoly  compounds  involves 
oxidation  reactions  [28].  The  heteropolyanions  comprise  polyanions  where  the  addenda 
atoms  are  substituted  by  d-electron  transition  metals.  A  number  of  reactions  have  been 
studied,  mainly  using  monosubstituted  heteropolyanions  in  different  solvents.  There  are  a 
variety  of  mechanisms,  depending  on  the  reaction.  A  recent  review  [28]  describes  many 
of  the  established  mechanisms,  which  have  been  elucidated. 

Three  classes  of  catalysis  have  been  identified  [27]  for  heterogeneous  systems 
(i.e.,  gas-solid  and  liquid-solid).  For  the  first  type,  surface  catalysis,  reactions  take  place 
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on  the  two-dimensional  surface  (outer  surface  or  pore  walls)  of  solid.  The  reaction  rate  is 
proportional  to  the  catalyst  surface.  The  second  class,  bulk  type  I  catalysis,  the  reactants 
are  absorbed  into  the  solid,  between  the  polyanions  but  not  in  the  anion.  Usually,  these 
molecules  substitute  waters  of  crystallization  or  expand  the  lattice  where  reactions  take 
place.  The  structure  of  the  polyanion  is  intact.  The  reaction  rate  is  proportional  to  the 
volume  of  the  catalyst  in  the  ideal  case.  The  rate  of  an  acid  catalyzed  reaction,  for 
instance,  is  proportional  to  the  total  number  of  acidic  groups  in  the  solid  bulk.  The  third 
class,  bulk  type  II,  occurs  for  some  oxidation  reactions  at  high  temperatures.  The 
principal  reaction  probably  proceeds  on  the  surface,  the  whole  solid  bulk  takes  part  in 
redox  catalysis  owing  to  the  rapid  migration  of  proton  and  electron-redox  reactants  into 
the  bulk.  The  rate  is  proportional  to  the  volume  of  catalyst  in  the  ideal  case.  Although 
the  three  classes  are  distinctly  different  from  each  other  in  the  ideal  case,  the  extent  of 
contribution  of  the  inner  bulk  of  the  solid  to  the  catalysis  depends  on  the  rate  of  the 
catalytic  reaction  relative  to  the  diffusion  of  the  reactant  and  product  molecules  in  the 
bulk  type  I  catalysis  and  on  the  rate  of  diffusion  of  redox  carriers  for  the  bulk  type  II 
catalysis. 

Heteropoly  acids  have  been  proven  to  be  of  value  in  fundamental  studies  and 
practical  applications.  The  utility  of  these  compounds  is  a  consequence  of  the  potential 
for  design  at  the  atomic  and  molecular  level.  Control  of  acidic  and  redox  properties, 
multifunctionality,  bulk-type  and  pseudoliquid  behavior,  and  the  unique  basicity  of  the 
polyanion,  makes  these  compounds  very  attractive  to  study. 
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The  work  presented  in  this  dissertation  is  concerned  with  acidity  characterization 
of  12-tungstophosphoric  acid  and  its  cesium  derivatives.  The  heteropoly  acid  H3PW12O40 
offers  a  unique  opportunity  to  be  studied  in  homogeneous  and  heterogeneous  system  due 
to  a  variety  of  media  in  which  this  acid  shows  solubiHties.  Progressive  substitution  of  the 
protons  of  H3PW12O40  by  cesium,  forms  compounds  of  the  type  CSXH3.XPW12O40  which 
have  very  Hmited  solubility.  The  acidities  of  the  cesium  compounds  are  compared  to  the 
parent  acid  in  heterogeneous  reaction  because  of  their  increasing  importance  in  catalytic 
reactions.  Moreover,  characterization  of  H3PW12O40  supported  on  solids  such  as  silica, 
carbon,  and  alumina  were  carried  out.  Comparison  with  the  unsupported  acid  allows  for 
an  investigation  of  the  effects  of  surface  area  on  the  acidity  of  the  solid.  Mathematical 
models  are  used  to  interpret  acidity  of  the  acids  in  homogeneous  and  heterogeneous 
systems  in  attempt  to  obtain  quantitative  understanding.  Thermodynamic  data  for  these 
compounds  may  be  used  for  comparing  catalytic  activity  and  acidity  for  different  acid 
catalyzed  reactions  with  the  heteropoly  acids  studied. 


CHAPTER  2 

ACIDITY  OF  H3PW,2O40  IN  HOMOGENEOUS  SYSTEM 

Introduction 

The  utility  of  heteropoly  acids  (HPAs)  as  catalysts  for  redox  and  acid  catalyzed 
reactions  has  stimulated  considerable  interest  in  these  compounds.  Our  interest  in  acid 
catalysis  led  us  to  investigate  the  acidity  of  these  materials.  The  acid  H3PW12O40 
(abbreviated  as  H3PW)  was  selected  for  study  because  it  is  reported  to  be  among  the 
strongest  acids  in  the  series,  it  is  commercially  available  as  a  crystalline  hydrate  and  it  has 
been  thoroughly  characterized  [32a-c]. 

The  1 2-tungstophosphoric  acid  belongs  to  the  Keggin  class  of  heteropoly 
compounds.  The  complete  structure  of  H3PWi204o  nH20  can  be  depicted  by  a  primary 
and  a  secondary  structure.  The  primary  structure  (Keggin  polyanion)  has  twelve  WOe 
octahedra  surrounding  a  central  PO4  tetrahedron.  The  twelve  W0(,  octahedra  are  better 
described  as  four  groups  of  three  edge-shared  units  of  W3O13  which  are  comer-shared  to 
each  other.  These  triplet  groups  actually  do  not  exist  independently,  but  they  are  the 
smallest  block  needed  to  build  the  entire  Keggin  anion.  There  are  four  kinds  of  oxygen 
atoms  in  this  polyanion:  4  Oj  (internal  oxygen,  connecting  P  and  W),  12  Oe  (edge-sharing 
oxygen,  connecting  W  atoms),  12  Oc  (comer-sharing  oxygen,  connecting  W3O13  units), 
and  12  Ot  (terminal  oxygen,  bonded  to  one  W  atom).  The  W  atoms  are  displaced  outward 
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(0.3-0.4  A)  from  the  center  of  the  W0(,  octahedron,  giving  the  W-Ot  bond  some  double 
bond  character.  The  Oi  is  located  inside  the  anion  and  the  W-Oj-P  bond  is  long  and  weak. 
The  whole  Keggin  anion  has  a  tetrahedral  symmetry,  and  considering  the  anion  as  a 
sphere,  it  has  a  diameter  of  approximately  10  A. 

The  secondary  structure  is  dependent  on  the  hydration  degree,  for  example  in  the 
(H3PWi204o-6H20)  the  polyanions  constitute  two  sets  of  three-dimensional  networks  in 
which  polyanions  are  connected  by  H(H20)2^  bridges.  The  two  sets  of  networks 
penetrate  each  other,  so  that  the  arrangement  of  polyanions  becomes  the  bcc  type.  This  is 
the  most  dense  secondary  structure  and  is  called  type  A.  Types  B,  C,  and  D,  which 
contain  a  greater  number  of  water  molecules,  are  also  known.  In  the  dehydrated 
compound,  a  proton  located  in  the  same  position  as  H^(H20)2  interacts  with  the  four  Ot 
atoms.  Since  all  terminal  oxygens  are  equivalent,  according  to  the  magic  angle  '^O  NMR 
spectrum,  the  proton  is  proposed  to  exchange  rapidly  between  the  four  Ot  atoms  [29a]. 
The  primary  and  secondary  structure  of  H3PW  is  illustrated  in  Figure  2-1. 

In  solution,  the  structure  of  H3PW  depends  upon  its  concentration  and  the  solvent 
in  which  it  is  dissolved.  Large  angle  X-ray  scattering  (LAXS)  analysis  of  concentrated 

2  1 

H3P W  solutions  (>  1 0"  mol  L" )  yields  radial  distribution  curves  which  can  be  compared 
with  those  calculated  from  the  known  crystallographic  dimensions  of  the  presumed  solute 
species.  This  method  has  shovm  [33a-d]  that  the  structiu-e  of  polyoxometalates  are 
maintained  when  dissolved  in  water.  Other  evidence  comes  from  similar  UV  spectra  of 
solutions  of  isostructural  compounds  (e.g.,  H3PW12O40,  H4SiWi204o,  and  H6[H2Wi204o]), 
from  IR,  Raman,  and  '^O  NMR  results  of  aqueous  and  non-aqueous  solutions  [32b]. 
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(b) 


Figure  2-1.  Structure  of  H3PW12O40:  (a)  Primary,  Shows  the  Four  Types  of  Oxygen 

Atoms  in  the  Polyanion:  (A)  P-0;  (B)  W=0,;  (C)  W-O-Wedge;  (D)  W-O-Wcomer-  [51] 

(b)  Secondary,  Displays  H(H20)2"'  Bridging  Four  Polyanions  for  the  Trihydrated  Acid 
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HaPW  and  many  of  its  salts  are  fully  dissociated  in  aqueous  solution  as  a  result  of 
the  large  anion  size,  combined  with  a  low  surface  density  of  the  negative  charge.  The 
heteropolyanions  are  weakly  solvated  and  behave  as  weak  ligands  [34].  H3PW  is  a  strong 
Bronsted  acid.  In  solution  at  pH  higher  than  1.5,  the  following  equilibrium  reactions 
occur: 

PW,204o^- +  3H2O  PWn039^-  +  W04^-  +  6H*  12W04^-  +  HP04^-  +23H^ 
Although  the  Keggin  structure  is  preserved  in  concentrated  aqueous  solutions, 
degradative  dissociation  occurs  in  diluted  solution  [32b]. 

Acidities  of  heteropoly  acids  in  solution  have  been  studied  by  conductivity  [35- 
36],  Hammett  acidity  constants  [37],  and  the  '^C  shift  of  mesityl  oxide  [38].  Based  on 
conductivity,  it  is  claimed  that  in  solvents  of  moderate  polarity  like  ethanol,  the  first  and 
second  proton  of  H3PW  are  completely  dissociated  while  in  acetic  acid,  H3PW  behaves  as 
a  monoprotic  acid  [32b,  38].  The  pKa  for  the  dissociation  of  the  protons  in  acetone  are 
1 .6,  3.0,  and  4.0.  The  pKg  values  in  ethanol  are  very  similar  to  those  in  acetone.  In  water, 
the  stepwise  three  step  dissociation  of  H3PW  is  not  detected  because  of  the  leveling 
influence  of  the  solvent.  However,  in  acetic  acid,  pKai  is  equal  to  4.77  [32b]. 

The  most  commonly  reported  criterion  of  acid  strength  for  H3PW  has  been 
Hammett  acidity  fianctions  (Ho).  H3PW  is  pretreated  under  different  thermal  conditions, 
transferred  to  a  solution  of  dry  benzene,  where  small  amounts  of  n-butylamine  is  added. 
Then,  the  solution  is  tested  with  a  few  drops  of  0.1%  solution  of  an  appropriate  indicator 
after  equilibration.  Each  indicator  has  a  strength,  according  to  its  dissociation  constant  in 
benzene.  Then,  the  acidity  is  determined  from  the  number  of  millimoles  of  n-butylamine 
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which  had  to  be  added  before  the  appropriated  color  change  occurred.  In  this  fashion,  it 
has  been  determined  [37]  0.90  mmol  g"'  at  Ho  -5.8  (strongest  sites),  and  3.09  mmol  g"'  at 
Ho  6.8  (weakest  sites)  for  H3PW  pretreated  at  320  V  (best  condition).  It  has  also  been 
reported  [29a,  29b]  the  measurement  of  Hammet  acidity  in  H2O,  H2O-ACOH,  H2O- 
dioxane,  HaO-EtOH,  and  H20-Me2CO  mixtures.  Other  studies  [39a-c]  suggested  that 
H3PW  is  a  superacid,  based  on  the  definition  that  a  superacid  is  an  acid  with  strength 
greater  than  that  of  100%  H2SO4,  i.e.,  a  value  of  Ho  <  -12. 

The  '^C  shift  of  mesityl  oxide  upon  reaction  with  H3PW  in  acetic  acid  is 
interpreted  to  indicate  that  the  three  protons  dissociate  simultaneously  and  not 
consecutively  [38].  Because  of  the  size  of  the  molecule,  the  negative  charge  of  the 
conjugate  base  is  claimed  to  be  delocalized  over  the  large  Keggin  structure,  and  there  is 
no  charge  destabilization  for  the  second  and  third  ionizations.  In  all  instances,  the 
reported  studies  suggest  that  the  acidity  of  H3PW  is  greater  than  HCIO4  and  H2SO4  [29a, 
29b,  38]. 

In  the  solid  state,  H3PW  possess  Bronsted  acidity,  which  is  claimed  to  be  stronger 
than  such  conventional  solid  acids  as  Si02-Al203,  H3P04/Si02,  and  HX  and  HY  zeolites 
[29a].  A  microcalorimetric  study  [40]  of  the  reaction  of  gaseous  ammonia  with  solid 
H3PW  was  carried  out  at  50  °C.  The  sample  activated  at  150  °C  shows  different  behavior 
from  one  activated  at  250  °C,  which  differs  [41]  from  a  sample  activated  at  300  °C  and 
titrated  at  150  °C.  Another  microcalorimetric  study  [42]  at  150  °C  claims  that  H3PW, 
treated  at  150  C  or  250  C  under  vacuum,  produces  the  anhydrous  acid,  whose  acidity  is 
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greater  than  that  of  the  anhydride  prepared  by  vacuum  drying  at  400  °C.  The  anhydride  is 
metastable  decomposing  into  constituent  oxides  [42]. 
H3PWi204o-nH20  -»  H3PW12O40  +  nHjO    (Anhydrous  acid) 
nH3PW,204o  ^  [PWi2038.5]n  +  l.SnHzO    (Anhydride  phase) 
[PW,203g5]n  ->  O.SnPzOs  +  12nW03  (Oxides) 

The  exotherms  from  the  gas-soUd  calorimetric  titration  of  studies  of  H3PW  with 
ammonia  are  complex.  Initially  enthalpies  of  47.8  kcal  mof'  are  reported  dropping  to  45 
kcal  mof'  before  0.5  equivalents  of  ammonia  are  added.  This  is  followed  by  a  gradual 
drop  to  39  kcal  mol"'  up  to  2.9  equivalents  and  a  sharp  drop  thereafter.  The  three  protons 
are  reported  to  have  comparable  acidity.  In  the  Cal-ad  study  [43]  of  HZSM-5,  gas-solid 
calorimetry  at  elevated  temperature  was  shovra  to  contain  both  dispersion  and  acid-base 
components  and  to  average  into  a  single  step  different  protons  whose  acidities  differ  by 
30  kcal  mol"'. 

Thermal  methods  such  as  TGA,  DTA,  TPR,  and  TPD  have  been  used  to 
characterize  the  stability  and  nature  of  the  sites  of  H3PW  [32c,  44-47].  These  data 
suggest  that  the  Keggin  structure  is  stable  up  to  450-500  "C.  Water  fi-om  the  structure  is 
considered  zeolitic  in  nature,  since  its  release  does  not  affect  the  crystal  lattice.  In 
addition,  desorption  of  pyridine  has  shown  that  the  maximum  uptake  of  pyridine  is  about 
three  per  Keggin  anion,  and  occurs  when  the  H3PW  is  treated  at  190-320  °C. 

IR  spectroscopy  was  also  utilized  to  characterize  the  Keggin  structure  and  to 
monitor  reactions  of  basic  probes  with  H3PW  [48-51].  Table  2-1  shows  the  vibrational 
bands  assigned  to  the  primary  structure  PWi204o^'.  Variable  amounts  of  water  do  not 
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cause  shifts  in  these  peaks.  A  peak  about  1715  cm''  indicates  the  presence  of  hydronion 
(H(H20)n^),  which  reacts  with  probes  such  as  NH3  or  Pyridine. 


Table  2-1.  Fingerprint  Bands  of  the  Keggin  Anion  for  H3PW12O40. 


Bond 


Frequency  (cm  ) 


(P  — O) 


1081 


(W  — O.) 

(W  — Oc  — W) 


984 


892 


( W— Oe  — W) 


794 


(O— P— O) 


595 


The  states  of  the  protons  and  phosphate  group  in  the  solid  H3PW  have  been 
studied  by  solid  state  NMR  spectroscopy  of  ^'P  and  'H  [52-53].  Decreasing  the  amount 
of  water  in  the  structure  reduces  proton  mobility,  and  leads  to  stronger  bonding  to  the 
polyanion.  This  is  indicated  by  broadening  and  weakening  of  a  sharp  proton  resonance  of 
H3PW-17H20  (6  =  7.5  ppm)  to  H3PW  dehydrated  (5  =  9.2  ppm).  The  ^'P  NMR  chemical 
shift  is  very  dependent  on  the  number  of  water  molecules  on  H3PW-nH20,  the  values 
being  for  n=6  (5  =  -15.1  to  -15.6  ppm),  and  n=0  (6  =  -11.1  to  -10.5  ppm).  It  is  suggested 
[29a]  that  the  protons  can  move  among  four  equivalent  positions,  W=0~H"'-0=W, 
attached  directly  to  oxygen  atoms  of  the  Keggin  units,  while  H(H20)2"'  is  connected  to  the 
heteropolyanion  by  hydrogen  bonding  at  terminal  oxygens. 

The  conflicting  reports  of  the  acidities  of  H3PW  have  prompted  further  study. 
With  a  three  step  equilibrium  possible,  solution  measurements  are  more  easily  interpreted 
than  calorimetric  studies  of  the  solid.  Enthalpies  are  measured  as  the  criterion  of  acid 
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strength,  and  they  are  expected  to  provide  a  temperature  independent  relative  acidity  scale 
for  the  temperature  ranges  used  in  acid  catalysis.  Work  from  this  laboratory  has 
emphasized  [54-56]  using  the  enthalpy  of  adduct  formation  in  poorly  solvating  solvents 
as  the  criterion  for  measuring  donor-acceptor  strength.  Most  of  the  acceptors  studied  to 
date  are  neutral  molecules  that  form  neutral  1:1  adducts.  The  strong  acid  H3PW,  on 
reaction  with  moderately  strong  bases,  forms  ionic  products,  and  the  acid  is  insoluble  in 
poorly  solvating  solvents.  These  properties  afford  an  opportunity  to  extend  the  enthalpy 
criterion  of  acid  strength  to  new  acid-base  reactions  but  require  a  reevaluation  of  the 
experimental  conditions  that  will  provide  the  most  meaningfiil  results. 

Reaction  of  one  molecule  of  base  with  a  monoprotic  acid,  HX,  is  described  by  the 
following  equilibria  [57]: 

Bsoi  +  H-Xsoi         B-H-Xsoi  ^  B-H^-X-soi         BH^oi  +  X'soi  (2.1) 

(I)  (11)  (III) 

where  (I)  is  the  hydrogen  bond  adduct,  (II)  is  an  intimate  ion  pair  with  HX  behaving  as  a 
Bronsted  acid  and  (III)  is  the  dissociated  or  solvent  separated  ion  pair  of  the  Bronsted 
acid,  HX.  In  the  solvent  separated  ion  pair,  the  base  and  at  least  one  solvent  molecule  is 
coordinated  to  the  proton.  Research  leading  to  the  Coordination  Model  [58]  for  non- 
aqueous solvents  has  shown  that  the  position  of  the  equilibrium  between  species  (I),  (II) 
and  (III),  in  different  pure  bases  as  solvent,  depends  on  the  basicity  and  polarity  of  the 
solvent  as  well  as  the  basicity  of  the  conjugate  base  X'.  Both  strong  basicity  and  large 
polarity  are  needed  to  displace  the  anion  from  the  acid  by  a  neutral  base. 
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To  establish  an  acidity  scale,  the  solvent  should  be  held  constant  to  minimize 
differences  from  non-specific  solvation  contributions  to  the  measured  enthalpy.  The  best 
solvents  for  studying  Bronsted  acids  are  non-protonic  and  weakly  basic  because  of  the 
leveling  effect.  Basic  solvents  drive  the  equilibrium  in  Equation  (2.1)  toward  species  (III) 
even  when  HX  is  a  weak  acid,  i.e.,  those  with  strong  conjugate  bases.  As  a  result,  the 
enthalpy  of  reaction  of  a  base  with  both  weak  and  strong  acids  in  basic  solvents  will  be 
similar  and  will  correspond  to  the  enthalpy  of  the  reaction  (where  S  is  a  solvent 
molecule): 

SzH^oi  +  X'soi  +  Bsoi  ^  BHS^oi  +  X-soi  +  S  (2.2) 
Slight  differences  in  -AH  may  arise  if  the  equilibrium  between  intimate  and  separated  ion 
pairs  (species  II  and  III  of  Equation  2.1)  differ.  Thus,  weakly  basic  solvents  must  be  used 
to  discriminate  the  relative  acidity  of  strong  acids. 

The  best  solvents  to  study  acidity  should  also  be  either  non  polar  or  only  slightly 
polar,  so  the  product  of  the  reaction  is  species  (II),  or  very  polar,  so  the  only  product  is 
species  (III).  In  the  studies  reported  here  acetonitrile  is  selected  as  the  weakly  basic,  polar 
solvent  because  solubility  limitations  of  H3PW  require  use  of  a  polar  solvent.  When  a 
strong,  weakly  polar  base  such  as  pyridine  is  added  to  H3PW  in  a  weakly  basic  polar 
solvent,  solvent  polarity  will  allow  pyridine  basicity  to  drive  the  reaction  towards  species 
(III)  displacing  X"  from  the  proton.  The  enthalpy  measured  will  be  the  net  of  displacing 
X'  from  the  H"^  by  pyridine  plus  the  net  non-specific  solvation  of  the  ionic  products  and 
reactants.  The  non-specific  solvation  is  expected  to  be  the  smaller  of  the  two 
contributions  and  relatively  constant  as  the  base  is  varied.  Thus,  in  weakly  basic  polar 
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solvents,  relative  basicities  can  be  obtained  from  the  enthalpies  without  varying 
contributions  from  the  extent  of  ionizations  of  (II)  to  (III).  When  di-  and  tri-protic  acids 
are  studied,  the  above  discussion  for  HX  pertains  to  each  step  of  the  reaction. 

With  H3PW,  a  series  of  reactions  vAth  different  bases  have  been  investigated  to 
determine  the  influence  of  the  base  on  the  three  steps  of  the  H3PW  equilibria.  The 
enthalpies  for  step  1  are  correlated  to  the  ECW  model  [56]  to  demonstrate  that  the 
basicity  trends  in  acetonitrile  toward  Bronsted  acids  parallel  those  toward  neufral  acids  in 
weakly  solvating  solvents.  In  order  to  compare  the  relative  acid  strength  of  H3PW  to 
other  acids  using  this  new  scale,  enthalpies  for  the  reaction  of  pyridine  with  a  variety  of 
acids  are  also  reported  in  CH3CN  solvent. 

Experimental 

Materials 

Elemental  and  TG  analysis  of  the  H3PW  (obtained  from  Aldrich)  revealed  16 
moles  of  water  per  mole  of  H3PW.  Acetonitrile  obtained  from  Fisher  was  purified  by 
drying  over  4A  molecular  sieves  for  24  hours,  and  then  distilled  over  P2O5.  Pyridine 
obtained  from  Fisher  was  distilled  over  CaHj.  Quinuclidine,  N-methylimidazole,  N,N- 
dimethylacetamide  and  triethylphosphine  (Aldrich),  quinoline  and  dimethylsulfoxide 
(Fisher)  were  used  as  supplied.  Sulftiric  acid  (97%  as  calculated  by  standardization),  p- 
toluenesulfonic  acid,  and  benzoic  acid  were  obtained  from  Fisher.  Triflic  acid  was 
obtained  from  Aldrich.   Trifluoroacetic  acid  and  m-chlorobenzoic  acid  were  obtained 
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from  Kodak.  Drying  and  standardization  procedures  for  the  acids  followed  literature 
guidelines.  Dried  samples  were  stored  in  a  container  with  4A  molecular  sieves. 

Thermal  Treatment  and  Calorimetric  Titration  of  H^PW 

H3PW  was  subjected  to  different  thermal  treatments.  It  was  heated  at  128  °C  and 
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162  C  in  a  vacuum  oven.  In  other  treatments,  it  was  heated  at  200  °C,  300  °C,  and  400  °C 
under  dry  air  flow  in  a  tube  furnace.  After  each  treatment,  the  H3PW  is  transported 
immediately  to  a  dry  box  and  stored  under  nitrogen. 

Samples  of  H3PW  were  weighed  and  transferred  to  a  vacuum  insulated 
calorimetric  cell,  containing  a  stir  bar.  Since  H3PW  is  very  hygroscopic,  these  weightings 
and  transfers  were  done  in  a  dry  box.  For  each  titration,  100  ml  of  acetonitrile  is  added  to 
the  cell.  A  calibrated  syringe,  filled  with  a  solution  of  known  concentration  of  base  (e.g., 
0.2  M),  is  inserted  into  the  cell  along  with  a  thermistor  and  a  heater  coil.  These 
operations  were  carried  out  in  an  inert  atmosphere  glove  bag.  The  thermistor  and  heater 
coil  were  cormected  to  an  electronic  circuit  bridge,  and  a  computer  [59].  The  thermistor 
was  calibrated  with  the  heater  coil  prior  or  immediately  after  each  titration.  A  set  of  12 
calibrated  brass  stops  were  used  for  the  syringe  injections.  After  each  injection  the  heat 
evolved  from  reaction  of  the  base  and  H3PW  is  measured.  The  heat  from  adding  the  base 
solution  to  the  solvent  is  negligible.  Each  titration  was  repeated  three  times  for  each  base 
studied.  Titrations  of  acids  other  than  H3PW  with  pyridine  were  performed  twice. 
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Spectral  Analysis 

A  Perkin-Elmer  Lambda  6  UV-vis  spectrophotometer  and  suprasil  quartz  cells  of 
0.1  cm  pathlength  were  employed.  A  Nicolet  Model  5PC  FTIR  spectrophotometer  and  a 
NaCl  cell  of  0.1  mm  pathlength  or  mull  in  a  NaCl  cell  were  employed. 

Thermal  Analysis 

The  DSC  and  TG  analyses  were  carried  out  in  a  TA  DuPont  2000  instrument 
(module  2910  DSC  and  module  2950  TGA)  at  rate  of  10  °C  min"'.  Isothermal  TGA  s 
were  run  for  2  hours.  The  samples  of  H3PW  were  taken  directly  from  the  bottle. 


Calculations 

The  heat  liberated  in  the  reaction  of  strong  bases  with  H3PW  in  acetonitrile  can  be 
calculated  according  to  the  following  set  of  equilibria,  giving  the  respectiye 
thermodynamic  parameters: 

H3PW  +  B         H2PW  +  BH^;        Ki,    AHi  (2.3) 
H2PW  +  B  —  HPW^-  +  BH^;       K2,    AH2  (2.4) 
HPW^-  +  B         PW^'  +  BlT;         K3,    AH3  (2.5) 
For  ease  of  representation,  the  acetonitriles  attached  to  the  protons  of  the  H3PW  species 
are  not  shown.  The  equilibrium  constants  are  represented  by  the  expressions: 

[H,PW][B]  '      IH,PW-][B]    ^    '      '     IHPW''][B]  ^  ' 
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where  the  brackets  indicate  concentration  of  the  species  (mol  L'').  The  ionic  strength  for 
all  solutions  in  this  work  was  constant  (1.39x10"^  mol  L"').  Therefore,  the  activity 
coefficients  are  kept  approximately  constants  and  the  concentration  equilibrium  constants 
calculated  here  are  close  to  the  true  thermodynamic  equilibrium  constants.  Transforming 
stepwise  equilibrium  constant  expressions  (Equations  2.6  to  2.8)  to  global  formation 
constants,  the  following  equations  are  obtained: 
_[H2PW-][BH'] 

'       [H,PW][B]  (^-^^ 

^^""^•^^^    [H,PW][Bf  (2.10) 
n  -KK  K  APW'-][BW]' 

P.-K^K,K,-  ^^^^^^^^^3  (2.11) 

The  mass  balance  for  this  system  is  described  by  the  total  concentration  (mol  L"')  of  acid 
([H3PW]t)  and  base  ([B]t): 

[HjPWJt  =  [H3PW]  +  [H2PW]  +  [HPW^-]  +  [PW^-]  (2. 12) 

[B]t  =  [B]  +  [BH"]  (2.13) 
Substituting  the  expressions  for  concentrations  of  HaPW,  HPW^",  and  PW^"  (based  on 
Equations  2.9  -  2.1 1)  in  Equation  (2.12),  yields: 

(2.14) 

[H,PW],  =  [H,PW]il  +  +  J3,  +  /?,  -M_ )  .2  15) 

"^'[BH^]    ^'[BH*f    ^'[BH'f'  ^ 
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"^'[5//^]    ^'[5//^]^    ^'[BH^f^  ^  ' 

The  total  calculated  heat  evolved  (h'cai)  from  that  system,  when  acid  reacts  with  base,  is 
given  by: 

heal  =  V{[H^PW-]M{,^  +[HPW'-]AH/  +[PW'-]AH/)  (2.17) 
where,  V  is  the  volume  of  the  titration  solution  (L),  and  AHi^  is  the  global  enthalpy  of 
deprotonation  of  H3PW  i.e.,  AH,^  =  AH,,  AHj^  =  AH.+AHj,  and  AHs"^  =  AH1+AH2+AH3, 
with  AHi  being  the  stepwise  formation  enthalpy.  Substituting  the  expressions  for 
concentrations  of  H2PW',  HPW^",  and  PW^'  (based  on  Equations  2.9-2.11)  in  Equation 
(2.17),  the  following  equations  are  obtained: 

[H,PW]  =  h'ca,/  V{p,  A///  +  /?2  A//, '  +  A//3 (2.1 9) 

Equations  (2.16)  and  (2.19)  show  the  relation  of  free  concentration  of  acid  with  the 
equilibrium  constants  and  enthalpies  of  deprotonation.  Setting  Equations  (2.16)  to  (2.19) 
equal  to  each  other,  and  rearranging  them,  Equation  (2.20)  is  obtained: 

[H,PW],V(P,  A///        A^^  T^^  A^^  r^ 

h\.  =  —   ^'[^//-]  ^'[BHn'      '     ^UbH^'  .2  20) 

(l.^,J^,^^J^,^^J^) 
^'[5//^]    "^'[BH'f  ^'[BH'f^ 

Considering  R  =  [B]/[BH^]  as  the  ratio  of  free  base  and  protonated  base  concentrations 
(mol  L"'),  the  simplified  expression  of  the  calculated  heat  is  given  by  Equation  (2.21): 
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A  computer  program  was  written  to  solve  these  equilibria  in  non  aqueous  solution 
(Appendix).  Kj,  K2,  K3,  the  total  concentration  of  H3PW,  and  the  total  base 
concentration  determine  the  concentration  of  all  the  species  in  solution  using  Equations 
(2.6  -  2.8).  The  ratio  between  free  base  and  protonated  base  is  obtained  assuming  an 
initial  value,  which  allows  the  calculation  of  concentrations  for  all  species.  The  routine 
adjusts  the  concentrations  through  a  mass  balance  of  the  protons  i.e.,  the  protons  on  BH* 
must  correspond  to  the  ones  lost  from  the  initial  acid.  The  calculation  is  repeated  until 
the  values  are  consistent.  Knowing  AHi,  AH2,  and  AH3  as  well,  Equation  (2.21)  can  be 
solved  for  h  cai-  Actually,  the  values  of  the  K  s  and  AH  s  are  unknown.  Using  initial 
estimates  of  these  parameters,  an  iterative  program  is  used  to  calculate  the  final  values 
which  give  the  minimum  of  the  non-weighted  sum  of  squares  of  the  residuals  (U)  of 
Equation  (2.22)  by  a  non-linear  routine. 

U(Pi,AHi^)  =  Z(hexp-h'cai)'  (2.22) 
The  method  of  steepest  descent  is  employed  to  find  the  best  values  of  Pi  and  AHj.  The 
first  derivatives,  5U/5pi  and  6U/5Hi  are  approximated  as  AU/Apj  and  AU/AHi  which  are 
easily  calculated  to  give  h  cai  (calculated  heat).  Using  h'cai  and  h  exp  (experimental  heat), 
the  minimum  of  the  fimction  (Equation  2.22)  is  obtained  (e.g.,  a  sample  data  fit  is  shown 
in  Table  2-3). 

In  order  to  evaluate  the  adequacy  of  the  proposed  model,  the  experimental  data 
were  repeated.  The  residuals  (h'exp  -  h'cai)  were  randomly  distributed  about  the  regression. 
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A  comparison  between  the  dispersion  (variance)  from  the  regression  and  one  obtained  by 
repHcation  is  made.  To  measure  the  dispersion  of  the  experimental  points  (the  external 
vanance,  a  ext),  the  sum  of  squares  of  the  residuals  is  obtained  and  divided  by  the  degrees 
of  freedom  (number  of  measured  points  minus  the  number  of  fit  parameters).  An  internal 
variance  (ct  int)  is  calculated  in  a  similar  manner,  but  the  sum  of  squares  is  based  on  the 
average  of  the  replicates  and  the  summation  is  done  point  to  point  without  adding  the 
heats.  Since  the  heats  are  added  up  for  analysis  purposes,  an  error  associated  with  each 
measurement  is  automatically  inserted  in  each  point.  The  system  HsPW-Pyridine,  for 
instance,  shows  a\xt=l -4903x10"^  and  a^nt=  1.1 405x10"^.  Assuming  that  a^im  is  constant 
for  these  calorimetric  measurements,  and  the  errors  are  normally  distributed,  both 
variances  can  be  compared  by  means  of  an  F  test.  The  ratio  a^ext/cT^im  gives  1 .30.  At  10% 
probability  for  the  F  distribution  (F=  3.05,  ^2  =  h""  6)  both  measures  of  dispersion  are 
equivalent.  This  test  of  the  goodness  of  fit  [60],  although  not  completely  rigorous  in 
nonlinear  modeling,  has  at  least  semiquantitative  value,  and  can  tell  if  the  model  provides 
a  good  to  fit  the  data.  Applying  the  F  test  to  all  systems  gave  similar  results  to  the 
example  shown. 

Another  method  to  evaluate  the  standard  deviations  for  the  K  s  and  AH  s  uses  a 
variance-covariance  matrix  of  parameters  [60],  which  gives  errors  as  standard  deviations. 
This  analysis  provides  standard  deviations  that  match  those  obtained  by  replicated  data 
within  ±5%  of  the  reported  value. 

For  all  the  data  in  Tables  2-2,  2-4  and  2-5,  the  values  of  the  parameters  are  based 
on  averaged  data  sets  of  duplicate  or  triplicate  experiments  and  the  errors  are  based  on  the 
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standard  deviations  of  the  variance-covariance  matrix.  Since  precision  exceeds  accuracy, 
interpretations  of  differences  less  than  ±0.5  kcal  mol"'  will  not  be  made  for  -AHi  and  less 
than  1-2  kcal  mol''  for  AH2  and  AH^. 

Results  and  Discussion 

Preparing  Anhydrous  H^PW 

As  with  any  polyprotic  oxyacid,  drying  is  complicated  by  the  conversion  of 
dihydroxy  functionality  that  gives  rise  to  the  Bronsted  acidity  to  an  anhydride  with  |x-oxo 
fiinctionality.  For  heteropoly  acids,  the  anhydride  has  been  shown  by  gas-solid 
calorimetry  to  be  a  much  weaker  acid  toward  ammonia  than  the  anhydrous  acid.  The 
hydrated  acid  is  expected  to  be  a  weaker  acid  than  the  anhydrous  acid  because  of  the 
endothermic  displacement  of  coordinated  water.  Thus  it  is  critical  for  acidity  and 
catalytic  studies  to  use  the  anhydrous  acid.  In  view  of  the  conflicting  reports  in  the 
literature  [61-63]  on  the  influence  of  pretreatment  on  acidity  and  reactivity,  we  decided  to 
determine  the  optimal  drying  procedure  to  produce  the  strongest  acid.  The  total  heat 
released,  equilibrium  constants  and  enthalpies  from  calorimetric  titrations  with  pyridine 
are  given  in  Table  2-2  for  samples  prepared  by  different  thermal  treatments. 

The  samples  in  Table  2-2  were  analyzed  by  UV-vis,  FTIR,  XRD,  DSC  and  TGA. 
No  differences  in  the  UV-vis  [64]  or  FTIR  [48-51]  spectra  of  acetonitrile  solutions  were 
observed  confirming  that  the  Keggin  structure  is  intact  in  all  samples.  Powder  XRD 
results  are  identical  except  for  the  sample  calcined  at  400  °C  where  a  slight  broadening  of 
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Table  2-2.  Results  of  Calorimetric  Titrations  (25±1°C)  with  Pyridine  for  Different  H3PW 
Samples  in  CH3CN  Solvent. 


Sample* 

(°C/hs) 

Total 

Heat  Evolved 
(cal) 

-AH," 
(i=l  2  2) 
fkcal  mol  '^ 

128/4 

6.187 

19.8 

1.6x10' 

15.3 

1.1x10^ 

Z.'tXlU 

128/8 

6.365 

20.5 

3.6x10^ 

9.4 

6.5x10' 

20  0 

s  9vi n' 

162/1 

5.548 

18.5 

4.6x10^ 

9.2 

1.2x10^ 

15  3 

4  9x1 fi' 

162/2 

5.768 

19.9 

2.2x10' 

9.3 

1.2x10' 

16  7 

9  9x10* 

162/4 

6.647 

21.0 

2.2x10' 

11.8 

2.4x10^ 

1 8  f> 
1  o.u 

A  9v 1 n' 

t.ZX  1  u 

200/2 

5.511 

18.6 

2.2x10' 

10.1 

5.8x10^ 

1  7  S 

Q  9 

300/2 

4.338 

4.5x10'* 

11.1 

3.2x10' 

3.0 

1.5x10^ 

300/  16 

3.827 

11.4 

4.8x10^ 

7.7 

8.5x10^ 

9.0 

5.5x10' 

400/2 

1.307 

8.2 

1.5x10^ 

0.3 

6.4x10^ 

3.3 

3.2x10"' 
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Table  2-2.  continued 


a.  Samples  of  H3PW  treated  at  different  Temperature  (°C)  /  Time  (hs).  Samples  at  128 
and  162°C  were  evacuated  and  those  at  200,  300,  and  400°C  were  calcined  under  air  flow. 
Two  samples  of  H3PW-16H20  (as  purchased)  were  titrated  with  pyridine  solution 
providing  -12.9  and  -5.4  kcal  mol"'  for  AHi  and  AH2  with  errors  of  3%  and  6% 
respectively,  and  6.8x10^  and  1.6x10^  for  Ki  and  K2  with  errors  of  10%  and  27% 
respectively. 

b.  Based  on  variance-covariance  matrix  analysis,  standard  deviation  are  1%,  2%)  and  5% 
of  AHi,  AH2,  and  AH3,  and  5%,  20%,  and  25%  of  Kj,  K2  and  K3  respectively. 
Conclusions  are  based  on  a  error  estimate  of  ±0.5  kcal  mof'  for  AHi  and  ±1-2  kcal  mol"' 
for  AH2,  and  AH3. 


'  
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the  lines  and  loss  of  intensity  occurs.  Similar  XRD  results  [46]  are  reported  for  a  sample 
of  H3PW  calcined  at  400  C. 

DSC  and  TG  analyses  give  two  endothermic  peaks,  corresponding  to  loss  of  water 
at  92  and  215  °C,  and  a  third  peak  at  about  550  °C  (Figure  2-2).  After  the  initial  weight 
loss,  isothermal  TGA's  at  150,  200,  300,  and  400  C  for  2  to  4  hours  do  not  show  any 
flirther  change  (Figure  2-3).  A  DTA  study  [44]  of  H3PW  reports  decomposition  at  about 
590  °C  forming  WO3.  TGA  [46]  and  TPD  [45]  studies  claim  that  decomposition  of 
H3PW  results  from  mass  loss  of  phosphorous,  as  P2O5,  occurring  at  500  °C  depending  on 
the  duration  of  heating.  Traces  of  green  and  yellow  powder  are  seen  on  the  surface  of 
H3PW  for  300/2,  300/16,  and  400/2,  (°C/hr),  samples.  The  color  changes  suggests 
formation  of  anhydrous  WO3  (yellow)  and  reduced  Wn03n.i  (blue)  which  accounts  for  the 
green  color.  When  WO3  is  heated  in  vacuum,  reduction  occurs  forming  a  variety  of 
tungsten  blue  materials  [1]. 

The  resuhs  of  Table  2-2  clearly  show  that  the  sample  treated  at  162°C  for  4  hours 
(162/4)  under  vacuum  is  the  strongest  acid,  which  should  be  the  anhydrous  acid. 
Reported  gas  phase  calorimetry  and  TGA  results  also  indicate  that  this  treatment  will 
produce  the  anhydrous  acid  [42].  The  FTIR  of  this  162/4  sample  compared  to  the 
purchased  sample  of  H3PW  I6H2O,  Figure  2-4,  shows  the  absence  of  the  characteristic 
broad  absorption  band  for  water  around  3500  cm"'. 

In  order  to  substantiate  that  162/4  produces  the  anhydrous  acid  with  the  data  in 
Table  2-2,  a  calorimetric  titration  of  162/4  with  water  as  base  in  CH3CN  was  carried  out 
giving  a  -AHi  of  4.2  kcal  mol"'  and  Ki  =  1.0  xlO^  The  large  K  indicates  that  any  excess 


Figure  2-2.  Thermal  Analysis  of  H3PW12O40  by:  (a)  DSC  in  air/He  at  10  °C  min 
TGA  in  air  at  10  °C  min''. 


Figure  2-3.  Isothermal  TGA  (300  °C)  of  H3PW12O40  in  air  at  10  °C  min"'. 
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Figure  2-4.  FTIR  of  H3PW  in  Mineral  Oil.  (— )  H3PW  Treated  at  162  °C  for  4  hs  under 
Vacuum.  (— )  H3PW  Taken  from  the  Bottle  (I6H2O).  Absorption  Bands  of  Mineral  Oil 
(Aldrich,  dried  on  3  A  Molecular  Sieves)  are  on  the  Typical  Range  of  2800-3000  cm''  and 
1400-1600  cm"'.  The  Preparation  of  the  Mull  was  Made  in  a  Dry  Box  Filled  with 
Nitrogen. 


42 

water  in  the  samples  of  Table  2-2  would  be  coordinated  to  H3PW.  Since  the  -AHi  values 
for  the  vacuum  dried  samples,  128/4,  162/1  and  162/2  (°C/hr)  are  smaller  than  162/4  by 
an  amount  less  than  the  4  kcal  mol''  found  for  the  monohydrate,  we  can  conclude  that 
128/4,  162/1,  and  162/2  are  partially  wet.  The  sample  128/8  is  almost  as  dry  as  162/4  but 
the  total  heat  evolved  indicates  a  very  small  amount  of  water  remains.  Flowing  air  is  not 
as  efficient  as  vacuum  for  drying  as  shown  by  the  200/2  versus  the  162/2  (°C/hr)  result. 
Since  the  flowing  air  samples  300/2,  300/16  and  400/2  (°C/hr)  are  expected  to  have  less 
water  than  200/2,  the  lower  enthalpies  suggest  some  degree  of  reduction  and  anhydride 
formation  consistent  with  the  observed  color  changes  on  the  surface  of  the  solid.  A  large 
decrease  in  the  total  heat  evolved  for  titration  of  sample  400/2  is  consistent  with  the 
decreased  acidity  reported  in  the  gas-solid  calorimetric  results  [42]  and  illustrates  the 
insensitivity  of  the  XRD  pattern  to  the  changes  occurring  on  drying.  The  calorimetric 
results  in  Table  2-2  would  mitigate  against  use  of  H3PW  as  an  acid  catalyst  at  or  above 
300°C. 

The  hydrated  H3PW  I6H2O  was  also  titrated  with  pyridine  giving  a  -AHi  of  12.9 
kcal  mol''.  Subtracting  -AHi  for  the  reaction  of  H3PW  with  water  (  4.2  kcal  mol"')  from 
21.0  kcal  mol"'  for  the  pyridine  reaction  with  anhydrous  H3PW  produces  -AHi  =16.8  for 
the  reaction  of  the  monohydrate  with  pyridine.  Solution  species  of  higher  water 
molecularity  in  H3PW  I6H2O  than  in  the  monohydrate  lead  to  a  more  endothermic  water 
displacement.  Other  complications  from  the  presence  of  water  are  evident.  With  all 
bases  used  in  this  study,  properly  dried  H3PW  leads  to  the  enthalpy  order  of  -AHi  >  -AH2 
<  -AH3  (vide  Table  2-4).  In  some  instances,  this  is  not  observed  in  Table  2-2.  Since  the 
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reactions  involved  in  steps  2  and  3  for  wet  samples  are  not  known,  the  possible  incorrect 
formulation  of  the  equilibria  precludes  a  more  detailed  analysis  of  the  results  in  Table  2- 
2.  These  experiments  accomplish  their  main  objective  by  showing  that  the  treatment 
162/4  produces  a  sample  that  evolves  the  largest  total  heat,  gives  the  highest  -AHi 
comparable  to  that  of  CF3SO3H  (vide  infra),  is  driest  and  most  likely  is  the  anhydrous 
acid.  Unless  noted  otherwise,  all  of  the  results  reported  in  the  next  sections  of  this  study 
were  obtained  on  H3PW  samples  dried  under  vacuum  at  162  °C  for  4  hours. 

Solution  Species 

The  reaction  of  H3PW  in  CH3CN  with  various  bases  was  monitored  by  UV-vis 
spectra.  The  absorption  maximum  for  H3PW  in  CH3CN  occurs  at  A,  266  nm.  The 
spectra  are  essentially  unchanged  when  excess  of  most  of  the  bases  over  that  required  for 
neutralization  of  the  three  protons  is  added  to  H3PW  in  acetonitrile  (e.g..  Figure  2-5). 
Quinuclidine  is  an  exception  to  this  behavior.  With  this  base,  the  spectra  are  essentially 
unchanged  up  to  a  1:1  ratio  of  base  to  H3PW.  A  slight  decrease  in  absorbance  but  no 
band  shift  is  observed  up  to  a  2:1  ratio.  At  the  ratios  of  quinuclidine  to  H3PW  greater 
than  2: 1 ,  a  band  shift  from  266  nm  to  25 1  nm  is  observed  leading  to  an  isosbestic  point  at 
255  nm  (Figure  2-6).  Different  specfral  behavior  of  H3PW  in  different  organic  solvents 
has  been  noted  [65-68].  'H,  '^O,  ^'P,  and  '^^W  NMR  and  electronic  absorption 
spectroscopy  [67],  of  complexes  of  [DP2H"^]3[PW^"]  (where  DP  =  N-methylpyrrolidinone 
(NMP),  1,1,3,3,-tetramethylurea  (TMU),  and  1 ,3-dimethyl-2-imidazolidinone  (DMEU)) 
attribute  the  changes  to  reversible  intermolecular  electron  donor-acceptor  complexes  in 


44 


Figure  2-5.  Spectra  of  Absorption  H3PW  with  N-methylimidazole  in  CH3CN  (b  =  0.1 
cm).  The  Concentration  of  H3PW  is  1.05  xlO'"*  mol  L"'  (the  Same  for  All  Solutions). 
From  Bottom  to  Top  {X  =  230  nm),  Concentration  of  N-methylimidazole  in  the  Solutions 
are:  0;  2.17,  4.34,  6.51,  and  8.68  xlO"^  mol  L"'  Respectively. 


! 
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Figure  2-6.  Spectra  of  Absorption  H3PW  with  Quinuclidine  in  CH3CN  (b  =  0.1  cm).  The 
Concentration  of  H3PW  is  1.57  xlO"*  mol  L"'  (for  All  Solutions).  From  Top  to  Bottom  (k 
=  230  nm),  Concentration  of  Quinuclidine  in  the  Solutions  are  Respectively:  1.67;  3.37, 
5.05,  6.75,  8.37,  10.12,  and  1 1.82  xlO"^  mol  L"'. 
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solution.  The  cause  of  the  change  with  quinuclidine  was  not  pursued  in  view  of 
predominant  interest  in  the  1:1  species  in  this  work,  vide  infra. 

Calorimetry 

The  measured  amounts  of  heat  (h  exp)  evolved  from  the  calorimetric  titration  of 
H3PW  with  pyridine  in  acetonitrile  are  shown  in  Table  2-3  along  with  the  calculated 
values,  h  cai,  from  the  data  fit  to  a  three  step  equilibrium.  Attempts  to  fit  the  data  to  a  one 
step  equilibrium  (three  protons  reacting  with  the  same  strength)  or  to  a  two  step 
equilibrium  (one  different  and  two  equivalent)  lead  to  larger  sums  of  the  squares  of  the 
deviations  and  meaningless  parameters. 

The  values  of  Ki,  K2  and  K3  for  the  pyridine  reaction  are  2100,  240,  and  41 
respectively.  The  corresponding  enthalpies,- AH i,-AH2,  and-AHs,  are  21.0,  11.8  and  18.6 
kcal  mol"'.  The  thermodynamic  data  in  Table  2-4  reveal  a  very  interesting  pattern  and 
each  step  will  be  interpreted  separately. 

The  first  step  of  the  equilibria  existing  in  solutions  of  H3PW  in  CH3CN  is  given 
by  Equation  (2.23): 

CH3CN  +  (CH3CN)3H3PW  (CH3CN)2H^  +  (CH3CN)2H2PW  (2.23) 
To  simplify  the  discussion,  coordination  number  two  is  assumed  for  the  proton,  anionic 
H2PW'  species  will  be  assumed  to  have  two  acetonitriles  coordinated,  and  no  attempt  will 
be  made  to  distinguish  intimate  ion-pairing  from  hydrogen  bonding,  i.e.,  species  I  and  II 
of  Equation  (2.1).  The  position  of  the  equilibrium  in  Equation  (2.23)  is  unknown  but  the 
pKa  of  perchloric  acid  in  acetonitrile  [69]  is  2  and  the  pKai  for  H3PW  in  acetone  is  1.6. 
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Table  2-3.  Heat  Evolved  (h'exp)"  upon  Addition  of  Pyridine  to  H3PW  in  CH3CN  Solvent 
and  Calculated  (h  cai)  from  Equilibrium  Analysis.  [HsPWJj  =  1.39  x  10"^M. 


Base  Added  X  10''  (M)       h'exp  (cal)     h'cai  (cal) 


2.43 

0.499 

0.504 

4.54 

0.932 

0.936 

6.74 

1.377 

1.372 

9.05 

1.827 

1.813 

13.0 

2.624 

2.621 

18.1 

3.325 

3.341 

23.5 

4.146 

4.147 

27.1 

4.690 

4.683 

31.7 

5.374 

5.372 

36.1 

6.011 

6.012 

40.6 

6.458 

6.459 

45.2 

6.647 

6.648 

a.  The  base  added  and  h  exp  values  are  summed  giving  the  total  base  and  the  total  heat 
after  each  addition  i.e.,  4.54x10"^  M  and  0.932  calories  are  the  sums  for  the  first  two 
additions. 
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Table  2-4.  Thermodynamic  Parameters"  for  HaPW-Base''  (25±1°C)  in  CH3CN  Solvent. 


Qi  Py  M-lm  Tep  Qu'  H2O  DMA  THF 

Ki  4.8x10^  2.1x10^  1.5x10^  S.SxlO'*  1.1x10^  1.0x10^  6.1x10^  9.2xl0' 

(±9.5)  (±4.0)  (±1.8x10^)  (±7.4x10^)  (±1.4x10^)  (±0.5)  (±5.5xlO')  (±1.2) 

K2  2.6x10^  2.4x10^  l.lxlO''  4.5x10^  9.2x10^  3.8x10^  2.4x10^  5.6xl0' 

(±3.9x10')  (±3.4)  (±4.8x10^)  (±1.7x10^)  (±5.2x10^)  (±6.8xl0')  (±3.1x10^)  (±1.7x10*) 


K3  5.2x10' 

4.1x10' 

6.5x10' 

7.5x10^ 

1.3x10^ 

7.1 

4.7 

(±1.3x10') 

(±5.8) 

(±7.1) 

(±9.9x10') 

(±4.6x10') 

(±0.7) 

(±2.6) 

-AHi  18.1 

21.0 

29.5 

32.2 

42.6 

4.2 

13.2 

11.4 

(±0.4) 

(±0.2) 

(±0.4) 

(±0.7) 

(±2.1) 

(±1.0) 

(±2.0) 

(±0.5) 

Atr  7  1 
-Ari2 

118 
11.0 

10.  J 

1  J.J 

1 0  A 

o.y 

3.4 

0.3 

(±0.2) 

(±0.3) 

(±0.5) 

(±0.3) 

(±0.4) 

(±0.5) 

(±0.5) 

(±0.1) 

-AH3  15.9 

18.6 

23.4 

18.0 

26.0 

2.1 

12.0 

(±0.7) 

(±0.5) 

(±0.8) 

(±0.8) 

(±1.1) 

(±2.6) 

(+1.9) 

-AGi  3.7 

4.5 

5.7 

6.2 

8.3 

4.1 

9.3 

2.7 

-AG2  3.3 

3.2 

5.6 

5.0 

8.1 

3.5 

8.7 

2.4 

-AGs  2.3 

2.2 

2.5 

3.9 

2.9 

1.2 

0.9 

-ASi  48.3 

55.1 

79.8 

87.1 

115.0 

0.3 

13.3 

29.2 

-AS2  12.9 

28.8 

36.0 

34.7 

35.1 

-8.8 

-18.0 

-7.0 

-AS3  45.4 

55.1 

70.0 

47.3 

77.4 

3.0 

37.2 
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Table  2-4.  continued 

a.  The  units  for  AHj  and  AGi,  and  are  kcal  mol"'.  The  unit  of  ASj  is  cal  Kelvin"'  mole''. 

b.  Bases  are  listed  as  Quinoline  (Qi),  Pyridine  (Py),  N-methylimidazole  (M-im), 
Triethylphosphine  (Tep),  Quinuclidine  (Qu),  N,N-Dimethylacetamide  (DMA),  and 
Tetrahydrofuran  (THF).  The  errors  are  standard  deviations  calculated  by  variance- 
covariance  matrix  analysis.  Deviations  for  AH2  and  AH3  reflect  cumulative  errors  on  the 
measured  heat  values  (h  exp)-  Since  the  equilibrium  constants  are  calculated  from  h  exp 
values,  small  errors  on  these  measurements  cause  large  errors  on  equilibrium  constants. 
Our  main  criterion  for  acidity  is  the  enthalpy  and  these  are  very  well  determined. 

c.  UV-vis  spectral  changes  indicate  this  data  may  reflect  a  more  complex  reaction. 
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Slight  ionization,  Equation  (2.23),  is  expected  in  acetonitrile  from  these  comparisons  and 
-AHi  for  the  reaction  with  pyridine  corresponds  predominantly  to  displacement  of 
(CH3CN)2H2PW"  to  form  CsHsNHCHaCN^  (Equation  2.24,  more  accurate  representation 
of  Equation  2.3).  Both  basicity  and  solvent  polarity  are  needed  to  dissociate  the  anion  or 
form  a  solvent  separated  ion-pair  [58]. 

(CH3CN)3H3PW  +  B         BHCHsCN^  +  (CH3CN)2H2PW  (2.24) 
It  must  be  take  into  account  that  Equation  (2.23)  as  well  as  the  second  and  the  third 
proton  dissociation  steps  of  H3PW  should  be  combined  with  the  expression  for  the 
protonation  of  the  base  (B  +  (CH3CN)2H*  BHCHsCN^  +  CH3CN)  in  order  to 

deduce  the  expressions  for  the  equilibrium  constant  defined  in  Equations  (2.3)  to  (2.5). 

The  second  and  third  ionization  steps  for  H3PW  in  acetonitrile  proceed  to  a  very 
small  extent  and  are  given  in  Equations  (2.25)  and  (2.26). 

CH3CN  +  (CH3CN)2H2PW"  (CH3CN)2HV  (CH3CN)HPW^-  (2.25) 
CH3CN  +  (CH3CN)HPW2-  (CH3CN)2H^  +  PW^-  (2.26) 
Reaction  of  the  second  proton  of  H3PW  with  pyridine  involves  reaction  of  the  species  in 
Equation  (2.25).  The  lower  -AH2  than  -AHi  for  all  bases  studied  (e.g.,  1 1.8  vs.  21.0  kcal 
mol"'  for  pyridine)  clearly  shows  that  step  2  involves  a  weaker  acid  than  step  1.  Had  the 
reactions  in  Equations  (2.23)  and  (2.25)  gone  to  completion  when  H3PW  is  added  to 
CH3CN,  -AH2  would  equal  -AHi  because  the  titration  reaction  would  involve 
displacement  of  CH3CN  from  (CH3CN)2H^  (Equation  2.2)  in  both  instances.  This 
resolves  the  literature  differences  by  showing  that  the  first  two  protons  removed  from 
H3PW  are  not  comparable  in  acidity.    Step  2  predominantly  involves  titration  of 
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(CH3CN)2H2PW",  shown  in  Equation  (2.25),  giving  a  smaller  -AH2  than  -AHi  for  all 
bases  studied.  The  question  then  arises,  why  is  a  large  K2  and  a  small  -AH2  observed  for 
all  bases  studied,  (Table  2-4)?  When  base  reacts  with  (CH3CN)2H2PW",  the  endothermic 
loss  of  CH3CN  molecules  could  occur  to  form  HPW^"  as  shown  in  Equation  (2.27). 
(CH3CN)2H2PW  +  B  BHCH3CN^  +  HPW^' +  CH3CN  (2.27) 
With  acetonitrile  being  bound  less  completely  or  not  at  all  to  the  weaker  acid  HPW^',  AS2 
increases  and  -AH2  decreases  for  the  second  step.  Equation  (2.27). 

Equation  (2.27),  a  more  accurate  representation  of  Equation  (2.4),  also  accounts 
for  the  smaller  K  and  larger  -AH  observed  for  step  3.  Since  endothermic  CH3CN 
displacement  is  not  involved  in  step  3,  a  larger  -AH  and  -AS  results  than  for  step  2. 
Equation  (2.28)  best  represents  the  equilibrium  given  in  Equation  (2.5)  for  the  third  step. 
CH3CN  +  B  +  HPW^"  BHCH3CN^  +  PW^"  (2.28) 
Note  that  in  contrast  to  step  1  where  one  pyridine  coordinates,  there  is  an  enthalpy 
contribution  from  one  molecule  of  pyridine  and  one  of  acetonitrile  coordinating  to  the 
displaced  proton  in  step  3.  With  weaker  donors,  step  3  could  involve  formation  of  a 
hydrogen  bonded  adduct  instead  of  PW^"  and  the  equilibria  would  be  incorrectly 
formulated  in  Equation  (2.5)  of  our  analysis.  Both  steps  2  and  3  may  be  more 
complicated  than  shown  above  because  of  reaction  of  BHCH3CN"^  with  added  base  to 
form  B2H''  instead  of  or  in  addition  to  the  reactions  in  Equations  (2.27)  and  (2.28). 
Complexes  of  heteropoly  acids  with  organic  molecules  with  formula  [(B2H"')3][PX^"]  (X 
=  Mo,  W)  have  been  prepared  [65-68].  X-ray  studies  of  these  complexes  show  that  the 
B2H"'  units  do  not  interact  with  the  polyoxoanions.  Spectral  studies  [65-68]  also  suggest 
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that  the  structure  of  these  complexes  in  solution  are  similar  to  those  in  the  solid  state. 
Because  of  the  complications  from  competing  reactions  and  extents  of  CH3CN 
displacement,  our  focus  will  be  on  step  1  for  acidity  comparisons  and  no  further 
interpretation  of  steps  2  and  3  will  be  offered.  Clearly,  our  analysis  has  accomplished  the 
primary  objective  by  demonstrating  that  the  three  protons  of  H3PW  have  very  different 
strengths. 

Reactions  of  H3PW  with  Other  Bases 

Calorimetric  titrations  were  carried  out  using  the  bases  quinoline  (Qi);  N- 
methylimidazole  (M-Im);  quinuclidine  (Qu);  triethylphosphine  (Tep);  tetrahydrofuran 
(THF);  water;  N,N-dimethylacetamide  (DMA);  and  dimethylsulfoxide  (DMSO).  Figure 
2-7  shows  the  calorimetric  isotherms  for  all  bases  titrated.  The  reactions  described  above 
for  pyridine  explain  the  consistent  pattern  of  Ki  >  K2  >  K3  and  -AHi  >  -AH3  >  -AH2  for  all 
bases.  For  THF,  the  total  heat  evolved  is  1.610  calories,  and  the  exotherm  virtually  stops 
when  base  is  added  above  a  2:1  mol  ratio  of  THF  to  H3PW.  Therefore,  K3  and  AH3  can 
not  be  determined  in  CH3CN  with  these  experimental  conditions,  and  the  data  are 
analyzed  with  a  two  step  equilibrium.  When  the  experimental  heats  (hexp)  from  the 
titration  curves  of  DMSO  are  analyzed  with  a  three  step  equilibrium,  meaningless 
equilibrium  constants  and  enthalpies  result.  A  reasonable  fit  is  obtained  with  a  2:1  model 
when  base  additions  up  to  a  2:1  mole  ratio  of  DMSO  to  H3PW  are  used.  This  limited 
data  is  fit  to  a  two  step  equilibrium  giving  values  of  Ki  =  1.6  (±0.6),  K2  =  6.6x10'^ 
(±4.2x10"^),  -AHi  =  22.5  (±0.3),  and  -AH2  =  9.2  (±4.8)  kcal  mol"'  respectively. 
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Figure  2-7.  Isotherms  of  the  Calorimetric  Titrations  of  H3PW12O40  with  Bases.  The 
Bases  are:  QuinoHne  (Qi),  Pyridine  (Py),  Tetrahydrofiiran  (THF),  N-methylimidazole  (M- 
Im),  Triethylphosphine  (Tep),  Quinuclidine  (Qu),  Dimethyl  sulfoxide  (DMSO),  and  N,N- 
dimethylacetamide  (DMA). 
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Figure  2-8  shows  an  excellent  linear  relation  of  log  Ki  to  -AHi  for  all  bases  except 
DMA,  water  and  DMSO.  DMA  and  water  lie  above  the  line  while  DMSO  falls  below. 
DMSO  and  DMA  also  show  a  poor  correlation  of -AHi  to  the  ECW  model.  Discussion  of 
the  complications  encountered  with  these  donors  and  limitations  on  the  interpretation  of 
these  enthalpies  will  be  treated  in  the  section  on  the  ECW  analysis.  The  large 
experimental  errors  for  K2  and  the  small  range  of  K3  values  preclude  meaningful  log  K 
vs.  -AH  plots  for  steps  two  and  three. 

Enthalpy  Analysis  with  the  ECW  Model 

The  ECW  approach  has  been  successfully  employed  [54-56]  to  correlate  many 
neutral  base  neutral  acceptor  adduct  formation  enthalpies  to  the  equation: 
-AH  =  EaEb  +  CaCb  +  W  (2.29) 
EaEb  and  CaCb  describe  electrostatic  and  covalent  contributions  for  the  donor  (B)- 
acceptor  (A)  reaction.  The  W  term  is  zero  for  simple  donor-acceptor  adduct  reactions  in 
poorly  solvating  media,  but  for  a  more  complex  reaction,  it  incorporates  any  constant 
contributions  of  a  particular  acid  that  is  independent  of  the  base  it  reacts  with.  Equation 
(2.29)  has,  for  the  most  part,  been  applied  to  1 : 1  adduct  formation  enthalpies.  Recently,  a 
successful  correlation  of  the  enthalpies  of  reaction  of  CF3SO3H  with  various  bases  [70]  to 
form  an  ion-pair  in  dicloroethane  has  been  reported  [71].  This  encouraged  us  to  attempt  a 
correlation  of  the  -AHi  values  for  the  ionic  products  of  the  reaction  of  a  Bronsted  acid 
with  the  bases  in  Table  2-4.  The  result  is  an  excellent  data  fit  to  Equation  (2.29)  shown 
in  Table  2-5  for  most  of  the  bases  studied  giving  Ea*  =  2.56  ±0.23,  Ca*  =  7.38  ±0.07,  and 
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Figure  2-8.  Log  Ki  versus  -AHi  (kcal  mol'^)  for  the  Reaction  of  H3PW  with  Various 
Bases  in  CH3CN  Solvent.  Bases  are:  Quinuclidine  (Qu),  Triethylphosphine  (Tep),  N- 
methyhmidazoie  (M-im),  Pyridine  (Py),  Quinoline  (Qi),  and  Tetrahydrofuran  (THF) 
respectively. 
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Table  2-5.  Results  of  ECW  Fits  for  Reactions  of  H3PW  -  Bases  in  CH3CN  Solvent." 


Base 

Cb/Eb 

-AHexp 

-AHecw'' 

Pyridine 

1.99 

21.0 

21.3 

Quinoline 

1.27 

18.1 

17.8 

N-methylimidazole 

4.24 

29.5 

29.9 

THF 

1.33 

11.4 

10.9 

Quinuclidine 

8.40 

42.6 

42.2 

Triethylphosphine 

19.75 

32.2 

32.1 

Acetonitrile 

0.43 

0.0 

0.07 

a.  Unit  for  AH  is  kcal  mol"'. 

b.  Calculated  with  2.56  (±0.22)  Eb  +  7.38  (±0.07)  Cb  -  9.36  (±0.0).  Values  of  Eb  and  Cb 
for  the  bases  listed  above  can  be  obtained  in  reference  [56].  Parameters  for 
triethylphosphine  are  Eb  =  0.28,  Cb  =  5.53. 
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W  =  -9.36.  The  Ea*  and  Ca*  parameters  characterize  the  acceptor  properties  of 
CHsCNH"^  in  acetonitrile.  The  W  value  incorporates  constant  net  solvation  energies  of 
products  and  reactants  and  the  endothermic  enthalpy  of  dissociation  of  the  anion  from 
(CH3CN)3H3PW. 

The  Ea*  and  Ca*  values  of  2.56  and  7.38  for  CHsCNH"^  compared  to  the  reported 
value  of  Ea*  and  Ca'  of  4.51  and  5.70  for  CF3SO3H  in  dichloroethane  suggest  that  the 
former  is  more  covalent  and  much  less  electrostatic  than  CF3SO3H.  With  different 
Ca*/Ea*  ratios  for  the  different  acids  involved,  the  enthalpies  of  reaction  for  a  series  of 
different  bases  reacting  with  H3PW  in  CH3CN  and  CF3SO3H  in  1 ,2-dichloroethane  are 
not  expected  to  plot  up  linearly. 

The  fit  of  the  H3PW  enthalpies  to  ECW  is  an  important  result  for  it  indicates  that 
the  same  base  properties  that  correlate  neutral  donor-acceptor  reactions  in  poorly 
solvating  solvents  also  determine  the  basicity  order  when  ionic  products  are  formed  from 
Bronsted  acids  in  acetonitrile.  The  acid  concentration  should  be  held  constant  to  establish 
basicity  orders  because  variation  in  the  extent  of  ion-pairing  with  acid  concentration  will 
influence  the  enthalpy  measured.  The  ECW  fit  shows  that  enthalpy  contribution  from  the 
basicity  and  polarity  of  acetonitrile  solvent  in  facilitating  ionization  of  HX  is  either  a 
constant  factor  or  is  linear  with  basicity.  As  a  result,  with  the  acid  concentration  constant, 
-AH  for  reaction  of  a  base  in  this  solvent  provides  a  measure  of  the  donor  strength  and 
substantiates  earlier  discussion  on  the  selection  of  acetonitrile  as  a  solvent  for  studying 
Bronsted  acidity. 
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The  enthalpy  data  for  DMSO,  H2O  and  DMA  do  not  obey  the  ECW  correlation 
and  reveal  some  interesting  conclusions  regarding  reactivity.  The  DMSO  reaction  could 
involve  bonding  to  sulfur  which  is  expected  for  acceptors  with  a  large  Ca'/Ea*.  With  the 
Eb  and  Cb  values  reported  [56]  for  DMSO  referring  to  oxygen  coordination,  even  an 
equilibrium  mixture  of  S  and  O  coordination  would  give  an  average  -AH  larger  than  that 
predicted  by  ECW.  Unfortunately  overlap  of  the  S=0  stretching  frequency  with  H3PW 
and  CH3CN  bands  precludes  testing  this  explanation.  Water  was  not  included  in  the 
ECW  analysis  because  Eb  and  Cb  parameters  are  only  available  for  the  gas  phase 
monomer  and  these  cannot  be  used  for  water  aggregates  in  solution.  Solutions  of 
heteropoly  acids  in  DMA  are  reported  to  be  photochemically  redox  active  [65-68].  The 
solution  resulting  from  our  titration  of  H3PW  with  DMA  has  a  bluish  tint  and  we  attribute 
the  discrepancy  with  ECW  to  enthalpy  contributions  from  redox  reactions. 

The  ECW  analysis  is  not  only  a  correlation  and  predictive  tool  but  can  provide 
thermodynamic  data  for  reactivity  that  carmot  be  measured  directly.  Using  Ea*  and  Ca* 
with  Eb  and  Cb  for  CH3CN,  the  enthalpy  of  coordinating  CH3CN  to  (CH3CN)Fr  in 
acetonitrile  is  calculated  to  be  9.4  kcal  mol'V  Comparing  this  quantity  to  the  W  value  of 
9.4  indicates  that  (CH3CN)H''  coordinates  acetonitrile  and  (CH3CN)2H2PW  with  equal 
strength,  so  displacement  of  the  anion  by  acetonitrile  would  have  a  zero  enthalpy 
contribution.  Adding  W  of  9.4  kcal  mol"'  to  -AHi  measured  in  acetonitrile,  gives  -AH  for 
reactions  of  bases  with  CH3CNH^  as  an  acid  in  acetonitrile. 
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An  Acidity  Scale  in  CHiCN  toward  Pyridine 

In  order  to  provide  a  comparison  of  relative  acidity,  the  enthalpies  of  reactions  of 
a  series  of  other  acids  with  pyridine  were  measured  in  CH3CN.  Results  of  the  calculated 
K  s  and  AH  s  are  presented  in  Table  2-6.  Since  the  position  of  the  ionic  equilibria  in 
Equation  (2.1)  will  be  concentration  dependent,  an  acid  concentration  of  1.39x10"^  M  is 
used  in  all  experiments  and  is  a  defined  quantity  for  our  acetonitrile  acidity  scale. 
Because  of  variation  in  ion-pairing,  differences  less  than  ±0.5  kcal  mol"'  are  not 
considered  significant.  The  -AHi  values  fiimish  an  acidity  order  : 
H3PW  >  CF3SO3H  >  P-CH3C6H4SO3H  ~  H2SO4  (97%) »  CF3COOH  »  CIC6H4COOH. 
A  plot  of  log  K]  vs.  -AH]  for  the  acids  studied  in  Table  2-6  is  shown  in  Figure  2-9.  The 
linear  firee  energy-enthalpy  relation  supports  a  similar  reaction  type  for  all  acids.  H3PW  is 
the  strongest  acid  that  can  be  measured  with  this  scale  because  an  acid  that  is  completely 
ionized  in  acetonitrile  will  exist  as  (CH3CN)2H^  and  it  was  shown  earlier  that  acetonitrile 
and  (CH3CN)2H2PW"  have  similar  enthalpies  of  displacement  from  (CH3CN)H^. 

An  enthalpy  for  the  reaction  of  pyridine  with  P-CH3C6H4SO3H  in  acetonitrile  is 
reported  [72].  At  low  base  concentration,  a  -AH  of  18.3  kcal  mof'  results  fi-om  a  limiting 
reagent  estimate  in  good  agreement  with  our  value  of  18.7.  When  the  limiting  reagent 
approach  is  used  with  a  ten  fold  excess  of  acid  an  enthalpy  of  18.6  kcal  mol"'  is  reported 
[72].  The  -AH  is  reported  to  decreases  with  base  concentration  and  this  is  attributed  to  a 
base  dependent  heat  of  reaction.  Our  magnitude  for  K  shows  this  behavior  is  expected 
from  incomplete  protonation  of  all  the  base  added  as  a  1:1  ratio  is  approached  and  the 
true  enthalpy  of  product  formation  is  base  concentration  independent. 
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Table  2-6.  Enthalpies  (kcal  mol'')  and  Equilibrium  Constants  (25±1°C)  for  the  Reaction 
of  Pyridine  with  Various  Acids"  in  CH3CN  Solvent.  [Acidji  =  1.39x10"^  M. 


Acid 

Ki 

K.2  K.3 

-AHi 

-AH2 

-AH3 

H3PW 

2.1x10^ 
(±4.0x10') 

2.4x10^  4.1x10' 
(±3.4x10')  (±5.8) 

21.0 
(±0.2) 

11.8 

(±0.3) 

18.6 

(±0.5) 

H2SO4'' 

1.5x10^ 
(±2.0x10') 

7.9 

(±5.8) 

18.4 

(±0.2) 

0.2 
(±0.5) 

CF3SO3H 

6.1x10^ 
(±5.0x10') 

19.9 
(±0.2) 

CH3C6H4SO3H 

2.1x10^ 
(±5.5x10') 

18.7 

(±0.2) 

CF3COOH 

2.8x10' 
(±1.2x10') 

14.2 

(±1.0) 

CIC6H4COOH 

2.1x10"' 
(±1.0x10') 

6.5 
(±1.9) 

a.  Standard  deviations  are  based  on  variance-covariance  matrix  analysis. 

b.  Concentrated  sulfuric  acid  used  for  preparation  of  standard  solution  of  H2SO4  in 
acetonitrile  was  97%  w/w,  calculated  by  standardization. 
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Figure  2-9.  Log  Ki  versus  AHi  (kcal  mol"')  for  the  Reaction  of  Pyridine  with  Various 
Acids  in  CH3CN  Solvent.  Acids  are:  1 2-tungstophosphoric  (H3PW),  triflic  (TRIP),  p- 
toluenesulfonic  (PTS),  sulfuric  (H2SO4),  trifluoroacetic  (TP A),  and  m-chlorobenzoic 
(CB)  respectively. 
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The  reaction  of  pyridine  with  triflic  acid  in  CH3CN  gives  an  enthalpy  of  19.9  kcal 
mol"'.  In  acetonitrile  the  enthalpy  of  reaction  of  (CHaCN)!!""  with  pyridine  is  calculated 
as  -30.4  kcal  mol''.  The  difference  in  the  measured  enthalpy  and  this  quantity  gives  10.5 
kcal  mol''  as  the  enthalpy  of  dissociating  the  CF3SO3'  from  CHsCNH""  in  acetonitrile 
compared  to  the  W  value  of  9.4  for  dissociating  (CH3CN)2H2PW".  Assuming  that  all 
acids  in  the  above  series  form  (CH3CN)H(C5H5N)'^,  similar  calculations  for  the  other 
acids  provide  an  enthalpy  scale  of  conjugate  base  strengths  toward  this  acid. 

An  extensive  compilation  of  the  dissociation  constants  (pKg)  in  acetonitrile  is 
available  [69].  For  the  acids  employed  in  this  work,  pKa  of  7.8,  8.7,  and  12.6  are  given 
for  sulfuric,  p-toluenesulfonic,  and  trifluoroacetic,  respectively.  The  pKa  trend  is  in  good 
agreement  with  our  enthalpy  trend  which  provides  a  direct  measure  of  bond  strength. 

Comparison  of  Solution  and  Solid  Calorimetric  Data 

Comparison  of  the  solution  enthalpies  with  those  from  gas-solid  calorimetry  is 
informative.  As  reported  [43],  the  gas-solid  data  contains  a  dispersion  component  not 
present  in  solution  leading  to  larger  measured  -AH  values  than  those  corresponding  to  the 
donor-acceptor  component.  Furthermore,  the  differential  heat  of  gaseous  ammonia 
adsorption  at  50  °C  for  the  solid  sample  treated'^  at  250  °C  is  claimed  to  show  that  all 
three  protons  of  the  solid  are  titrated  producing  the  same  -AH  of  s  35  kcal  mol"'.  There  is 
no  solvent  leveling  or  displacement  in  the  gaseous  ammonia  reaction  but  this  can  not 
account  for  three  protons  of  equivalent  acidity.  Unfortunately  equilibrium  constants  are 
not  measured  in  reported  gas-solid  calorimetric  studies  [43].  If  Ki,  K2,  and  K3  decrease 
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slightly  for  the  gas  phase  reaction,  simultaneous  reaction  of  all  three  protons  would  occur 
leading  to  averaged  heat  evolution  and  an  average  -AH  of  35  kcal  mol"'.  The  -AH  could 
in  reality  correspond  to  the  average  of  large,  medium  and  small  enthalpies  with  sunilar 
K's.  Slightly  decreasing  equilibrium  constants  would  account  for  a  -AH  slightly  higher 
than  35  with  0.1  moles  of  NH3  to  1  mole  of  H3PW  added  and  the  small  -AH  of  «  20  kcal 
mol"'  based  on  a  limiting  reagent  assumption  for  the  remaining  weakest  site  up  to  a  3:1 
NH3  to  H3PW  mole  ratio.  Titration  of  samples  with  varying  water  content  are  reported'' 
[40]  but  as  described  above  (Table  2-2)  interpretation  of  these  measurements  is  difficult. 
Calorimetric  titration  with  ammonia  [42]  at  150  °C  on  samples  activated  at  150  and  200 
°C  are  also  reported.  Reaction  is  practically  complete  for  2.89  moles  of  added  base  with 
the  differential  heat  decreasing  fi-om  48  to  36  kcal  mol"'  during  the  titration.  At  this 
higher  temperature  the  step  equilibrium  constants,  for  steps  with  decreasing  enthalpies, 
will  be  closer  in  magnitude  than  at  lower  temperatures.  Without  a  knowledge  of  the  three 
step  equilibrium  constants,  the  interpretation  of  the  results  from  a  gas-solid  enthalpy 
titration  is  ambiguous  because  the  measured  heat  evolved  cannot  be  resolved  into 
contributions  from  the  different  steps  [43]. 

Solution  calorimetry  unequivocally  shows  the  three  protons  are  not  strong  acids  of 
equal  strength  and  the  discussion  offered  in  this  work  is  relevant  to  measuring  acidity  in 
homogeneous  systems.  Solid  acids  have  the  conjugate  base  anion  as  part  of  the  crystal 
lattice  and  proton  transfer  from  a  Bronsted  site  produces  an  intimate  ion-pair.  The  lattice 
energy  of  the  ion-pair  compared  to  that  of  the  solid  acid  contributes  to  -AH.  Thus,  any 
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direct  comparison  of  enthalpies  for  homogeneous  and  heterogeneous  systems  should  take 
into  account  the  differences  in  the  reactions  involved. 

Conclusions 

This  study  provides  thermodynamic  parameters  for  the  three  step  equilibria  of 
H3PW  reacting  with  a  series  of  bases  in  acetonitrile  solution.  The  results  indicate  that  the 
three  protons  are  not  of  equal  strength  as  proposed  from  gas-solid  calorimetry 
measurements.  Only  one  strong  proton  per  molecule  of  acid  is  furnished  for 
stoichiometric  or  catalytic  reactions  in  any  solvent  of  basicity  equal  to  or  less  than 
CH3CN.  Consistent  with  literature  claims  of  superacidity,  the  enthalpies  in  acetonitrile 
show  that,  as  a  monoprotic  acid  toward  pyridine,  H3PW  is  comparable  to  triflic  acid  and 
is  stronger  than  sulfuric  acid,  and  p-toluenesulfonic  acid. 

The  enthalpies  of  reactions  for  the  first  protonation  process  were  successfiilly 
correlated  to  the  ECW  model,  leading  to  substantial  covalent  character  (Ca*  /Ea*  =  2.88) 
in  the  donor-acceptor  reactions  of  the  species  (CH3CN)H^.  With  a  database  of  more  than 
one  hundred  bases  in  the  ECW  correlation,  enthalpies  of  reaction  of  all  these  bases  with 
H3PW  can  be  predicted.  Of  even  more  significance  is  the  ability  of  the  solvation 
minimized  ECW  parameters  for  neutral  donor  acceptor  interactions  to  predict  enthalpies 
for  the  reaction  of  Bronsted  acids  in  acetonitrile  solvent.  This  justifies  arguments 
presented  for  using  enthalpies  in  acetonitrile  as  a  Bronsted  acidity  scale. 

The  dehydration  of  H3PW  is  complicated  by  decomposition,  and  it  must  be  dried 
carefully  in  order  to  obtain  the  anhydrous  acid.  A  procedure  for  accomplishing  this  is 
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reported.  The  demonstrated  instability  of  H3PW  at  300  °C  precludes  its  use  as  an  acid 
catalyst  at  this  or  higher  temperatures. 

In  any  acid  catalyzed  reaction  occurring  in  the  condensed  phase,  the  reactivity  will 
depend  upon  the  number  of  protons  provided  as  well  as  the  strength  of  interaction  of  the 
acid  with  the  substrate.  Hammett  indicators,  which  provide  a  measure  of  the  number  of 
acid  species  supplied,  do  not  measure  the  strength  of  the  acid-substrate  bond. 
Calorimetry  at  low  base  concentration,  assuming  complete  complexation,  can  provide  a 
measure  of  the  strength  of  the  interaction  for  monoprotic  acids  or  when  Ki>K2  for 
diprotic  acids.  The  determination  of  both  K  and  -AH  furnishes  a  complete 
characterization  of  the  strength  and  amount  of  acid  species.  With  a  linear  log  K  vs.  -AH 
plot  (Figure  2-9),  the  trend  in  amount  and  strength  are  the  same  for  H3PW  in  acetonitrile. 
Thus,  the  calorimetric  measurement  of  -AH  for  different  acids  in  acetonitrile  solvent 
provides  a  temperature  independent  Bronsted  acidity  scale  of  relevance  to  homogeneous 
catalytic  reactions  in  which  unusual  entropy  effects  are  absent.  We  emphasize  that  this  is 
a  one  parameter  scale  and  should  be  applied  in  a  maimer  consistent  with  the  limitations 
[73]  of  such  a  scale. 


CHAPTER  3 

ACIDITY  OF  H3PW12O40  AND  CSXH3.XPW12O40  IN  HETEROGENEOUS 

SYSTEM 

Introduction 

Research  in  heteropoly  acids  has  focused  attention  in  their  industrial  use  as 
homogeneous  and  heterogeneous  catalysts  for  redox  and  acid  catalyzed  reactions  [32a]. 
H3PW12O40  (abbreviated  as  H3PW)  is  the  strongest  heteropoly  acid,  and  it  has  been 
characterized  recently  as  a  monoprotic,  superacid  in  solution  by  calorimetry  [74] 
resolving  conflicting  reports  from  NMR  [52],  IR  [48-51],  UV-vis  [64],  thermal  methods 
(TG,  DSC,  TPR,  TPD)  [44-45,  47],  microcalorimetry  with  ammonia  absorption  [40-42], 
Hammett  acidity  constants  [37],  conductivity  [36],  and  the  '^C  shift  of  mesityl  oxide  [38]. 

The  crystal  structure  of  H3PW12O40  [4a,  4b]  is  demonstrated  with  X-ray  and 
neutron  diffraction  [75]  to  be  a  hexahydrate.  The  structure  can  be  divided  into  primary 
(anion)  and  secondary  (protons  and  water)  parts.  The  complete  structure  of  the 
phosphotungstic  acid  has  been  described  in  more  details  in  chapter  2. 

Recently  CSXH3.XPW12O40  (abbreviated  as  CsxHs.xPW)  derivatives  have  gained 
attention  because  of  their  claimed  superacidity  and  shape-selectivity  in  alkylations 
reactions  [76-78].  Unlike  H3PW,  the  CsxHs.xPW  salts  have  limited  solubility  in  water  or 
organic  solvents  and  are  used  as  heterogeneous  catalysts.  They  show  higher  catalytic 
activity  and  higher  thermal  stability  in  some  reactions  than  the  parent  acid  [27].  In 
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CS3PW,  the  Cs"^  occupies  the  positions  of  H502^  in  the  hydrate  [27]  having  the  same 
cubic  structure  of  hexahydrated  H3PW,  but  with  a  smaller  lattice  constant.  The  structure 
of  the  other  CsxHa.xPW  are  not  completely  established.  The  proposed  structure  of  cesium 
complexes  with  x<  3  are  reported  to  be  mixtures  of  H3PW  and  CsxHa.xPW.  A  study  of 
the  precipitates  [79]  from  the  reaction  of  H3PW  with  CS2CO3  evaporated  to  dryness 
together  with  the  solution  concludes  that  in  the  range  of  0  <  x  <  2  are  mainly  CsiHPW. 
Compounds  prepared  in  the  range  of  2  <  x  <  3  are  CS3PW,  although  Cs2  5H0  5PW  is 
considered  a  mixture  of  CS2HPW  and  CS3PW.  Reported  XRD  spectra  for  CsxH3.xPW  (x 
=  0,  1,  2,  and  3)  [80]  claims  that  CSH2PW  is  a  mixture  of  H3PW  and  CS2HPW  which  is 
transformed  in  a  nearly  homogeneous  acidic  salt  after  thermal  treatment  of  the  precipitate 
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at  300  C.  The  lattice  constants  values  for  those  cesium  salts  indicate  all  of  them  belong 
to  cubic  system.  The  XRD  pattern  is  very  similar  for  CS2HPW  and  CS3PW,  but  distinct 
enough  from  H3PW  and  CSH2PW.  The  CsxH3.xPW  (x  =  2,  2.5,  and  3)  were  studied  by 
^'P  NMR  [81].  The  ^'P  spectrum  of  CS2HPW  has  four  peaks,  which  were  interpreted  as  a 
mixture  of  four  species  (x  =  0,  1,2,  and  3).  The  results  for  CS2.5H0  5PW  claimed  that  this 
compound  is  a  mixture  of  CS2HPW  and  CS3PW.  Supported  H3PW  on  CS3PW  presents 
the  same  ^'P  spectrum  as  CS2.5H0  5PW  [81],  accounting  to  the  exchange  between  H"^  and 
Cs^  during  the  heating  process.  In  all  those  studies  [79-81]  the  cesium  derivatives  were 
prepared  similarly  [77],  though  the  thermal  treatment  before  analysis  of  the  solids  were 
not  always  the  same. 

The  tertiary  structures  of  these  cesium  salts  are  important.  The  tertiary  structure  is 
composed  of  the  size  of  the  primary  and  secondary  particles,  distribution  of  protons  and 
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cations,  and  pore  structure  [27].  Porosity  is  important  in  catalysis  [80]  and  leads  to 
shape-selectivity  [82].  There  is  a  remarkable  change  of  surface  area  in  going  from  H3PW 
to  CS3PW.  H3PW  has  a  surface  area  of  5-6  m'^  g  ',  which  decreases  slightly  as  the  cesium 
content  is  increased  up  to  CSH2PW  (1  m^  g'').  Above  two  cesiums,  the  surface  area 
increases  rapidly,  for  instance,  CS2.5H05PW  is  reported  to  have  135  m^  g'  of  area  [77]. 
Micropores  form  in  these  solids  [77,  82]  when  the  proton  is  substituted  by  cesium. 
Adsorption  studies  with  different  probes  (e.g.,  nitrogen,  benzene,  etc.)  [77]  give  the  pore 
sizes  for  Cs2.2H0.8PW  as  6.2-7.5  A  and  those  for  Cs2.5H0.5PW  are  larger  than  8.5  A. 
Microporosity  studies  of  CSxHs.xPW  carried  out  by  '^^Xe  NMR  [82]  indicate  slightly 
smaller  pores. 

Studies  of  the  microporosity  of  monovalent  sahs  of  the  heteropoly  compounds 
have  been  reported  [83-86]  using  nitrogen  adsorption-desorption  and  X-ray  diffraction 
analysis.  The  solids  are  reported  to  contain  interstitial  voids  set  apart  from  each  other  by 
the  terminal  oxygen  atoms  of  the  polyanions.  Widening  of  the  interstices  occurs  as  the 
size  of  the  cation  increases.  These  studies  demonstrate  that  the  pore  sizes  can  be 
controlled  in  heteropoly  acids  substituting  the  proton  with  different  cations. 

In  an  earlier  report  [74],  thermodynamic  data  describing  the  acidity  of  H3PW  was 
measured  in  acetonitrile  solution  by  calorimetric  titration.  The  Keggin  structure  of 
crystalline  H3PW  is  maintained  in  solution  [2].  This  chapter  aims  to  characterize  the 
acidity  of  solid  H3PW  and  the  solid  CsxH3.xPW  salts  in  their  reactions  with  pyridine  in 
cyclohexane  using  combined  measurements  from  a  calorimetric  titration  and  adsorption 
isotherm  measurement,  Cal-ad  [87,  59,  43].  A  multiple-site  Langmuir  equation  is  used  to 
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solve  both  measurements  simultaneously  and  calculates  the  number  of  sites  (nO, 
equilibrium  constants  (Ki),  and  enthalpies  (AHj)  for  different  acid  sites  on  these  solids. 
Thus,  the  Cal-ad  analysis  provides  a  total  characterization  of  solid  acid  strength  and 
furnishes  enthalpies  for  different  sites  to  provide  a  temperature  independent  acidity  scale. 
Differences  in  the  reactivity  of  solid  H3PW  and  acetonitrile  solutions  of  this  acid  are 
discussed. 

Experimental 
Materials  and  Preparation  of  Cesium  Derivatives 

Elemental  and  TG  analysis  of  the  H3PW  (purchased  from  Aldrich)  revealed  16 
moles  of  water  per  mole  of  H3PW.  Cyclohexane  obtained  from  Fisher  was  purified  by 
drying  over  4A  molecular  sieves  for  24  hours,  and  then  distilled  over  P2O5.  Pyridine 
obtained  from  Fisher  was  distilled  over  CaH2.  Dried  samples  were  stored  in  a  container 
with  4A  molecular  sieves.  CsxHs-xPW  (x=l,  2,  2.5,  and  3)  was  prepared  by  adding 
predetermined  amounts  of  a  0.2  mol  L"'  aqueous  solution  of  CS2CO3  to  a  0.08  mol  L'' 
solution  of  H3PW  at  a  rate  of  1  mL  min'  at  room  temperature.  A  white  precipitate 
formed  in  the  solutions.  The  solutions  were  left  to  stand  overnight  at  room  temperature. 
Then,  according  to  the  literature  [77],  the  solutions  were  evaporated  to  dryness  at  40  °C. 
It  is  interesting  to  note  that  as  the  cesium  content  increases  in  the  compounds,  the  product 
appears  whiter.  The  extent  of  proton  exchange  by  cesium  is  in  close  agreement  with  the 
expected  stoichiometric  value,  as  determined  by  ICP-MS  (Table  3-1). 
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Table  3-1.  Chemical  Analysis"  and  Surface  Area''  of  CsxHa.xPW  Compounds. 

Sample        Cs  per  K.U.  (x)     Surface  Area  (m^  g"') 
CSH2PW        03  i 
CS2HPW         1.8  63 
Cs2.5H0.5PW     2.3  119 
CS3PW  2.7  92 

a.  Cesium  contend  was  determined  by  ICP-MS  after  dissolving  the  compounds  in  0.01 
mol  L''  NaOH  solutions. 

b.  Siu-face  area  values  were  obtained  in  reference  78. 
Drying  Procedures 

H3PW  was  dried  at  162  °C  for  4  hours  under  vacuum  [74].  The  CsxHs.xPW  dried 
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at  200  C  for  4  hours  under  vacuum  gave  higher  heats  of  reaction  with  pyridine  in 
cyclohexane  than  at  300  °C. 

The  apparatus  for  drying  H3PW  and  CsxHs.xPW  consisted  of  a  glass  fixed  bed 
reactor  with  two  stopcocks  at  the  ends  to  close  the  system.  An  external  small  opening 
under  the  reactor  bed  contains  a  wire  thermocouple  to  control  and  measure  the 
temperature  of  the  reactor.  The  sample  is  weighed,  placed  into  the  reactor,  and 
evacuated.  The  reactor  is  filled  with  argon  and  evacuated  two  to  three  times  during  the 
drying  period  to  remove  equilibrium  water  vapor  from  the  reactor.  After  4  hours,  the 
sample  is  cooled  to  room  temperature,  the  reactor  is  filled  with  argon  and  transported  to  a 
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inert  atmosphere  box.  The  samples  are  stored  in  vials  in  a  dry  box  until  used.  The  same 
apparatus  has  been  successfully  used  to  dry  zeolites  in  this  laboratory  [43]. 

Spectral  Analysis 

A  Perkin-Elmer  Lambda  6  UV-vis  spectrophotometer  and  quartz  cells  of  1 .0  cm 
pathlength  were  employed  for  the  adsorption  experiment.  Infrared  spectra  were  obtained 
with  a  Nicolet  Model  5PC  FTIR  spectrophotometer  in  nujol  (Aldrich)  mulls  on  NaCl 
plates.  Mulls  of  dried  samples  were  prepared  inside  the  inert  atmosphere  box.  The 
determination  of  cesium  was  measured  for  '^^Cs  by  inductively  coupled  plasma  mass 
spectrometry  with  a  Finningan  MAT  SOLA  ICP-MS.  The  powder  X-ray  diffraction  data 
were  obtained  with  a  Siemens  General  area  detector  (XRDS)  with  x-ray  tube  power  at  50 
kV  and  40  mA  (adding  time  of  10000  seconds). 

The  ^'P  NMR  spectra  of  H3PW,  HsPW-Py,  and  CsxHs-xPW  compounds  were 
taken  on  a  Unity  500  NMR  using  a  Varian  MAS  probe  tuned  to  202.640  MHz.  Samples 
were  loaded  into  silicon  nitride  rotors.  Packing  of  the  rotor  with  the  solids  was 
accomplished  by  filling  the  rotor  2/3  of  the  way  full  with  the  sample,  and  then  packing 
the  remainder  with  silica  gel.  This  prevented  the  samples,  specially  Cs2.5H0.5PW,  from 
sticking  to  the  inside  of  the  rotor,  and  allowed  the  rotor  cap  to  be  removed  after 
acquisition.  The  magic  angle  was  tuned  using  K^^Br.  A  recycle  delay  of  2  seconds  was 
used.  The  acquisition  time  was  set  to  0.02  s,  with  the  FID  no  longer  showing  the 
evolution  of  signal  after  0.01  s.  Signals  were  indirectly  referenced  to  85  wt.%  phosphoric 
acid.  One  hundred  plus  transients  were  obtained. 
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Calorimetric  Titrations  and  Adsorption  Measurements 

For  solid  acids,  a  pyridine  solution  is  added  to  a  slurry  of  the  solid  in  cyclohexane, 
and  the  measurement  of  heat  evolved  and  the  equilibrium  amount  of  base  in  solution  is 
determined  in  two  independent  experiments.  The  calorimetric  data  generates  an  isotherm 
of  total  heat  evolved  versus  the  total  moles  of  base  added.  Calorimetry  does  not  account 
for  the  amount  of  base  in  solution,  but  the  adsorption  measurement  provides  this 
information  producing  an  isotherm  of  moles  of  base  adsorbed  versus  moles  of  base  in 
solution.  It  is  important  that  the  concentrations  of  base  added  and  the  ratio  (V/g)  between 
the  volume  of  solution  (V)  and  the  mass  of  solid  used  (g)  be  the  same  in  both  calorimetric 
and  adsorption  experiments  to  interchange  the  data.  Samples  of  H3PW  or  CSxHa-xPW,  are 
weighed  and  transferred  to  a  vacuum  insulated  calorimetric  cell  containing  a  stir  bar.  For 
each  titration,  100  mL  of  cyclohexane  is  added  to  the  cell.  A  calibrated  syringe,  filled 
with  a  solution  of  known  concentration  of  pyridine  (e.g.,  0.1  mol  L"')  is  inserted  into  the 
cell  along  with  a  thermistor  and  a  heater  coil.  All  these  operations  are  carried  out  in  an 
inert  atmosphere  glove  bag.  The  thermistor  and  heater  coil  are  cormected  to  an  electronic 
bridge,  and  a  computer  [59].  The  thermistor  is  calibrated  with  the  heater  coil  prior  to  or 
immediately  after  each  titration.  A  set  of  12  calibrated  brass  stops  are  used  for  the 
syringe  additions.  After  each  addition  the  heat  evolved  fi-om  reaction  of  the  base  with  the 
samples  is  measured.  The  heat  from  adding  the  base  solution  to  the  solvent  is  negligible. 
Each  titration  was  repeated  two  times. 

The  adsorption  experiments  are  carried  out  using  1  g  of  solid,  weighed  inside  the 
dry  box,  and  added  to  a  sealed  three  neck  round  bottom  flask  with  100  mL  of 
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cyclohexane.  The  same  syringe  and  calibrated  brass  stops  used  in  the  calorimetric 
experiment  are  used  for  addition  of  pyridine  in  this  experiment.  After  each  addition  of 
base,  at  least  3  min  is  allowed  for  equilibrium  to  be  reached.  Then,  a  1  mL  sample  of 
solution  is  removed  from  the  flask  and  placed  into  a  quartz  cell  of  1  cm  pathlength  and  1 
mL  volume.  One  milliliter  of  cyclohexane  is  added  back  into  the  flask  in  order  to 
maintain  a  constant  volume.  The  absorbance  of  pyridine  is  measured  at  251  nm  to 
determine  its  equilibrium  concentration  in  solution.  Since  the  amoimt  of  pyridine  added 
is  known,  the  amount  of  base  adsorbed  by  the  solid  is  calculated  by  difference.  Each 
adsorption  experiment  was  repeated  two  times. 

Cal-ad  Calculations 

The  Cal-ad  analysis  is  based  on  simultaneously  solving  data  from  calorimetric  and 
adsorption  experiments  [87,  59,  56,  43].  The  model  employs  multiple-site  equilibrium 
for  the  different  acid  sites  present  in  a  solid  acid.  A  Langmuir-type  equation  is  summed 
over  all  sites  of  the  solid,  according  to  Equation  (3.1): 

h/g  =  I  (niK,[B]/(l  +  Ki[B]))AHi  (3.1) 
where  h  is  the  sum  of  the  heat  evolved  (cal)  from  the  calorimetric  tifration;  g  is  the  mass 
of  solid  (g);  [B]  is  the  concentration  of  base  in  solution  at  equilibrium  (mol  L'*);  nj  is  the 
number  of  each  different  site;  Kj  is  the  equilibrium  constant  for  each  site;  and  AHj  is  the 
enthalpy  of  reaction  for  each  site  on  the  solid. 

A  Simplex  routine  is  used  to  calculate  the  parameters.  The  data  adjusted  to  the 
model  consist  of  the  mol  adsorbed  and  equilibrium  concentration  in  solution  (adsorption), 
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and  the  sum  of  the  heat  and  equilibrium  concentration  in  solution  (calorimetry).  As  the 
calorimetric  titration  does  not  measure  the  equilibrium  concentration  in  solution,  this 
information  is  obtained  through  a  polynomial  series,  according  to  this  expression: 
(V/g)[T]  =  x"[B]"  +  x"-'[B]"-'  +  x"-2[B]"-'  +  x"-3[B]"-3  +  ...  (3.2) 
where  (V/g)  is  the  ratio  of  volume  of  solution  (L)/mass  of  solid  (g);  [T]  is  the  total 
concentration  of  base  added  for  each  addition  in  the  calorimetric  experiment  (mol  L"'); 
[B]  is  the  equilibrium  concentration  of  base  in  solution  (mol  L"').  The  polynomial  is  one 
degree  higher  than  the  number  of  sites  (e.g.,  for  two  sites  a  third  degree  polynomial  is 
used  as  shown  in  Equation  3.2).  The  x"  terms  have  the  following  meaning  for  a  two  site 
adjustment:  x^  =  KiK2(V/g);  x^  =  n,  K1K2  +  nj  K,K2  -  (V/g)KiK2[T]  +  (V/g)K,  + 
(V/g)K2;  and  X  =  n,K,  +  n2K2  -  (V/g)Ki[T]  -  (V/g)K2[T]  +  (V/g).  Estimates  of  ns  and 
K  s  are  obtained  from  a  Langmuir  analysis  of  the  adsorption  data  using  Equation  (3.3): 
STB  =  Zn,Ki[B]/(l+Ki[B])  (3.3) 
where  StB  is  the  total  number  of  moles  of  base  adsorbed  per  gram  of  solid.  Cal-ad  has 
been  successfully  applied  for  calculation  of  K  s ,  n  s,  and  AH  s  for  the  solid  acids  silica  gel 
[59],  ZSM-5  [43],  and  sulfated  zirconia  [88]. 

Results  and  Discussion 

Drying  the  Solids 

As  pointed  out  in  a  previous  report  [74]  the  dehydration  of  H3PW  to  the 
anhydrous  acid  is  complicated  by  decomposition  into  constituent  oxides.  We  have  shown 
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[74]  that  treatment  of  the  H3PW  at  160-200  °C  for  4  hours  under  vacuum  resuhs  in  the 
anhydrous  acid. 

The  drying  of  the  CsxHa.xPW  compounds  was  investigated  using  the  maximum 
total  heat  released  in  the  titration  with,  pyridine  as  a  criterion  of  dryness  [74]. 
Temperatures  of  200  and  300  °C  for  4  hours  under  vacuum  were  used  and  the  former  gave 
better  results,  Table  3-2.  All  cesium  salts  were  dried  at  200  °C  for  4  hours  under  vacuum. 

Cal-ad  Analysis  of  Solid  HjPW 

The  results  of  the  calorimetric  titration  and  adsorption  experiments  for  pyridine 
reacting  with  a  slurry  of  H3PW  in  cyclohexane  are  shown  in  Figure  3-1.  Cyclohexane 
was  chosen  as  the  solvent  because  it  is  a  non-specific  interacting  solvent  in  which  H3PW 
is  insoluble.  The  heat  evolved  from  the  reaction  indicates  that  not  all  of  the  1.04 
millimoles  of  protons  in  one  gram  of  H3PW  are  titrated.  After  0.24  mmol  of  pyridine  are 
added  (about  24%  of  the  moles  of  protons  in  a  gram  of  solid),  heat  evolution  practically 
stops  and  the  adsorption  of  pyridine  by  the  solid  also  reaches  a  plateau. 

Three  possible  equilibria  were  evaluated  and  the  results  are  presented  in  Table  3- 
3.  The  Cal-ad  data  were  analyzed  by  assuming  a  one  site  fit,  i.e.,  all  of  the  protons 
titrated  have  the  same  strength  and  equilibrium  constant.  The  resuhs  obtained  are:  AH  = 
-25.6  ±1.1  kcal  mol'',  K  =  3.5x10''  ±8.1x10^  and  n  =  0.185  ±0.005  mmol  g'.  This  fit 
shows  a  higher  sum  of  squares  and  higher  errors  than  a  fit  to  two  sets  of  non-equivalent 
protons.  The  two  site  fit  produces:  AHi  =  -32.7  ±0.3  kcal  mol"',  Ki=  3.7x10^  ±2.0xl0\ 
ni  =  0.079  ±0.002,  AH2  =  -19.6  ±4.8  kcal  mol"',  K2=  2.9x10^  ±1.0x10^  112=  0.159  ±0.056 
mmol  g"'.  The  deviations  in  the  experimental  and  calculated  heat  evolved  in  the  titration 
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Table  3-2.  Experimental  Heat  Evolved  (h'exp)  in  Cyclohexane  Upon  Addition  of  Pyridine" 
to  H3PW  and  CS2HPW  for  Different  Thermally  Treated  Solids^ 


— I  

h  exp  (cal) 


mmol  py  added 

H3PW  (162/4) 

H3PW  (200/4) 

CS2HPW  (200/4) 

CS2HPW  (300/4) 

0.0123 

0.365 

0.363 

0.299 

0.221 

0.0230 

0.693 

0.688 

0.540 

0.391 

0.0341 

1.029 

1.009 

0.754 

0.558 

0.0458 

1.381 

1.338 

0.862 

0.674 

0.0690 

2.076 

1.992 

1.018 

0.806 

0.0920 

2.704 

2.628 

1.152 

0.928 

0.1197 

3.279 

3.315 

1.294 

1.063 

0.1382 

3.613 

3.684 

1.390 

1.141 

0.1615 

3.946 

4.062 

1.467 

1.208 

0.1842 

4.198 

4.376 

1.538 

1.270 

0.2074 

4.421 

4.568 

1.598 

1.335 

0.2311 

4.618 

4.651 

1.658 

1.397 

a.  The  millimoles  of  pyridine  added  and  the  heat  evolved  values  are  summed  along  the 
titration. 

b.  Thermal  treatments  refer  to  temperature  (°C)  and  time  (hs)  which  the  solid  was  heated 
under  vacuum. 
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Figure  3-1 .  Isotherms  of  Calorimetry  and  Adsorption  for  Reaction  of  H3PW  with  Pyridine 
in  Cyclohexane. 
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Table  3-3.  Cal-ad  Results  for  Three  Types  of  Fit  for  H3PW  Reacting  with  Pyridine  in 
Cyclohexane. 


1  site  fit 

2  site  fit 

3  site  fit 

ni  (mmol  g"') 

0.186  ±0.005 

0.079  ±  0.002 

0.081  ±0.002 

K,  (M-') 

3.5x10'' ±8.1x10^ 

3.7x10^  ±2.0x10^ 

3.6x10^  ±2.3x10^ 

-AHi  (kcal  mof') 

25.6  ±  1.1 

32.7  ±0.3 

32.1  ±0.3 

n2  (mmol  g'') 

0.16  ±0.05 

0.08  ±  0.05 

K2  (M"') 

2.9x10^  ±  1.0x10^ 

3.0x10^  ±2.6x10^ 

-AH2  (kcal  mol"') 

19.6  ±4.8 

25.0  ±5.6 

n3  (mmol  g'') 

0.09  ±0.08 

K3  (M-') 

2.5x10^  ±  1.9x10^ 

-AH3  (kcal  mol"') 

15.1  ±  12.2 

U* 

1.60x10"' 

6.59x10"^ 

1.65x10"^ 

a.  The  quantity  U  is  the  sum  of  squares  of  residuals,  calculated  for  calorimetry  as: 

X  (h  exp'h  calc)  • 
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are  close  to  the  sum  of  the  squares  of  the  experimental  error  (1.65x10"^).  The  two  site 
model  also  gives  a  value  for  the  sum  of  the  sites  that  is  in  good  agreement  with  the 
estimate  of  the  moles  of  pyridine  reacted  from  the  calorimetric  exotherm.  There  is  a  four- 
fold improvement  in  the  sum  of  squares  for  the  three  site  fit  compared  to  the  two  but  the 
improvement  in  the  deviations  in  the  data  fit  are  not  significant  compared  to  experimental 
error.  Within  experimental  error  a  three  site  fit  cannot  be  distinguished  from  the  two  site 
fit.  This  leads  to  a  choice  of  the  two  site  fit  and  the  next  concern  is  with  the  assignment 
of  the  sites. 

The  number  of  surface  protons  reported  [89]  for  H3PW,  0.008  mmol  g"',  agrees 
with  our  estimate  from  its  surface  area  (5  m^  g'')  and  the  diameter  of  an  assumed 
spherical  Keggin  anion  (10  A).  This  is  an  order  of  magnitude  less  than  calculated  sitel. 
Thus,  in  addition  to  reacting  with  the  protons  on  the  surface,  pyridine  also  penetrates  into 
and  partially  reacts  with  those  in  the  bulk  solid.  This  result  is  consistent  with  reports  [90] 
that  the  absorption  of  polar  molecules  by  solid  H3PW  lead  to  catalytic  reactions  in  the 
bulk  of  the  crystal  as  well  as  on  the  surface.  The  solid  behaves  like  a  concentrated 
solution  and  this  behavior  is  described  as  a  pseudoliquid  phase  [11].  Non-polar 
molecules  are  not  absorbed. 

In  this  titration,  not  all  of  the  protons  react.  In  order  to  determine  if  incomplete 
reaction  is  caused  by  forming  a  layer  of  product  on  the  surface  that  could  not  be 
penetrated  by  pyridine,  a  sample  of  the  same  calcined  solid  was  ground  to  finer  particles 
inside  the  inert  atmosphere  box  and  titrated  with  pyridine.  An  identical  isotherm  resulted 
demonstrating  that  restricting  diffusion  from  a  surface  constraining  is  not  a  problem. 
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Pyridine  does  not  access  certain  protons  in  the  structure.  The  reaction  with 
pyridine  may  require  a  slight  lattice  expansion.  In  order  to  understand  how  an  expansion 
in  the  structure  produce  a  layered  material,  it  is  necessary  keep  a  clear  image  of  the 
complete  structure  of  H3PW.  Figure  3-2(a)  illustrates  the  unit  cell  of  H3PW.  This  unit 
cell  consists  of  a  total  of  two  anions,  8  anions  at  the  vertices  of  a  cube  shared  by  eight 
cubes  and  1  unshared  anion  in  the  center.  Six  of  the  protons  are  positioned  on  the  center 
of  faces  and  shared  by  two  units  cells.  Twelve  of  protons  bridge  Keggin  ions  and  are 
located  on  12  edges.  The  edge  protons  are  shared  by  four  unit  cells  giving  the  six  protons 
needed  to  charge  compensate  the  two  Keggin  anions  to  produce  the  molecular  formula 
H3PW.  The  quantity  of  protons  titrated  (24%),  suggests  that  on  average  one  layer  is 
neutralized  out  of  every  four.  In  a  given  layer,  four  of  the  protons  neutralized  in  the  unit 
cell  are  shared  with  four  unit  cells  and  the  face  proton  is  shared  by  two.  This  gives  rise  to 
a  2:1  ratio  of  different  type  protons  labeled  edge  and  face  in  Figure  3-2(b).  This  ratio 
corresponds  to  the  ratio  of  sites  calculated  from  Cal-ad  with  enthalpies  of  -32.7  ±0.3 
consisting  of  0.08  mmol  g''  and  those  with  enthalpies  of -19.6  ±4.8  kcal  mof'  consisting 
of  0.16  mmol  g''.  Thus  we  assigned  the  strongest  sites  to  the  face  protons.  It  is  important 
to  realize  that  all  protons  on  the  H3PW  structure  are  equivalent,  as  can  be  observed  in  the 
extended  structure.  The  two  sites  arise  when  pyridine  starts  to  react  with  the  protons  and 
therefore  it  creates  a  non  symmetric  environment  of  the  remained  protons. 

The  Cal-ad  results  with  pyridine  differ  from  microcalorimetry  of  gaseous 
ammonia  adsorption  [40-42].  With  ammonia  the  three  protons  on  H3PW  react  and  the 
enthalpy  decreases  from  -48  to  -39  kcal  mol"'  over  the  titration.  The  ammonia  enthalpies 
are  larger  for  they  contain  both  donor-acceptor  and  dispersion  contributions.  In  contrast 
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Edges  only     Faces  Only     Kegginlons  K=Comer Keggin ions 

H3PW-ABA  type,  bcc  K=Center  Keggin  ion 
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Figure  3-2.  Structure  of  H3PW12O40:  (a)  Simple  Unit  Cell;  (b)  Extended  Structure. 
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to  all  protons  reacting,  Cal-ad  reaction  of  H3PW  with  pyridine  indicate  that  only  24%  of 
the  available  protons  react  with  pyridine  and  these  have  different  strengths  (-32.7  and  - 
19.6  kcal  mol"').  The  differences  in  numbers  of  protons  reacted  is  attributed  to  the  small 
size  of  ammonia  which  enables  it  to  occupy  positions  occupied  by  water  in  the  hydrate 
without  appreciable  change  in  lattice  dimensions. 

Results  of  FTIR.  ^'P  MAS  NMR  and  XRD  for  H.PW 

FTIR  spectra  of  samples  produced  from  the  reaction  of  H3PW  with  pyridine 
contain  unchanged  characteristic  Keggin  bands.  The  region  between  1700-1450  cm"'  is 
informative  for  it  contains  bands  at  1540  cm"'  and  1485  cm"'  which  have  been  assigned  to 
pyridinium  ion  [91].  Coordinated  pyridine  or  hydrogen-bonded  pyridine  shows  a  band  at 
1612  cm"'  [91].  The  pyridinium  band  is  clearly  observed  (Figure  3-3)  for  H3PW,  but  only 
a  small  band  at  1612  cm"'  is  shown.  The  relative  ratio  based  upon  the  height  of  the  bands 
at  1612  cm''  to  that  at  1540  cm"'  is  0.2.  This  ratio  is  different  that  obtained  by  Cal-ad 
calculations  from  nj  over  ni  (2.0).  Thus,  site  2  can  not  be  considered  as  hydrogen-bonded 
pyridine  to  H3PW,  which  is  consistent  with  the  high  enthalpy  for  that  site  (AH2  =  -19.6 
kcal  mol'').  Both  protons,  assigned  as  site  1  and  2,  in  H3PW  are  Bronsted  type. 

The  ^'P  MAS  NMR  spectrum  of  H3PW  contains  a  single  sharp  peak  at  -12.05 
ppm  relative  to  an  external  85  wt.%  H3PO4  reference  (Figure  3-4a).  This  compares 
favorably  with  the  literature  reports  of  -10.9  ppm  [81],  -11.1  ppm  [52]  and  -12.4  ppm 
[92].  The  spectrum  of  the  product  of  the  pyridine  titration  gave  a  single  resonance  at  - 
16.53  ppm  (Figure  3-4b).  This  is  a  surprising  resuh  because  only  24%  of  the  protons  in 
H3PW  have  reacted  but  is  not  unexpected  in  view  of  the  mobility  of  the  protons.  The 
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Figure  3-3.  FTIR  Spectra  of  H3PW  (a)  and  CS3PW  (b)  after  Reaction  with  Pyridine 
Cyclohexane. 
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Figure  3-4.  ^'P  MAS  NMR  Spectra  of:  (a)  H3PW  (5  =  -12.05  ppm);  (b)  1.04  mmol 
H3PW-O.24  mmol  Pyridine  (6  =  -16.53  ppm). 
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single  peak  suggests  that  the  remaining  protons  not  titrated  are  not  symmetrically  located 
between  the  unit  cells  of  the  average  four  layer  units  but  are  displaced  in  the  bridges 
toward  the  anionic  layer.  In  this  way  the  negative  charge  formed  by  reaction  of  the 
protons  with  pyridine  in  a  given  layer  is  distributed  over  the  structure  making  the  electron 
density  on  all  the  anions  similar.  Spectra  of  ^'P  NMR  for  1 .04  mmol  of  H3PW  reacting 
with  0.08  and  0.12  mmol  of  pyridine  also  show  a  single  peak  at  -17.49  ppm  suggesting 
the  protons  are  mobile  in  the  structure  (Figure  3-5). 

The  XRD  spectrum  of  H3PW  gives  a  pattern  similar  to  the  reported  in  the 
literature  [78]  (Figure  3-6).  The  main  peaks,  related  to  the  strongest  reflections,  are 
present  at  29  =  10.40,  20.90,  and  25.40.  The  broader  peaks  indicates  loss  of  crystallinity, 
as  reported  for  H3PW  dried  at  different  conditions  [45,  93].  The  XRD  pattern  of  the 
product  of  the  reaction  of  1.04  mmol  (1  gram)  of  H3PW  with  0.24  mmol  of  pyridine 
shows  various  new  reflections  (Figure  3-7).  The  most  intense  peak  now  occurs  at  20  = 
25.57.  The  peak  at  29  =  10.40  split  up  in  two  less  intense  shoulders  at  20  =  10.30  and 
10.52.  Two  new  reflections  at  29  =  9.27  and  9.98  are  present  in  the  powder  pattern. 
These  results  confirms  formation  of  a  different  compound.  No  extensive  change  was 
observed  in  the  XRD  pattern  when  1.04  mmol  of  H3PW  was  reacted  with  0.08  mmol  of 
pyridine  (site  1  of  Cal-ad),  except  that  a  more  crystalline  spectrum  was  obtained  Figure  3- 
8).  The  main  peaks  with  the  intensities  in  the  XRD  spectra  of  those  compounds  are 
presented  in  Table  3-4. 
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Figure  3-5.  ^'P  MAS  NMR  Spectra  of:  (a)  1.04  mmol  H3PW-O.O8  mmol  Pyridine  (6  = 
17.49  ppm);  (b)  1.04  mmol  H3PW-O.I2  mmol  Pyridine  (5  =  -17.49  ppm). 
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Figure  3-6.  XRD  Pattern  of  H3PW  (dried  at  162  °C/  4  hs/  vacuum). 
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Figure  3-7.  XRD  Pattern  of  1 .04  mmol  H3PW-O.24  mmol  Pyridine. 
600 


500 


400 


c  300 
3 
o 
o 


200 


100 


new  sample  3:  hjPW+n^Py 


— 1  ■  1  .  1  .  1  .  1  ,  r— 

10         15         20         25         30  35 

2-Theta 


Figure  3-8.  XRD  Pattern  of  1.04  mmol  H3PW-O.O8  mmol  Pyridine. 
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Table  3-4.  XRD  Patterns  of  H3PW  and  Products  of  Reaction  with  Pyridine. 


2G  Angle  (degree) 

d-spacing  (A) 

Intensity  (V) 

H3PW  (162  C/  4  hs/  vacuum) 

10.40 

8.51 

61.5 

14.90 

5.95 

6.3 

18.10 

4.90 

7.0 

on  on 

4.25 

10.3 

23.20 

3.83 

7.6 

25.40 

3.51 

32.7 

29.00 

3.08 

7.5 

33.50 

2.68 

7.6 

1.04  mmol  H3PW-O.24  mmol  Pyridine 

9.27 

9.53 

8.7 

9.98 

8.86 

19.7 

10.30 

8.58 

27.3 

10.40 

8.51 

60.0 

10.52 

8.41 

42.9 

1 9  1  "J 

18. 1 J 

A  OQ 

4.oy 

36.2 

21.00 

4.23 

9.41 

23.30 

3.82 

8.2 

25.57 

3.48 

142.5 

29.20 

3.06 

20.6 

1.04  mmol  H3PW-O.O8  mmol  Pyridine 

10.40 

8.51 

47.8 

1  8  1  n 

18.  lU 

A  on 

1  n  1 
19.1 

20.90 

4.25 

7.4 

23.30 

3.82 

10.4 

25.60 

3.48 

84.9 

29.20 

3.06 

20.77 
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Cal-ad  Acidity  of  Cs^H^^yPW  Compounds 

Calorimetric  and  adsorption  isotherm  curves  for  the  reaction  of  CsxHa.xPW  sahs 
with  pyridine  in  cyclohexane  are  shown  in  Figures  3-9  and  3-10  respectively.  The  results 
of  the  Cal-ad  analysis  are  reported  in  Table  3-5.  Fits  of  the  combined  calorimetric  and 
adsorption  data  on  all  systems  were  attempted  for  both  one  and  two  processes.  For  all 
solids,  the  one  site  fit  gave  large  errors  in  the  resulting  K,  n,  and  AH  values  as  well  as  the 
calculated  heats  and  moles  adsorbed.  Meaningful  parameters  with  improved  errors  and  a 
decrease  in  the  sum  of  the  squares  of  the  fit  residuals  result  when  two  sites  are  employed. 
FTIR  spectra  of  cesium  salts,  except  CS3PW,  after  reaction  with  pyridine  also  supports 
this  choice,  since  they  indicate  both  Bronsted  and  hydrogen-bonded  pyridine  bands  (vide 
infra).  The  second  site  present  on  the  analysis  of  CS3PW  is  determined  because  formation 
of  pores  in  this  compound,  which  is  responsible  for  pyridine  physisorbed.  It  is  important 
also  to  keep  in  mind  the  results  of  cesium  analysis  (Table  3-1)  when  interpreting  the 
numbers  of  sites  and  moles  of  protons  reacted  obtained  by  Cal-ad. 

The  CsHiPW  titration  with  pyridine  involves  about  0.10  mmol  g"'  of  acid  sites 
(16%  of  the  total  0.66  mmol  of  protons  present  in  one  gram  of  solid).  With  CsiHPW, 
about  0.12  mmol  g"'  react  with  pyridine  (38%  of  the  total  0.32  mmol  of  protons  present  in 
one  gram  of  solid).  There  are  0.16  mmol  of  protons  in  one  gram  of  Cs2.5H0.5PW,  and  they 
all  react  with  pyridine.  In  CS3PW,  there  are  0.09  mmol  g''  of  the  residual  protons  in  the 
solid.  The  higher  number  of  second  sites  formed  for  CssPW  is  attributed  to  physisorbed 
pyridine.  CS3PW  has  the  same  structure  of  H3PW  [27],  but  with  a  higher  surface  area. 
The  presence  of  micropores  in  CS3PW  allows  some  pyridine  to  be  adsorbed  in  the  solid. 
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Figure  3-9.  Calorimetric  Isotherms  of  CsxHs.xPW  Salts  Reacted  with  Pyridine  in 
Cyclohexane. 
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2.5e-4 


Figure  3-10.  Adsorption  Isotherms  of  CsxHa.xPW  Salts  Reacted  with  Pyridine 
Cyclohexane. 
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Table  3-5.  Cal-ad  Results  for  Reactions  of  CsxHa.xPW  Compounds  with  Pyridine  in 
Cyclohexane. 


L.S3r  W 

til  /^mmnl  o~^\ 

Ll\   lllllliUl  ^  1 

n  099  +0  001 

n  097  4-0  001 

0  0'\'J  4-0  oo*; 

U.U4j  ±U.UUi 

1  -Sxio^  i"?  1x1  n'' 

1  9vio^  +1  Rvin^ 

S  OyI  0^  4A  7y1  O"* 

7  1  -vl  0**  4-1  Ox/ 1  a3 
Z.IXIU  ±1.ZX1U 

ATT/"!  1 

-AHi(kcal  mol  ) 

20.0  ±0.3 

28.3  ±0.6 

27.4  ±0.5 

12.7  ±0.5 

n2  (mmol  g'') 

0.08  ±0.02 

0.09  ±0.03 

0.14  ±0.09 

0.15  ±0.15 

K7  CM'') 

3  6x10^  ±9  9x10^ 

1  5x10^  ±S  7x10^ 

4  6x  1 0"*  ±3  1 X 1 0" 

1  1 V 1 0^  49  9v  1 

-AH2(kcal  mol"') 

9.3  ±2.3 

15.1  ±2.7 

6.4  ±4.6 

1.4  ±1.8 

Total  " 

0.66 

0.32 

0.16 

0.09 

(mmol  g"') 

reacted  " 

0.10 

0.12 

0.19 

0.19 

(mmol  g'') 

a.  This  amount  is  based  on  the  stoichiometric  formula  of  each  compound,  except  CS3PW 
where  the  actual  formula  was  used. 

b.  This  amount  is  based  upon  the  number  of  sites  (ni  +  n2)  obtained  by  Cal-ad. 
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The  residual  protons  of  CS3PW  are  mainly  on  the  surface,  which  are  homogeneously 
distributed  after  thermal  treatment  of  the  solid. 

The  enthalpies  are  informative  when  compared  to  the  enthalpies  of  reaction  of 
H3PW  with  pyridine.  For  CS3PW,  the  magnitude  of  the  AHi  (-12.7  kcal  mol"')  is 
comparable  to  the  hydrogen  bond  site  for  reaction  of  pyridine  with  silica  gel  [59]  (AHi  =  - 
12.6  kcal  mol"').  The  IR  spectrum  (Figure  3-3)  does  not  contain  the  band  assigned  to 
pyridinium  ion  but  only  the  hydrogen-bonded  adduct  band  is  observed  at  1614  cm"'  [91], 
confirming  the  nature  of  this  site  for  CS3PW. 

For  CSH2PW,  the  AHi  =  -20.8  kcal  mof'  is  close  to  the  value  of  AH2  obtained  for 
H3PW  (-19.6  kcal  mol"').  The  FTIR  spectrum  of  this  solid  (vide  infra)  after  reaction  with 
pyridine  shows  a  very  small  peak  at  1 540  cm"'  (pyridinium  ion),  which  is  assigned  to  the 
strongest  site.  A  second  band  at  1614  cm"'  corresponds  to  hydrogen-bonded  pyridine  to 
the  solid  CSH2PW.  Compared  to  FTIR  spectrum  of  H3PW  reacted  with  pyridine  (Figure 
3-3  ),  CsHaPW  shows  a  larger  band  for  hydrogen-bonded  pyridine.  As  discussed  before, 
this  pyridinium  band  is  consistent  to  be  assigned  to  the  strongest  site.  Having  the  same 
enthalpy  of  site  2  for  H3PW,  this  CSH2PW  compound  has  fewer  strong  sites  responsible 
for  the  higher  enthalpy  obtained.  This  is  compatible  with  the  Cal-ad  resuks, 
demonstrating  site  1  for  CsHaPW  has  the  same  nature  of  second  site  for  H3PW. 

The  enthalpy  result  for  site  1  on  CsHaPW  indicates  that  the  site  1  protons  in 
H3PW  were  completely  neutralized  by  cesium  and  there  is  no  accessible  H3PW  in  the 
solid.  Only  1 6%  of  the  available  protons  in  the  solid  are  titrated  which  is  lower  than  the 
24%  titrated  for  H3PW.  Considering  that  0.34  mmol  g"'  of  solid  are  cesium  atoms  and 
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0.10  mmol  g"'  are  the  total  protons  reacted  with  pyridine,  0.44  mmol  g"'  of  the  total  1.04 
mmol  g"'  reacted  from  the  H3PW.  These  42%  of  the  original  1.04  mmol  of  protons  from 
H3PW  correspond  to  a  ratio  of  4:1  between  the  protons  positioned  at  the  edges  and  on  the 
faces.  That  ratio  is  close  to  the  one  obtained  by  Cal-ad  from  site  2  to  site  1  (3.6).  This 
suggests  that  reaction  with  pyridine,  on  average,  occurs  in  one  out  of  every  two  layers  in 
the  solid.  Coupled  with  the  ^'P  NMR  results  (vide  infra)  which  shows  no  free  H3PW,  we 
conclude  that  the  reacted  layer  is  mainly  cesium  atoms  which  combined  with  the  low 
surface  area  prevents  access  to  most  of  the  protons  in  the  solid.  The  lower  enthalpy  of 
reaction  for  site  2  (9.3  kcal  mol'')  is  attributed  to  hydrogen-bonded  pyridine,  as  confirmed 
by  the  presence  of  a  large  band  at  1614  cm''  in  the  FTIR  spectrum. 

The  enthalpies  of  reaction  of  pyridine  with  CS2HPW  together  with  the  quantity  of 
acidic  protons  (ni  =  0.03  mmol  g'')  are  significantly  larger  than  CSH2PW.  The  larger  ni 
is  attributed  to  the  tertiary  structure  which  leads  to  formation  of  pores  and  therefore 
higher  surface  area.  The  strongest  acid  site  for  CS2HPW  (AHi  =  -28.3  kcal  mol"')  is  a 
weaker  acid  than  the  first  site  of  H3PW.  This  can  be  qualitatively  evidenced  by  a  smaller 
band  at  1540  cm"'  for  CS2HPW  than  for  H3PW.  However,  this  band  is  a  little  larger  than 
the  one  obtained  by  CSH2PW,  which  is  concordant  with  the  higher  enthalpy  of  CS2HPW 
(vide  infra). 

The  ^'P  NMR  (vide  infra)  shows  only  one  peak  (-14.95  ppm)  for  CS2HPW  and  no 
signal  at  -12.05  ppm  (free  H3PW).  A  previous  report  [81]  based  on  ^'P  NMR  results, 
pointed  out  that  CS2HPW  is  a  mixture  of  four  different  species  of  CsxH3.xPW  (x=  0,  1,2, 
and  3).  According  to  our  results,  this  compound  is  a  homogeneous  solid  with  the  protons 
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being  homogeneously  distributed  in  the  structure.  The  preparation  from  aqueous  solution 
might  generate  a  series  of  equilibria  in  solution  involving  up  to  four  species,  but  after 
thermal  treatment  ^'P  NMR  shows  the  cesium  atoms  are  homogeneously  distributed  in 
the  solid  CS2HPW.  Considering  that  thermal  treatment  for  this  solid  was  the  same  that 
one  reported  (200  C  under  vacuum)  [81],  the  differences  demonstrates  that  structurally 
different  solids  were  obtained.  This  fact  is  confirmed  by  the  XRD  spectrum  (vide  infra) 
which  shows  small  differences  from  the  one  reported  [80]. 

The  magnitude  of  the  enthalpy  values  of  AHi  (-28.3  kcal  mol"')  and  AH2  (-15.1 
kcal  mol"')  for  CS2HPW  show  this  solid  is  a  weaker  acid  than  H3PW  and  is  not  a  simple 
mechanical  mixture  of  the  H3PW  and  the  other  cesium  salts.  The  IR  spectrum  (vide 
infra)  also  confirms  the  presence  of  a  Bronsted  site  (1540  cm"')  and  a  weaker  hydrogen- 
bonded  site  (1612  cm"'). 

The  distribution  of  the  protons  in  the  structure  is  considered  from  the  0.68  mmol 
of  cesium  per  gram  of  CS2HPW  and  the  total  protons  reacted  (0.12  mmol  g"'  of  solid). 
Thus  0.80  mmol  of  protons  per  gram  reacted  from  1.04  mmol  g"'  of  H3PW.  The  77%  of 
the  original  protons  of  H3P W  correspond  to  a  ratio  of  3 : 1  between  protons  positioned  at 
the  edges  and  on  the  faces.  That  ratio  is  close  to  that  obtained  by  Cal-ad  from  site  2  to 
site  1  (3.3).  This  suggests  that  on  average,  reaction  of  CsaHPW  with  pyridine  occurs  in 
four  out  of  every  five  layers  in  the  solid. 

It  is  observed  that  the  strong  acid  site  in  all  these  solids  is  only  a  fraction  of  the 
protons  available,  but  this  fraction  increases  with  the  amount  of  cesium  atoms.  The 
strength  decreases  in  the  order  H3PW  >  CS2HPW  >  CSH2PW.  This  order  is  not  expected 
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on  the  basis  of  charge  on  the  conjugate  base  and  suggests  that  the  lattice  energy  changes 
that  accompany  reaction  with  pyridine  make  a  significant  contribution  to  the  enthalpy 
measured. 

Solid  CS2  5H0.5PW  is  reported  to  be  a  mixture  of  CS2HPW  and  CS3PW  [81].  Our 
results  from  ^'P  NMR  pointed  out  that  only  one  type  of  compound  is  present  (vide  infra). 
The  increased  amount  of  acidic  sites  going  from  CsHaPW  to  Cs2.5H0.5PW  indicates  that 
exposure  of  protons  are  markedly  dependent  of  the  surface  area  for  these  cesium 
heteropoly  compounds.  The  formation  of  pores  are  also  well  evidenced  from  the 
adsorption  of  pyridine  on  Cs2  5H0  5PW  (Figure  3-10).  As  the  surface  area  increases,  more 
protons  can  be  accessed  and  titrated. 

FTIR  spectrum  of  Cs2  5H0  5PW  (vide  infra)  shows  a  larger  band  for  the  pyridinium 
ion  (1540  cm"')  compared  to  the  CS2HPW  and  CSH2PW,  in  this  order  respectively.  This 
agrees  to  the  Cal-ad  results  where  the  Cs2  5H0.5PW  is  as  sfrong  as  CS2HPW,  and  more 
protons  are  accessible.  Our  Cal-ad  results  suggests  that  every  single  layer  of  the 
Cs2  5H0.5PW  salt  reacts  with  pyridine.  This  can  be  obtained  considering  that  0.87  mmol 
of  cesium  per  gram  of  Cs2.5H0.5PW  plus  the  total  protons  reacted  (about  0.17  mmol  g''  of 
solid)  leads  to  1 .04  mmol  of  protons  per  gram  reacted  from  1 .04  mmol  g"'  of  H3PW.  The 
100%  of  the  original  1.04  mmol  of  protons  from  H3PW  correspond  to  a  ratio  of  3:1 
between  protons  positioned  at  the  edges  and  on  the  faces.  That  ratio  is  close  to  that 
obtained  by  Cal-ad  from  site  2  to  site  1  (2.6).  Thus,  all  layers  of  Cs2  5H0  5PW  have  been 
titrated  with  pyridine. 
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The  magnitude  of  the  AHi  for  Cs2  5H0.5PW  (-27.5  kcal  mol"')  as  well  as  the  higher 
number  of  acidic  protons  (m  =  0.05  mmol  g"'),  compared  to  only  0.008  mmol  g"'  for 
surface  protons  on  H3PW,  show  why  this  compound  is  more  effective  in  some  alkylation 
reactions  than  H3PW.  With  presence  of  pores  in  its  structure  and  higher  surface  area  than 
H3PW,  Cs2  5H0.5PW  have  more  surface  protons  available  to  catalytic  reactions. 

Results  of  FTIR,  XRD  and  ^'P  MAS  NMR  of  CsvH^.yPW  Compounds 

FTIR  spectra  of  CsxHs-xPW  salts  were  taken  before  and  after  reaction  of  the 
compounds  with  pyridine.  A  comparative  analysis,  presented  in  Table  3-6,  concludes  the 
Keggin  anion  is  practically  not  altered  by  substitution  of  by  Cs"^,  as  observed  for 
and  Al^*  salts  [51].  Reaction  with  pyridine  does  not  change  the  characteristic  Keggin 
bands.  However,  the  region  between  1700-1450  cm"'  shows  the  pyridinium  ion  through 
the  presence  of  the  bands  at  1540  cm"'  and  1485  cm"'  [91].  Coordinated  or  hydrogen- 
bonded  pyridine  shows  a  band  at  1612  cm"'  [91].  Those  bands  were  observed  (Figure  3- 
1 1)  for  all  cesium  salts  except  CS3PW  which  presents  only  a  wide  band  at  1612  cm"'.  The 
ratio,  based  on  the  relative  height  of  the  bands,  between  the  hydrogen-bonded  pyridine 
band  (1612  cm"')  to  the  pyridinium  band  (1540  cm"')  are  3.9  for  CsHaPW,  3.1  for 
CS2HPW,  and  2.4  for  Cs2.5H0.5PW.  Compared  to  Cal-ad  analysis  these  ratio  of  the  bands 
support  the  choice  of  two  site  fit  for  the  cesium  salts,  since  the  ratio  between  n2  to  ni  are 
3.6  for  CSH2PW,  3.3  for  CS2HPW,  and  2.6  for  CS2  5H0  5PW.  Based  on  these  results  the 
nature  of  second  site  for  CSH2PW,  CS2HPW  and  Cs2  5H0  5PW  are  assigned  as  hydrogen- 
bonded  pyridine  in  the  solid,  while  the  strongest  site  1  is  a  Bronsted  type  site. 
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Table  3-6.  Comparative  Infrared  Bands  (cm'')  of  the  Keggin  Anion  for  H3PW  and 
CS2HPW. 


Group 

T  T    TfcH  78 

H3PW 

H3PW  dry 

H3PWI6H2O 

HsPW-Py 

CS2HPW' 

Cs2HPW-Py 

P-0 

1081 

1082 

1076 

1080 

1078 

1078 

W-Ot 

984 

991 

978 

978 

984 

984 

W-Oc-W 

892 

895 

916 

891 

889 

889 

W-Oe-W 

794 

801 

734 

806 

799 

799 

a.  Values  obtained  for  H3PW6  H2O  from  reference  [51]. 

b.  Values  obtained  for  H3PW  dried  at  162  °C/4  hs/vacuum  from  reference  [74]. 

c.  CS2HPW  dried  at  200  C/4  hs/vacuum.  The  other  cesium  salts  have  the  same  bands. 
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Figure  3-11.  FTIR  Spectra  of  CSH2PW  (a),  CsjHPW  (b)  and  Cs2.5H0.5PW  (c)  after 
Reaction  with  Pyridine  in  Cyclohexane. 
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The  XRD  patterns  for  the  cesium  salts  were  recorded.  The  results  are  presented  in 
Figure  3-12(a-d).  It  is  observed  that  the  main  peaks  present  on  the  free  H3PW  are  shifted 
toward  higher  angles  in  CsHaPW  and  CsiHPW.  This  result  is  consistent  with  the  smaller 
cubic  lattice  parameters  obtained  for  these  compounds  in  the  literature  [80].  However, 
Cs2.5H0.5PW  and  CS3PW  shows  shifts  to  lower  angles,  compared  to  H3PW,  which 
indicates  these  prepared  compoimds  may  be  different  from  the  ones  reported  in  the 
literature  [80].  In  fact,  no  XRD  pattern  is  reported  for  Cs2.5H0.5PW  in  the  literature.  The 
spectra  display  only  small  differences  going  from  CSH2PW  to  CS3PW,  but  a  carefiilly 
analysis  of  the  reflections  in  Table  3-7  demonstrates  they  are  distinct,  although  they 
probably  have  similar  structures. 

Spectra  of  ^'P  MAS  NMR  were  taken  on  all  cesium  derivatives.  The  results  are 
presented  in  Figure  3-13.  Contrary  to  the  literature  [81],  the  spectra  show  only  a  single 
peak  for  each  cesium  containing  compound.  There  is  an  approximately  linear  shift  to 
lower  field  in  the  ■'^P  NMR  spectra  going  from  H3PW  to  CS3PW.  Exchange  of  protons 
for  cesium  atoms  causes  increased  electronic  density  at  the  phosphorous  core  because  the 
protons  are  less  mobile  in  these  structures.  The  lower  mobility  of  the  protons  in  the 
cesium  containing  compounds  causes  the  observed  shielding  shift.  These  single  peaks 
suggest  that  the  protons  are  located  in  random  positions  in  the  structure  of  each  of  the 
cesium  derivatives.  This  result  conflicts  with  another  report  [81]  but  is  consistent  with  all 
other  measurements  performed  in  these  solids.  No  evidence  of  mixtures  of  the  cesium 
compounds  with  free  H3PW  could  be  assigned  in  this  report. 
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Figure  3-12.  XRD  Pattern  of:  (a)  CsHjPW;  (b)  CsjHPW. 


Figure  3-12.  XRD  Pattern  of:  (c)  Cs2.5H0.5PW;  (d)  CS3PW. 
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Table  3-7.  XRD  Patterns  of  CsHzPW,  CS2HPW,  CS2  5H0  5PW,  and  CS3PW. 


26  Angle  (degree) 

d-spacing  (A) 

Intensity  (V) 

CsHzPW 

1  n  Ac\ 

O.D  1 

11 A 

18.20 

4.87 

14.6 

21.00 

4.23 

5.5 

23.30 

3.82 

7.7 

25.60 

3.48 

56.4 

29.30 

3.05 

16.4 

CS2HPW 

lU.oU 

Q  1^ 
O.J  J 

40. U 

15.10 

5.87 

7.8 

18.60 

4.77 

7.6 

21.10 

4.21 

12.3 

25.80 

3.45 

23.8 

29.60 

3.02 

3.8 

105 


-7.  continued 
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Figure  3-13.  ^'P  MAS  NMR  Spectra  of:  (a)  CsH2PW(5  =  -12.53  ppm);  (b)  CsiHPW  (6 
-14.95  ppm). 
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Figure  3-13.  ^'P  MAS  NMR  Spectra  of:  (c)  CS2  5H0.5PW  (5  =  -15.51  ppm);  (d)  CsjPW  (6 
=  -17.01  ppm). 
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Comparison  of  Acidity  of  H^PW  in  Homogeneous  and  Heterogeneous  System 

Acidity  of  H3PW  was  measured  in  solution  of  CH3CN  using  calorimetry  [74].  In 

homogeneous  solution  of  a  solvent  of  equal  or  lower  basicity  than  acetonitrile,  H3PW 

behaves  as  a  triprotic  acid  with  the  first  proton  stronger,  but  only  partially  dissociated. 

According  to  ECW  analysis  [74],  the  first  proton  has  an  enthalpy  of  -30.4  kcal  mol"' 

when  it  reacts  with  pyridine  (corrected  for  the  enthalpy  of  displacement  of  CH3CN) .  The 

value  obtained  in  this  work  (-32.7  kcal  mol"')  is  in  good  agreement  with  the  former  study. 

This  demonstrate  that  the  use  of  enthalpies  to  compare  neutral  donor-acceptor  reactions 

and  ionic  products  are  possible  in  low  polarity  solvents,  like  acetonitrile,  as  well  as  non 

polar  solvents,  like  cyclohexane.  This  is  only  possible  if  the  contribution  of  solvation  of 

the  polar  solvent  can  be  factor  out.  Using  the  ECW  model  [74]  an  estimate  of  the  net 

solvation  energies  of  products  and  reactants  together  with  the  endothermic  enthalpy  of 

dissociation  of  the  anion  from  (CH3CN)3H3PW  was  obtained. 

As  pointed  out  [74],  in  solution  of  acetonitrile  an  equilibrium  mixture  between  the 

completely  dissociated  ion  pair  and  hydrogen  bonding  species  is  possible  for  the  H3PW 

equilibria,  especially  for  the  second  and  third  steps  of  its  dissociation: 

Bsoi  +  H-Xsoi  — ^  B-H-Xsoi         B-H^-X-soi         BH^oi  +  X'soi  (3.4) 

(a)  (b)  (c) 

The  solvent  may  influence  the  association,  proton-transfer,  and  dissociation  steps.  The 

main  factor  which  determines  the  position  of  a  Bronsted  acid-base  equilibrium  in 

Equation  3.4  is  the  differential  solvation  of  the  covalent  and  the  ionic  hydrogen-bonded 

complexes  (a)  and  (b).  For  H3PW  in  CH3CN,  second  and  third  proton  dissociations  have 

shown  possible  mixtures  between  species  (b)  and  (c)  [74]. 
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In  cyclohexane  much  less  complications  are  related  to  the  solvent.  H3PW  is  not 
soluble  in  this  solvent.  Cyclohexane  is  a  non  specific  interacting  solvent.  It  has  a  low 
dielectric  constant,  and  therefore  is  expected  formation  of  only  ion  pairs  in  solution 
(products  b  and  c  from  Equation  3.4).  H3PW  is  a  strong  Bronsted  solid  acid  and 
formation  of  ion-pair  product  in  cyclohexane  is  expected,  where  the  proton  is  completely 
transferred  to  the  base  (product  b  in  Equation  3.4). 

Comparison  of  the  second  enthalpy  in  cyclohexane  (AH2  =  -19.6  kcal  mol'^)  with 
that  obtained  in  CH3CN  (AH2  =  -21.2  kcal  mol''),  considering  the  displacement  of  the 
anion  (CH3CN)2HPW^-  to  be  equal  to  -9.4  kcal  mol"'  (the  same  order  of  magnitude  of  the 
dissociation  of  (CH3CN)2H2PW")  shows  a  close  agreement  between  enthalpy 
measurements  in  solution  and  in  the  solid.  The  enthalpy  of  the  second  site  in  solid  H3PW 
is  considerably  reduced  over  that  of  the  first.  This  could  result  from  an  endothermic 
contribution  from  expanding  the  lattice  of  H3P  W.  Titration  of  the  protons  on  the  edges 
of  the  cube  expand  the  structure,  while  the  strong  acid  site  on  the  face  of  the  unit  cell  is 
accessed  in  the  solid  structure.  In  solution,  the  lower  enthalpy  is  attributed  to  the  higher 
energy  bonding  of  the  proton  to  the  polyanion  resulting  from  the  dissociation  of  the  first 
proton. 

Conclusions 

This  work  provides  thermodynamic  constants  for  H3PW  and  CsxH3.xPW  (x  =  1,  2, 
2.5,  and  3)  reacting  with  pyridine  in  cyclohexane  solution.  The  Cal-ad  results  indicate 
that  the  protons  of  H3PW  in  the  solid  behaves  like  in  solution,  with  all  having  different 
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strengths.  This  achievement  disagrees  with  previous  claims  by  gas-solid  calorimetry  that 
all  protons  on  solid  H3PW  are  equivalent.  Similar  to  solutions  of  low  basicity  solvents 
such  as  CH3CN,  where  H3PW  behaves  as  a  strong  monoprotic  acid,  in  solid  state  it  can 
donate  protons  with  different  levels  of  acidity.  Restrictions  of  size  and  polarity  of  the 
donor  molecule  and  polarity  of  medium  can  influence  the  effectiveness  of  the  proton 
transfer.  Two  kinds  of  protons  with  different  strengths  are  assigned  for  solid  H3PW, 
according  to  its  cubic  unit  cell:  face  protons  (site  1)  and  edge  protons  (site  2).  The 
enthalpies  for  reaction  of  the  most  acidic  proton  (site  1)  with  pyridine  (AHi  =  -32.7  kcal 
mol"')  and  site  2  (AH2  =  -19.6  kcal  mol"')  are  in  good  agreement  with  the  value  obtained 
in  solution  of  acetonitrile  (AHi  =  -30.4  kcal  mol  '  and  estimated  AH2  =  -21.2  kcal  mol'') 
respectively. 

The  incomplete  neutralization  of  the  total  available  protons  of  H3PW  with 
pyridine  in  cyclohexane  is  explained  by  formation  of  a  layered  product.  On  average,  one 
out  of  every  four  layers  of  H3PW  is  neutralized  out  by  pyridine.  This  result  is  confirmed 
by  Cal-ad  through  the  ratio  of  the  number  of  sites  for  site  2  (assigned  to  protons  on  the 
edge)  to  site  1  (assigned  to  protons  on  the  face). 

Measurements  of  enthalpy  for  the  strongest  protons  in  each  compound  lead  to  the 
following  order  of  acidity:  H3PW  >  CS2HPW  =  Cs2  5H0  5PW  >  CSH2PW  >  CS3PW.  The 
increased  acidity  with  the  surface  area  in  the  cesium  series  confirms  Cs2.5H0.5PW  as  a 
strong  catalyst,  having  an  acidity  comparable  to  CS2HPW.  CS3PW  presents  only  residual 
protons  and  the  concentration  of  protons  becomes  very  low  and  so  does  the  acidity  as 
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measured  by  Cal-ad.  The  CSH2PW  compound  has  a  smaller  surface  area  than  H3PW 
which  explains  its  comparatively  lower  acidity. 

The  ^'P  MAS  NMR,  XRD,  FTIR  and  Cal-ad  results  confirm  that  CsxHj.xPW 
compounds  are  not  simple  mixtures  of  H3PW  with  the  exchanged  cesium  derivatives. 
The  positions  of  proton  and  cesium  atoms  in  the  salt  structures  are  random  and 
homogeneously  distributed.  Based  on  the  Cal-ad  results  for  the  relative  amount  of  sites  1 
and  2  in  the  cesium  derivatives  an  attempt  to  distribute  the  proton  and  cesimn  atoms  in 
the  structure  is  proposed. 

Reaction  with  pyridine  of  the  cesium  derivatives  shows  incomplete  neutralization 
of  the  total  available  protons  for  each  compound,  except  for  CS2.5H0  5PW  and  CS3PW. 
For  Cs2.5H0.5PW,  all  protons  available  react  with  pyridine,  which  concludes  that  every 
single  layer  in  the  solid  is  titrated.  CS3PW  has  only  residual  protons  and  all  react  with 
pyridine.  For  CsHaPW,  one  out  of  every  two  layers  reacts  with  pyridine,  while  for 
CS2HPW  four  out  of  every  five  layers  reacts  with  pyridine.  This  proposition  is  supported 
by  the  Cal-ad  results  as  well  as  the  other  spectroscopic  methods  used. 


CHAPTER  4 

ACIDITY  OF  H3PW12O40  SUPPORTED  ON  SILICA,  ALUMINA,  AND  CARBON 

Introduction 

Keggin-type  heteropoly  acids  (HP As)  have  been  extensively  studied  because  they 
possess  strong  acidic  and  redox  properties  [32].  The  different  acidity  and  relatively  high 
thermal  stability  of  the  heteropoly  acids  allows  them  to  be  used  in  several  acid-catalyzed 
reactions  in  homogeneous  solution,  liquid-solid  and  gas-solid  heterogeneous  reactions.  It 
has  been  demonstrated  that  HP  As  are  much  more  effective  in  acid  catalysis  than  protonic 
mineral  and  organic  acids  such  as  sulfuric  and  p-toluenesulfonic  acids  [31]. 

HP  As  have  been  supported  on  different  solids  to  increase  the  surface  to  find  a 
better  catalyst  for  reactions  of  non  polar  molecules  with  heteropoly  compounds.  It  is 
known  that  HPAs  can  absorb  polar  molecules  into  their  bulk  leading  to  a  pseudoliquid- 
phase  type  of  catalysis.  Non-polar  molecules  react  only  on  the  surface  or  between  surface 
layers  of  the  crystals  [11].  Because  the  surface  area  of  HPAs  are  usually  very  low, 
increasing  the  surface  area  of  supported  HPAs  could  produce  a  catalyst  with  higher 
activity  for  reactions  involving  non  polar  substrates. 

Silica  gel  is  the  most  often  used  carrier  for  HPAs.  Si02-H3PW  supported  has 
been  characterized  by  XRD  and  TPD  of  ammonia  [94],  microcalorimetry  with  ammonia 
adsorption  [41,  95],  and  NMR  [92,  96].  No  distinct  diffraction  peaks  appear  on  the  XRD 
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spectrum,  but  when  loading  is  higher  than  20  wt.%  a  broad  peak  occurs  at  20  =  10°, 
characteristic  of  H3PW.  Mechanical  mixtures  of  H3PW  and  Si02  show  this  peak  with 
amounts  of  heteropoly  acid  as  low  as  5  wt.%  [94].   According  to  TPD  of  ammonia 
absorption  [94],  the  strength  of  Si02-H3PW  supported  is  independent  of  amount  of  the 
heteropoly  acid,  but  the  total  acidity  increases  with  the  H3PW  content.  TPD  [94]  and 
microcalorimetry  [41,  95]  with  NH3  indicate  that  H3PW  supported  on  SiOi  is  weaker  than 
the  pure  acid.  The  nature  of  interaction  between  H3PW  and  Si02  has  been  elucidated  by 
'H  and  ^'P  MAS  NMR.  No  difference  in  the  'H  NMR  spectra  of  Si02  and  Si02-H3PW 
supported  is  observed  up  to  50  wt.%  of  H3PW,  except  that  the  Si-OH  line  intensity  drops 
as  loading  of  H3PW  increases  [92].  It  is  concluded  that  at  concentrations  below  50  wt.% 
of  H3PW  over  Si02  the  proton  structure  is  different  from  that  crystalline  H3PW.  This  is 
also  confirmed  by  the  shift  in  the  ^'P  NMR  spectra  for  3-20  wt.%  of  H3PW  supported  on 
Si02  [92,  96].   The  interaction  of  Si02  with  H3PW  is  attributed  to  the  formation  of 
(=SiOH2"^)(H2PWi204o")  species  [96].    At  concentrations  up  to  20  wt.%  of  H3PW 
supported,  the  surface  species  are  isolated,  while  at  higher  than  50  wt.%  content  crystals 
of  free  H3PW  are  formed  on  the  silica  surface  [92].  Various  forms  of  the  polyanions  are 
observed  on  the  silica  surface  using  TEM  (transmission  electron  spectroscopy).  Isolated 
molecules  and  clusters  of  50  to  500  A  in  size  are  identified  [97],  depending  on  the  HPA 
loading. 

Alumina  is  usually  considered  as  unsuitable  as  a  support  because  it  decomposes 
the  Keggin  structure  of  HP  As  [27,  11].  This  decomposition  is  controversial  for 
H3PM012O40  [98,  99].  Studies  by  ^'P  NMR  indicate  that  the  interaction  between  H3PW 
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and  Y-AI2O3  is  not  strong  and  the  Keggin  structure  is  maintained.   If  Y-AI2O3-H3PW  ' 

supported  is  heated  above  200  °C,  the  ^'P  signal  broadness  suggesting  decomposition  of 

the  polyanion  structure  [98].  Comparison  of  the  catalytic  activity  of  H3PW  supported  on 

y-Al203  and  Si02  for  cumene  cracking,  toluene  disproportionation  and  n-hexane 

isomerization  indicate  that  centers  of  different  acid  strength  are  present  on  both  catalysts 

with  the  same  loading  of  heteropoly  acid  [100].  It  is  claimed  [100]  that  the  interaction  of 

H3PW  with  alumina  forms  aluminum  salts,  for  contents  up  to  10  wt.%  of  the  acid.  The 

aluminum  salts  of  H3PW  are  reported  [50]  to  contain  both  Lewis-acid  centers 

(coordinatively  unsaturated  aluminum  cations)  and  Bronsted-acid  sites.  The  IR  spectrum 

of  alumina-H3PW  supported  shows  bands  of  very  low  intensity  with  no  distinct  peaks  in 

the  region  of  Keggin  vibrations  [100]. 

Activated  carbon  entraps  HP  As  in  pores  forming  insoluble  compounds  that  do  not 
leach  any  acid  from  the  catalyst  support.  The  catalyst  can  be  used  for  liquid-solid  and 
gas-solid  organic  reactions  in  polar  media.  Leaching  of  HP  As  has  been  tested  by 
extraction  of  the  polyacids  with  water  and  methanol  [101].  The  fraction  of  HPA 
dissolved  out  of  the  carbon  support  (Calgon  F-300)  depends  upon  the  initial  loading  and 
the  solvent  employed,  but  a  final  concentration  of  7-14  wt.%  is  left  after  about  20  hs  of 
extraction.  Characterization  of  carbon-HsPW  supported  performed  by  microcalorimetry 
with  NH3  adsorption  [102]  shows  that  interaction  of  H3PW  with  activated  carbon  is 
strong.  H3PW  supported  on  carbon  (Darco)  is  a  weaker  acid  than  the  pure  acid,  and  the 
heat  is  dependent  upon  the  polyanion  loading.  It  is  claimed  [103]  that  supporting  H3PW 
on  various  oxides  (e.g.,  silica)  does  not  cause  such  a  drastic  decrease  in  acidity.  It  is 
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proposed  [104]  that  proton  transfer  occurs  from  the  heteropoly  acid  to  the  basic  groups 
(e.g.,  alcohols,  amines)  on  the  carbon  support.  Nonetheless,  results  of  microcalorimetry 
are  interpreted  assuming  only  coulombic  interactions  of  H3PW  with  the  surface  groups  on 
carbon  [103].  Potentiometric  titrations  and  ^'P  MAS  NMR  experiments  [103]  of  H3PW 
supported  on  carbon  confirms  no  degradation  of  the  Keggin  structure  when  it  is  deposited 
on  carbon  black. 

The  objective  of  this  research  is  to  characterize  the  acidity  of  H3PW  supported  on 
different  carries,  such  as  silica  gel,  alumina,  and  carbon  by  Cal-ad.  Cal-ad  determines  the 
number  of  sites  (nj),  equilibrium  constants  (Ki),  and  enthalpies  (AHj)  for  different  acid 
sites  on  these  solids.  Comparison  with  the  Cal-ad  results  of  free  H3PW  are  discussed. 

Experimental 
Materials  and  Preparation  of  HiPW  Supported 

Elemental  and  TG  analysis  of  the  H3PW  (purchased  from  Aldrich)  revealed  16 
moles  of  water  per  mole  of  H3PW.  Cyclohexane  obtained  from  Fisher  was  purified  by 
drying  over  4 A  molecular  sieves  for  24  hours,  and  then  distilled  over  P2O5.  Pyridine 
obtained  from  Fisher  was  distilled  over  CaHi.  Dried  samples  were  stored  in  a  container 
with  4A  molecular  sieves.  The  support  was  Davisil®  silica  gel  grade  62,  obtained  from 
Aldrich,  with  particles  within  60-200  mesh  and  a  surface  area  of  260  m^  g''.  Brokman® 
activated  acidic  aluminum  oxide  obtained  from  Fisher,  was  used  with  particles  within  80- 
200  mesh  and  a  surface  area  of  155  m^  g''.  The  support  supplied  by  Rohm  and  Haas  was 
Ambersorb®  572  with  100-200  mesh  particles  and  wdth  a  surface  area  of  1160  m^  g'. 
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The  H3PW  supported  catalysts  were  prepared  by  impregnation.  Aqueous  solutions  of 
H3PW  in  0.1  mol  L"'  of  HCl  were  prepared  using  concentrations  that  depend  upon  the 
desired  loading  of  the  support  (e.g.,  20  wt.%  H3PW).  The  support  was  added  to  a  volume 
of  solution  needed  to  form  a  suspension.  The  suspension  was  stirred  with  a  magnetic 
bar/stir  plate  and  evaporated  in  a  oil  bath  at  80  °C  to  dryness.  Then  the  solid  was  ground 
to  fine  particles  and  dried  in  a  reactor  under  vacuum.  No  pretreatment  of  supports  was 
carried  out,  except  for  Ambersorb®  572  which  was  washed  with  a  solution  of  0.1  mol  L"' 
of  HCl  for  6  hours  and  then  dried  at  100  °C  overnight  under  vacuum,  before  impregnation 
with  H3PW.  All  H3PW  supported  samples  were  dried  under  vacuum  at  200  °C  for  6 
hours,  except  carbon-H3PW  supported  which  was  dried  at  100  °C  for  12  hours. 

FTIR,  Calorimetric  and  Adsorption  Measurements 

Calorimetric,  FTIR,  and  adsorption  measurements  were  performed  on  the  H3PW 
supported  samples  according  to  the  procedures  described  in  Chapter  3.  The  analysis  by 
the  Cal-ad  method  was  also  explained  in  Chapter  3. 

Results  and  Discussion 
Cal-ad  Analysis  of  Silica-H^PW  Supported 

Calorimetric  titrations  with  pyridine  were  performed  on  H3PW  supported  on  silica 
with  different  loadings  of  the  polyacid.  Loadings  of  H3PW  ranging  from  8-24  wt.%  were 
studied.  The  level  of  coverage  to  produced  a  monolayer  in  SiOi  with  a  surface  area  of 
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300  g"'  is  about  20  wt.%  of  H3PW.  The  results  of  average  enthalpy  (calculated  based 
on  limiting  reagent  assumption)  for  reaction  of  these  solids  are  presented  in  Table  4-1 . 

Table  4-1.  Results  of  Calorimetric  Titrations  (25  ±  1  °C)  of  SiOi-HaPW  supported  with 
Pyridine. 

wt.%  H3PW    mmol  HsPW/g  of  solid    mmol  H^g  of  solid"    -AHavg  (kcal  mor')" 
24.5  0.084  0.252  27 

19.5  0.068  0.204  26 

15.3  0.053  0.159  24 

7.6  0.026  0.078  21 

a.  Value  based  on  the  amount  of  H3PW. 

b.  The  average  enthalpy  is  based  on  the  first  two  or  three  additions  of  pyridine  in  the 
calorimetric  titration. 

This  semi-quantitative  result  shows  that  the  interaction  of  H3PW  v^th  Si02  is  strong 
(AHavg  =  -26  kcal  mof')  and  almost  independent  of  the  heteropoly  acid  content  in  the 
range  fi-om  15-24  wt.%  of  H3PW.  Loading  8  wt.%  of  H3PW  gives  lower  enthalpy 
(AHavg  =  -21  kcal  mol'').  Based  upon  these  observations,  full  Cal-ad  analysis  was  made 
on  the  Si02-H3PW  supported  24.5  wt.%  sample.  The  calorimetric  and  adsorption 
isotherms  of  its  reaction  with  pyridine  in  cyclohexane  are  shown  in  Figure  4-1. 

The  Cal-ad  data  were  analyzed  assuming  the  presence  of  two  sites  and  three  sites. 
Three  site  fit  gives  a  meaningless  result  for  site  3  with  high  standard  deviation  errors  for 
the  parameters.  The  two  site  fit  gives  better  standard  deviations  on  the  parameters  and 
has  the  sum  of  squares  of  residuals  of  the  same  order  of  magnitude  as  the  experimental 
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Figure  4-1.  Isotherms  of  Calorimetry  and  Adsorption  of  Si02-H3PW  supported  (24.5 
wt.%)  Reacted  with  Pyridine  in  Cyclohexane. 
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error.  Table  4-2  reports  the  thermodynamic  parameters  for  the  interaction  with  pyridine 
of  H3PW  supported  on  siHca. 

The  enthalpy  of  interaction  with  pyridine  for  site  1  (AHi  =  -27.9  kcal  mol"')  is 
smaller  than  the  enthalpy  of  pure  H3PW  (AHi  =  -32.7  kcal  mol''),  but  the  number  of  sites 
are  slightly  higher  (0.10  mmol  g"')  for  H3PW  supported  on  silica  than  for  the  free  acid 
(0.08  mmol  g'').  The  second  site  enthalpy  is  10.1  kcal  mol''  and  a  number  of  sites  of  0.31 
mmol  g'',  while  free  H3PW  shows  AH2  =  -19.6  kcal  mol"'  with  na  =  0.16  mmol  g"'. 
These  results  suggest  that  H3PW  is  well  dispersed  on  silica  gel.  The  total  protons 
available  from  H3PW  (0.255  mmol  g'')  are  titrated  as  shown  by  the  sum  of  ni  and  n2 
(0.41  mmol  g"').  A  strong  interaction  between  H3PW  and  silica  is  indicated  by  'H  and 

T  1 

P  NMR  [92,  96].  At  concenfrations  up  to  20  wt.%  the  species  on  the  surface  are 
considered  isolated  [92].  Therefore  it  can  be  concluded  that  reaction  occurs  between 
H3PW  and  Si02  forming  a  weaker  acid  than  the  free  acid.  The  FTIR  spectnmi  (Figure  4- 
2)  shows  that  the  Keggin  structure  is  practically  intact  when  H3PW  is  supported  on  silica. 
After  reaction  with  pyridine  the  formation  of  pyridinium  ion  is  indicated  by  the  band  at 
1540  cm"'  [91].  Hydrogen-bonded  pyridine  is  also  indicated  by  a  band  at  1620  cm''. 

Cal-ad  Analysis  of  Altrniina-H^PW  Supported 

Cal-ad  results  for  alumina-H3PW  supported  are  shown  in  Table  4-2  while 
isotherms  from  calorimetry  and  adsorption  experiments  are  given  in  Figure  4-3.  The 
values  obtained  for  site  1  (AHi  =  -22.6  kcal  mol"',  ni  =0.15  mmol  g'')  and  site  2  (AH2  =  - 
12.0  kcal  mol'',  n2  =  0.20  mmol  g'')  demonstrate  this  is  a  weaker  solid  acid  than  H3PW. 
However,  the  quantity  of  sites  are  considerable  higher  than  that  of  pure  H3PW.  This 
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Table  4-2.  Cal-ad  Results  for  Si02''-H3PW  supported  (24.5  wt.%),  Al203''-H3PW 
supported  (19.4  wt.%),  and  Carbon'-HsPW  supported  (20.1  wt.%)  with  Pyridine  in 
Cyclohexane. 


H3PW/S1O2 

H3PW/AI2O3 

H3PW/Carbon 

ni  (mmol  g'') 

0.102  ±  0.005 

0.152  ±  0.007 

0.129  ±0.002 

Ki  (M  ) 

6.7x10  ±  1.3x10 

1        1  a5  I       1    1  a3 

1.2x10  ±2.3x10 

1.2x10  ±  3.2x10 

-Aril  (Kcal  mol  ) 

TT  Q  -1-  n  7 

11  .y  ±  U.  / 

ZZ.6  ±  U.z 

22.4  ±0.1 

(KCdLV  moi  ) 

7  Q 

/  .7 

O.y 

y./ 

A&i  (cai  Js.  mol  ) 

-00.  / 

^0  7 

/IT  0 

n2  (mmol  g"') 

0.31  ±0.06 

0.20  ±  0.04 

0.21  ±0.11 

K2  (M"') 

3.8x10'' ±7.2x10^ 

4.9x10^  ±8.7x10^ 

8.6x10^  ±4.8x10'' 

-AH2  (kcal  mol"') 

10.1  ±2.9 

12.0±  1.1 

6.0  ±4.9 

-AG2  (kcal  mol"') 

6.2 

5.0 

6.7 

AS2  (cal  K"'  mol"') 

-12.8 

-23.4 

2.4 

Total      (mmol  g'')" 

0.255 

0.202 

0.209 

a.  Davisil®  silica  gel  grade  62  with  surface  area  of  260  m^  g"'. 

b.  Brokman®  activated  acidic  aluminum  oxide  with  surface  area  of  155  m^  g"'. 

c.  Ambersorb®  572  carbonaceous  adsorbent  with  surface  area  of  1 160  m^  g"'. 

d.  Total  millimoles  of     is  based  on  the  H3PW  content. 


Figure  4-2.  FTIR  Spectrum  of  Si02-H3PW  supported  (24.5  wt.%). 
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Figure  4-3.  Isotherms  of  Calorimetry  and  Adsorption  of  AI2O3-H3PW  supported  (19.4 
wt.%)  Reacted  with  Pyridine  in  Cyclohexane. 
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suggests  that  H3PW  is  better  dispersed  on  alumina  than  siUca,  even  though  the  later 
carrier  has  a  higher  surface  area.  Other  evidence  of  good  dispersion  comes  from  the  total 
acidic  protons  presents  from  the  H3PW  loading.  About  75%  of  the  total  0.202  mmol  of 
protons  per  gram  of  solid  are  titrated  as  site  1 .  The  nature  of  this  site  is  quite  different 
from  pure  H3PW  or  alumina.  The  strength  and  quantity  of  sites  for  Y-AI2O3  [105]  are  AHi 
=  -15.2  kcal  mol'',  m  =  0.16  mmol  g"'  and  AH2  -  -6.7  kcal  mol'',  112  =  0.25  mmol  g"'  for 
site  1  and  2  respectively.  These  results  imply  that  reaction  of  H3PW  with  alumina 
produces  new  sites  stronger  than  pure  alumina  but  weaker  than  the  heteropoly  acid  non- 
supported. The  FTIR  spectrum  (Figure  4-4)  of  H3PW  supported  on  alumina  shows  that 
the  Keggin  bands  are  not  well  defined  and  shifts  of  the  bands  from  the  parent  acid 
suggests  some  degree  of  decomposition  of  the  heteropoly  acid  on  the  surface  of  alumina. 
After  reaction  of  H3PW  supported  on  alumina  with  pyridine,  FTIR  spectrum  shows  two 
bands,  one  corresponding  to  pyridinium  ion  at  1540  cm"'  and  one  assigned  to  weakly 
coordinated  pyridine  at  1630  cm"'. 

Cal-ad  Analysis  of  Carbon-HjPW  Supported 

Ambersorb®  572  (A-572)  load  with  20.1  wt.%  of  H3PW  reacts  with  pyridine  to 
give  the  calorimetric  and  adsorption  isotherms  shown  in  Figure  4-5.  The  results  of  Cal- 
ad  (Table  4-2)  shows  0.13  mmol  per  gram  of  solid  and  an  enthalpy  of -22.4  kcal  mol"'  for 
site  1 .  The  acid  H3PW  on  A-572  is  much  weaker  than  H3PW  (AHi  =  -32.7  kcal  mol"'  and 
ni  =  0.08  mmol  g"').  The  amount  of  sfrongest  sites  shows  that  H3PW  is  well  dispersed  on 
the  surface  of  the  carbon.  Although  the  supported  acid  is  weaker  than  the  pure  H3PW, 
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Figure  4-4.  FTIR  Spectrum  of  AI2O3-H3PW  supported  (19.4  wt.%) 
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Figure  4-5.  Isotherms  of  Calorimetry  and  Adsorption  of  Carbon-HsPW  supported  (20.1 
wt.%)  Reacted  with  Pyridine  in  Cyclohexane. 
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77%  of  the  0.202  mmol  of  protons  per  gram  of  solid  are  titrated  as  site  1 .  This  shows  an 
increase  in  relation  to  pure  H3PW  where  only  8%  of  1.04  mmol  of  protons  per  gram  are 
titrated  with  pyridine  as  site  1 .  For  reactions  demanding  solid  acids  with  enthalpies  about 
25  kcal  mol"'  H3PW  supported  on  carbon  is  a  viable  choice.  The  Cal-ad  results  help  to 
explain  why  some  reactions  show  better  activity  for  carbon-HsPW  supported  than  the 
unsupported  acid.  The  higher  surface  area  leads  to  a  higher  number  of  relatively  strong 
sites.  The  carbonaceous  adsorbent  used  (A-572)  has  been  characterized  using  gas-solid 
adsorption  with  nonpolar  gases  [106].  This  material,  a  macroreticular  sulfonated 
polystyrene  resin  which  is  pyrolyzed,  does  not  show  any  measurable  acidity.  Titration  of 
the  carbon  with  pyridine  yield  no  detectable  heat.  Therefore,  all  acidic  properties  of  the 
H3PW  supported  on  Ambersorb  572  carbon  is  provided  by  the  H3PW.  The 
thermodynamic  parameters  for  the  second  site  of  carbon-H3PW  supported  (AH2  =  -6.0 
kcal  mol'^  and  n2  =  0.21  mmol  g'')  correspond  to  weak  hydrogen-bonding  site. 
Physisorption  of  pyridine  by  H3PW  supported  on  A-572  can  be  seen  in  the  adsorption 
curve  (Figure  4-5).  This  contributes  to  the  low  enthalpy  of  site  2. 

Conclusions 

This  work  reports  thermodynamic  constants  for  H3PW  supported  on  silica, 
alumina,  and  carbon  reacting  with  pyridine  in  cyclohexane  solution.  Loadings  of  H3PW 
up  to  25  wt.%  on  these  carriers  decrease  the  acidity  of  the  free  acid.  However,  the  Cal-ad 
results  show  that  H3PW  is  well  dispersed  on  these  supports.  The  number  of  the  strongest 
sites  increases  compared  to  the  free  H3PW.   This  trend  does  not  follow  parallel  the 
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surface  area  of  the  different  supports,  as  expected.  The  amount  of  the  strongest  sites 
increases  in  the  order:  HaPW/alumina  >  HaPW/carbon  >  HaPW/silica,  while  the  surface 
area  is  carbon  »silica  >  alumina. 

The  acidity  order  is  consistent  with  the  literature  results  where  H3PW  >  silica- 
H3PW  supported  >  alumina-HaPW  supported  =  carbon-HaPW  supported.  This  acidity 
order  reflects  a  stronger  interaction  of  H3PW  with  alumina  and  carbon  than  with  silica. 
FTIR  spectrum  of  H3PW  supported  on  alumina  confirms  the  interaction  of  the  Keggin 
anion  with  the  silica  and  alumina  supports.  Shifts  of  the  Keggin  bands  and  even  loss  of 
the  bands  in  the  FTIR  spectrum  are  observed  for  H3PW  supported  on  alumina.  This 
explains  the  decrease  in  the  enthalpy  of  H3PW/AI2O3  catalyst  compared  to  H3PW. 


APPENDIX 


COMPUTER  PROGRAM  FOR  CALCULATION  OF  CALORIMETRIC  DATA  IN 

HOMOGENEOUS  SYSTEM 
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c 

C  FORTRAN 

C  PROGRAM:  CALORl. FOR 
C 
C 

C  PURPOSE: 

C  THIS  PROGRAM  WILL  CALCULATE  FROM  CALORIMETRIC  TITRATION  DATA 

C  THE  ENTHALPIES  OF  DEPROTONATION,  THE  EQUILIBRIUM  CONSTANTS  AND 

C  THE  RESPECTIVE  STANDARD  DEVIATIONS  FOR  SYSTEMS  OF  ACIDS  WITH 

C  ONE,  TWO,  OR  THREE  PROTONS  IN  NON  AQUEOUS  SOLUTIONS. 

Q*****Af*}lf*tAfAfMftfif*itMf********t***********************  **************** 

c 

C  REAL  VARIABLES: 

C 

C  B     (R*4)    =  BASE  CONCENTRATION 

C  BH    (R*4)    =  PROTONATED  BASE  CONCENTRATION 

C  DEFIND  (R*8,6)  =  DEFINED  VARIABLES  IN  THE  ORDER  Kl ,  K2,  K3, 

C  WTOT,  BTOT,  H. 

C  RESULT  (R*8,6)  =  CONCENTRATIONS  IN  ORDER  B,BH,WH3,WH2,WH1,WH0 
C 

Q*************************************************************************** 

C  AUTHOR:  TOM  MAIER  02/27/97  INITIAL  PROGRAM  (Version  1 .5) 

Q************************************************************** ************* 

C  THE  VERSION  OF  THIS  PROGRAM  FOR  IBM/PC  WAS  ADAPTED  BY  JOSE  A.  DIAS 

C  ON  MARCH  1997  USING  A  MICROSOFT  FORTRAN  POWER  STATION  VERSION  1 .0 

C  FOR  WINDOWS  3.1. 
C 

C  A  FILE  WITH  INSTRUCTIONS  ACCOMPANY  THE  DOCUMENTATION  OF  THIS  PROGRAM. 

Q*******1f*1f************************************************************^:^::tifif 


IMPLICIT  NONE 

REAL*8  Y(IOOO),  BN0N(6),  X(1000,6),  Z(IOOO),  VAR,  YH(IOOO) 
REAL*8  BASE(IOOO),  D(  1 000,6),  DY(IOOO),  HEATIN(IOOO),  VOL(IOOO) 
REAL*8  Kl,  K2,  K3,  HI,  H2,  H3,  SVART,  ACID(IOOO),  SSQE,  STD 
REAL*8  DUMMY,  KlIN,  K2IN,  K3IN,  HIIN,  H2IN,  H3IN,  STDE(6) 
INTEGER  I,  ICON,  IN(6),  ITMAX,  J,  NERR,  NSIG,  ND,  NPTOT,  NPF,  NPV 
INTEGER  NFIRST,  NLAST,  NPR,  IFIX(6),  INF(6),  CASE 
CHARACTER*5  CHAR*80,  VARY  I  (6),  VARY2(6),  VARY3(6) 
CHARACTER*80  INFILE,  OUTFILE 


DATA  IN  /  1,2,  3,  4,  5,6/ 

DATA  VARYl /'  KIV  HIV    V    V    V  '/ 
DATA  VARY2  / '  Kl','  K2  ','  HI  ','  H2  ','    ','  '/ 
DATA  VARY3  / '  Kl','  K2  ','  K3  ','  HI  ','  H2  ','  H3  '/ 


1  PRINT*," 

PRINT  *,'  THERE  ARE  THREE  CASES  WHICH  THIS  PROGRAM  WILL  HANDLE:' 
PRINT  *,' ' 

PRINT*,'  1  ONE  PROTON  ON  THE  ACID' 
PRINT*,'      AH  +  B  =  A  +  BH  ' 
PRINT  *,' ' 

PRINT  *,'  2  TWO  PROTONS  ON  THE  ACID' 
PRINT*,'      AH2+  B  =  AH+  BH  ' 
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PRINT  *;      AH  +  B  =  A  +  BH  • 

PRINT  *,'      AH2  +  2B  =  A  +  2BH  (TOTAL  REACTION)' 
PRINT  *; ' 

PRINT*;  3  THREE  PROTONS  ON  THE  ACID' 

PRINT  *;      AH3  +  B  =  AH2  +  BH  ' 

PRINT  *;      AH2  +  B  =  AH  +  BH  ' 

PRINT*,'      AH  +  B  =  A  +  BH  ' 

PRINT  *;      AH3  +  3B  =  A  +  3BH  (TOTAL  REACTION)' 

PRINT  *,' ' 

PRINT  *;  WHICH  CASE  DO  YOU  WANT  TO  USE?' 
READ(5,*)  CASE 

IF  (CASE  .LT.  1  .OR.  CASE  .GT.  3)  GO  TO  I 
NPR  =  0 
PRINT  *; ' 
C  OPEN  THE  INPUT  FILE. 
NPR  =  0 
WRITE  (6,3) 

3  FORMAT  (SX.-Name  of  Input  File?') 
READ  (5,'(A80)')  INFILE 
OPEN(UNIT=8,FILE=INFILE,STATUS='OLD') 
READ  (8,*)  NFIRST,  NLAST 

DO  I=NFIRST,NLAST 
J=I-NFIRST  +  1 

READ  (8,  *)  ACID(J),  BASE(J),  HEATIN(J),  VOL(J) 
END  DO 

ND  =  NLAST  -  NFIRST  +  I 
CLOSE  (8) 

C  THE  DATA  I  AM  USING  HAS  FOUR  COLUMNS  IN  THE  ORDER,  ACID,  BASE, 
C  HEAT  MEASURED,  AND  TOTAL  VOLUME  OF  THE  SOLUTION. 
C  THE  UNITS  ARE  MOL/L,  MOL/L,  CALORIES,  AND  LITERS  RESPECTIVELY. 
C  ONE  CAN  USE  DIFFERENT  UNITS  IF  CORRECTLY  DIMMENSIONED. 

C  OPEN  THE  OUTPUT  FILE. 

WRITE  (6,4) 

4  FORMAT  (5X, 'Name  of  Output  File?') 
READ  (5,'(A80)')  OUTFILE 

OPEN  (UNIT=43,FILE=OUTFILE ) 

C  START  THE  CALCULATIONS 

GO  TO  (20,  30,  40)  CASE 
20  PRINT*," 
PRINT  *,'  ENTER  AN  INITIAL  EST  OF  Kl' 
READ  (5,*)K1 
PRINT  *,' ' 

PRINT  *,'  ENTER  AN  INIT  EST  OF  ENTHALPY,  HI' 
PRINT  *,'  HI  IS  FOR  THE  REACTION  AH  +  B  ->  A  +  BH' 
PRINT  *,' ' 
READ  (5,*)  HI 
PRINT  *,' ' 

PRINT  *,'  HERE  ARE  THE  VALUES  YOU  ENTERED,  IN  THE  ORDER  YOU' 
PRINT  *,'  ENTERED  THEM.  IF  YOU  WANT  A  TERM  TO  VARY  FN  THE' 
PRINT  *,'  REGRESSION,  PUT  A  1  UNDER  THE  NUMBER.  IF  YOU  WANT' 
PRINT  *,'  IT  TO  STAY  THE  SAME  (NO  CHANGE  BY  THE  REGRESSION), ' 
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PRINT  *;  PUT  A  0  UNDER  THE  NUMBER.' 
PRINT  *; ' 

WRITE(6,FMT='(2X,6D10.2)')  Kl,  HI 
DO  1=3,6 

IFIX(I)  =  0 

END  DO 
READ  (5,*)  IFIX(l),  IFIX(2) 
GO  TO  50 
30  PRINT*," 

PRINT  *,'  ENTER  AN  INITIAL  EST  OF  Kl  AND  K2' 
READ  (5,*)K1,K2 
PRINT  *,' ' 

PRINT  *,'  ENTER  AN  INIT  EST  OF  EHTHALPIES  HI  AND  H2' 

PRINT  *,'  HI  IS  FOR  THE  REACTION  AH2  +  B  ->  AH  +  BH' 

PRINT  *,'  H2  IS  FOR  THE  REACTION  AH2  +  2B  ->  A  +  2BH' 

PRINT  *,' ' 

READ  (5,*)  H1,H2 

PRINT  *,' ' 

PRINT  *,'  HERE  ARE  THE  VALUES  YOU  ENTERED,  IN  THE  ORDER  YOU' 
PRINT  *:  ENTERED  THEM.  IF  YOU  WANT  A  TERM  TO  VARY  IN  THE' 
PRINT  *,'  REGRESSION,  PUT  A  1  UNDER  THE  NUMBER.  IF  YOU  WANT' 
PRINT  *,'  IT  TO  STAY  THE  SAME  (NO  CHANGE  BY  THE  REGRESSION), ' 
PRINT  *,'  PUT  A  0  UNDER  THE  NUMBER.' 
PRINT  *,' ' 

WRITE(6,FMT='(2X,6D10.2)')  Kl,  K2,  HI,  H2 
IFIX(5)  =  0 
IFIX(6)  =  0 

READ  (5,*)  IFIX(l),  IFIX(2),  IFIX(3),  IFIX(4) 
GO  TO  50 
40  PRINT*," 

PRINT  *,'  ENTER  AN  INITIAL  EST  OF  Kl,  K2,  K3' 
READ  (5,*)K1,K2,  K3 
PRINT  *,' ' 

PRINT  *,'  ENTER  AN  INIT  EST  OF  ENTHALPIES  HI,  H2,  H3' 
PRINT  *,'  HI  IS  FOR  THE  REACTION  AH3  +  B  ->  AH2  +  BH' 
PRINT  *,'  H2  IS  FOR  THE  REACTION  AH3  +  2B  ->  AHl  +  2BH' 
PRINT  *,'  H3  IS  FOR  THE  REACTION  AH3  +  3B  ->  A  +  3BH' 
PRINT  *,' ' 

READ  (5,*)HI,H2,  H3 
PRINT  *,' ' 

PRINT  *,'  HERE  ARE  THE  VALUES  YOU  ENTERED,  IN  THE  ORDER  YOU' 
PRINT  *,'  ENTERED  THEM.  IF  YOU  WANT  A  TERM  TO  VARY  IN  THE' 
PRINT  *,'  REGRESSION,  PUT  A  1  UNDER  THE  NUMBER.  IF  YOU  WANT' 
PRINT  *,'  IT  TO  STAY  THE  SAME  (NO  CHANGE  BY  THE  REGRESSION), ' 
PRINT  *,'  PUT  A  0  UNDER  THE  NUMBER.' 
PRINT  *,' ' 

WRITE(6,FMT='(2X,6D10.2)')  Kl,  K2,  K3,  HI,  H2,  H3 
READ  (5,*)  IFIX 
50  K1IN  =  K1 
K2IN  =  K2 
K3IN  =  K3 
H1IN  =  H1 
H2IN  =  H2 
H3IN  =  H3 
NPF  =  0 
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NPV  =  0 
DO  1=1,6 

IF  (IFIX(I)  .EQ.  0)  THEN 
NPF  =  NPF  +  1 
INF(NPF)  =  I 
END  IF 

IF  (IFIX(I)  .NE.  0)  THEN 
NPV  =  NPV  +  1 
IN(NPV)  =  I 
END  IF 
END  DO 
DO  1=1, NPF 
IN(NPV+I)  =  INF(I) 
END  DO 
1002  FORMAT(6A5) 
IF  (CASE  .EQ.  1)THEN 
IF  (NPV  .GT.  0)  THEN 
PRINT  *,' ' 

PRINT  *,'  THESE  ARE  THE  TERMS  WHICH  WILL  VARY:' 
PRINT  *,' ' 

WRITE(6, 1 002)  (VARY  1  (IN(I)),I=  1  ,NPV) 
END  IF 

IF  (NPF  .GT.  0)  THEN 
PRINT  *,' ' 

PRINT  *,'  THESE  ARE  THE  TERMS  WHICH  WILL  BE  FIXED:' 
PRINT  *,' ' 

WRITE(6, 1 002)  (VARY  1  (IN(NP V+I)),I=  1  ,NPF) 
END  IF 
END  IF 

IF  (CASE  .EQ.  2)  THEN 
IF  (NPV  .GT.  0)  THEN 
PRINT  *,' ' 

PRINT  *,'  THESE  ARE  THE  TERMS  WHICH  WILL  VARY:' 
WRITE(6, 1 002)  (VARY2(IN(I)),I=  1  ,NPV) 
PRINT  *,' ' 
END  IF 

IF  (NPF  .GT.  0)  THEN 
PRINT  *,' ' 

PRINT  *,'  THESE  ARE  THE  TERMS  WHICH  WILL  BE  FIXED:' 
PRINT  *: ' 

WRITE(6, 1 002)  (VARY2(IN(NPV+I)),I=1  ,NPF) 
END  IF 
END  IF 

IF  (CASE  .EQ.  3)  THEN 
IF  (NPV  .GT.  0)  THEN 
PRINT  *,' ' 

PRINT  *,'  THESE  ARE  THE  TERMS  WHICH  WILL  VARY:' 
PRINT  ' 

WRITE(6, 1 002)  (VARY3(IN(I)),I=  1  ,NPV) 
END  IF 

IF  (NPF  .GT.  0)  THEN 
PRINT  *,' ' 

PRINT  *,'  THESE  ARE  THE  TERMS  WHICH  WILL  BE  FIXED:' 
PRINT  *; ' 

WRITE(6, 1 002)  (VARY3(rN(NPV+I)),I=  1  ,NPF) 
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END  IF 
END  IF 
PRINT  *; ' 

PRINT  *;  HOW  MANY  ITERATIONS  DO  YOU  WANT  THE  REGRESSION  TO  TRY?' 

PRINT  *,'  IF  YOU  ENTER  0,  IT  WILL  DO  NO  REGRESSION,  BUT  IT  WILL' 

PRINT  *;  DO  THE  CALCULATIONS  WITH  YOUR  INPUT  VALUES.' 

PRINT  *;  IF  YOU  ENTER  A  NUMBER  GREATER  THAN  0,  IT  WILL  DO  THAP 

PRINT  *;  NUMBER  OF  ITERATIONS.' 

PRINT  *; ' 

READ  (5,*)  ITMAX 

PRINT  *,' ' 

NPTOT  =  6 

NSIG  =  3 

SVART=  1000000. 

DO  1=1, ND 

Y(I)  =  HEATIN(I) 

END  DO 

BN0N(3)  =  0.0 

BN0N(4)  =  0.0 

BN0N(5)  =  0.0 

BN0N(6)  =  0.0 

IF  (CASE  .EQ.  1)THEN 

BN0N(1)  =  K1 

BN0N(2)  =  H1 

END  IF 

IF  (CASE  .EQ.  2)  THEN 
BN0N(1)  =  K1 
BN0N(2)  =  K2 
BN0N(3)  =  H1 
BN0N(4)  =  H2 
END  IF 

IF  (CASE  .EQ.  3)  THEN 
BN0N(1)  =  K1 
BN0N(2)  =  K2 
BNON(3)  =  K3 
BN0N(4)  =  H1 
BN0N(5)  =  H2 
BN0N(6)  =  H3 
END  IF 
PRINT  *,'  ' 
NPR  =  0 

CALL  CHEAT(BNON,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 

SSQE=0.0 

DO  1=1, ND 

YH(I)=Z(I) 

DUMMY=  Y(I)-Z(I) 

SSQE=SSQE+DUMMY*DUMMY 

END  DO 

IF  (NPV.GT.O  .AND.  ITMAX.GT.O)  CALL  NLREG(NPTOT,ND,X,Y,D,BNON, 

*ICON,Z,YH,DY,NPF,NPV,IN,ITMAX,NSIG,VAR,BASE,ACID,VOL,CASE, 
■^STDE) 

PRINT  *,'  • 

IF  (ICON.GT.  0)  PRINT  *,'  NLREG  CONVERGED',ICON,'  ITERATIONS' 
IF  (ICON.EQ.-l)  PRINT  *,'  NLREG  STOPPED,  TOO  MANY  ITERATIONS' 
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IF  (IC0N.EQ.-2)  PRINT  *,'  NLREG  STOPPED,  XPX  MATRIX  SINGULAR' 
IF  (IC0N.EQ.-3)  PRINT  *,'  NLREG  STOPPED,  NO  BETTER  SOLUTION' 

IF  (CASE  .EQ.  1)THEN 
Kl  =  BNON(l) 
HI  =  BNON(2) 
END  IF 

IF  (CASE  .EQ.  2)  THEN 
Kl  =BN0N(1) 
K2  =  BN0N(2) 
HI  =BN0N(3) 
H2  =  BNON(4) 
END  IF 

IF  (CASE  .EQ.  3)  THEN 

Kl  =BNON(l) 

K2  =  BNON(2) 

K3  =  BN0N(3) 

HI  =  BNON(4) 

H2  =  BN0N(5) 

H3  =  BN0N(6) 

END  IF 

NPR  =  0 

CALL  CHEAT(BNON,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 

SSQE=0.0 

DO  1=1, ND 

YH(I)=Z(I) 

DUMMY=  Y(I)-Z(I) 

SSQE=SSQE+DUMMY'»DUM\fY 

END  DO 

VAR  =  SSQE/(ND-6) 
STD  =  DSQRT(VAR) 
PRINT  *,'  BEST  FIT  FOUND:' 
PRINT*,'  ' 

1000  FORMAT(A34,  2D  14.4) 

1001  FORMAT(34X,  A31) 
1003  FORMAT(/A34,2D14.4) 

IF  (CASE  .EQ.  1)THEN 

WRITE(6,1001)'    BEST  FIT    INITIAL  VALUE' 
WRITE(6,1000)'  AH  +  B->A  +  BH  K1=',K1,K1IN 
WRITE(6,1000)'  AH  +  B->A  +  BH  ENTHALPY  =',  H 1 ,  H 1  IN 
WRITE(6, 1 003) '  SUM  OF  SQUARES  ERROR        =',  SSQE 
WRJTE(6,1000)'  VARIANCE  =',  VAR 

WRITE(6,1000)'  STANDARD  DEVIATION  =',  STD 

WRITE(6, 1 003) '  STANDARD  ERROR  OF  K 1         =',  STDE(  1 ) 
WRITE(6, 1 000) '  STANDARD  ERROR  OF  ENTHALPY    =',  STDE(2) 

WRITE(43 , 1 00 1 )  •    BEST  FIT    INITIAL  VALUE' 
WRITE(43,1000)'  AH  +  B->A  +  BH  K1=',K1,K1IN 
WRITE(43,1000)'  AH  +  B ->  A  +  BH  ENTHALPY  =',H1,  HI  IN 

WRITE(43,1003)'  SUM  OF  SQUARES  ERROR  =',SSQE 

WRITE(43,1003) '  STANDARD  ERROR  OF  Kl         =',  STDE(l) 
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WRITE(43,1003) '  STANDARD  ERROR  OF  ENTHALPY  =',  STDE(2) 
END  IF 

IF  (CASE  .EQ.  2)  THEN 

WRITE(6, 1 00 1 ) '    BEST  FIT    INITIAL  VALUE' 
WRITE(6,1000)'  AH2+  B->AH1+  BH  K1=',K1,KIIN 
WRITE(6,1000)'  AHl  +  B->  A  +  2BH      K2  =',  K2,  K2IN 
WR1TE(6,I003)'  AH2+  B->AH1+  BH  ENTHALPY     HI,  H UN 
WRITE(6,1000)  •  AHI  +  2B  ->  A  +  2BH  ENTHALPY  =',  H2,  H2IN 
WRITE(6,I003)'  SUM  OF  SQUARES  ERROR        =',  SSQE 
WRITE(6,1000)'  VARIANCE  =',  VAR 

WRITE(6,1000) '  STANDARD  DEVIATION  =',  STD 

WRITE(6, 1 003)  •  STANDARD  ERROR  OF  K 1         =',  STDE( I ) 
WRITE(6,1000)'  STANDARD  ERROR  OF  K2         =',  STDE(2) 
WRITE(6,1003)  •  STANDARD  ERROR  OF  ENTHALPY  1  STDE(3) 
WRITE(6,1000)  •  STANDARD  ERROR  OF  ENTHALPY  2  =',  STDE(4) 

WRITE(43,1001)'    BEST  FIT   INITIAL  VALUE' 
WRITE(43,1000)'  AH2  +  B->AH1+  BH      Kl  =',K1,  KlIN 
WRITE(43,I000) '  AHl  +  B  ->  A  +  BH      K2  =',K2,  K2IN 

WRITE(43,1003)'  AH2  +  B->AH1+  BH  ENTHALPY  =',H1,  HI  IN 
WRITE(43,1000) '  AH2  +  2B  ->  A  +  2BH  ENTHALPY  =',H2,  H2IN 

WRITE(43,1003)'  SUM  OF  SQUARES  ERROR  =',SSQE 

WRITE(43, 1 003) '  STANDARD  ERROR  OF  K I  =',  STDE(  1 ) 
WRITE(43 , 1 000) '  STANDARD  ERROR  OF  K2         =',  STDE(2) 

WRITE(43,1003) '  STANDARD  ERROR  OF  ENTHALPY  1  =',  STDE(3) 
WRITE(43,1000) '  STANDARD  ERROR  OF  ENTHALPY  2  =',  STDE(4) 
END  IF 

IF  (CASE  .EQ.  3)  THEN 

WRITE(6, 1 00 1 ) '    BEST  FIT    INITIAL  VALUE' 
WRITE(6,1000) '  AH3  +  B  ->  AH2  +  BH      Kl  =',  Kl,  KlIN 
WRITE(6,1000) '  AH2  +  B  ->  AHI  +  BH      K2  =',  K2,  K2rN 
WRITE(6,1000) '  AHl  +  B  ->  A  +  BH      K3  =',  K3,  K3IN 
WR1TE(6,1003)'  AH3+  B ->  AH2  +  BH  ENTHALPY  =',  HI,  HI  IN 
WR1TE(6,1000) '  AH3  +  2B  ->  AHl  +  2BH  ENTHALPY  =',  H2,  H2IN 
WRITE(6,1000) '  AHl  +  3B  ->  A  +  3BH  ENTHALPY  =',  H3,  H3IN 
WRITE(6, 1 003) '  SUM  OF  SQUARES  ERROR         =',  SSQE 
WRITE(6,1000)'  VARIANCE  =',  VAR 

WRITE(6,I000)'  STANDARD  DEVIATION  =',  STD 

WRITE(6,1003) '  STANDARD  ERROR  OF  Kl         =',  STDE(l) 
WRITE(6, 1 000) '  STANDARD  ERROR  OF  K2         =',  STDE(2) 
WRITE(6,1000)'  STANDARD  ERROR  OF  K3         =',  STDE(3) 
WRITE(6,1003) '  STANDARD  ERROR  OF  ENTHALPY  1   =',  STDE(4) 
WR1TE(6,1000)'  STANDARD  ERROR  OF  ENTHALPY  2  =',  STDE(5) 
WRITE(6, 1 000) '  STANDARD  ERROR  OF  ENTHALPY  3  =',  STDE(6) 

WRITE(43 , 1 00 1 )  '    BEST  FIT    INITIAL  VALUE' 
WRITE(43,1000) '  AH3  +  B  ->  AH2  +  BH      Kl  =',K1,  KlIN 
WRITE(43,1000) '  AH2  +  B  ->  AHl  +  BH      K2  =',K2,  K2IN 
WRITE(43,1000) '  AHl  +  B  ->  A  +  BH      K3  =',K3,  K3IN 
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WRITE(43,1003)'  AH3  +  B ->  AH2  +  BH  ENTHALPY  =',H1,  HI  IN 
WRITE(43,1000) '  AH3  +  2B  ->  AHl  +  2BH  ENTHALPY  =',H2,  H2IN 
WRITE(43, 1000) '  AH3  +  3B  ->  A  +  3BH  ENTHALPY  =',H3,  H3IN 

WRITE(43,1003)'  SUM  OF  SQUARES  ERROR  =',SSQE 

WRITE(43,1003) '  STANDARD  ERROR  OF  Kl  =',  STDE(l) 
WRITE(43 , 1 000) '  STANDARD  ERROR  OF  K2  =',  STDE(2) 
WRITE(43,1000)'  STANDARD  ERROR  OF  K3         =',  STDE(3) 

WRITE(43,1003) '  STANDARD  ERROR  OF  ENTHALPY  1  =',  STDE(4) 
WRITE(43, 1 000) '  STANDARD  ERROR  OF  ENTHALPY  2  =',  STDE(5) 
WRITE(43, 1 000) '  STANDARD  ERROR  OF  ENTHALPY  3  =',  STDE(6) 
END  IF 

IF  (ITMAX  .GT.  0) 

*WRITE(43,1003)'  VARIANCE  =',  VAR 

IF  (ITMAX  .GT.  0) 

* WRITE(43 , 1 000) '  STANDARD  DEVIATION  STD 
IF  (ITMAX  .EQ.  0) 

*WRITE(43,FMT='(A32,D14.4)') '  NO  REGRESSION  WAS  DONE.  ' 
WRITE(43,FMT='(//,A32,/)') '  RECONSTRUCTION  OF  INPUT  DATA: ' 
WRITE(43,FMT='(A36,A13  )•)'  POINT    MEASURED  PREDICTED' 

*  ,  •  DIFFERENCE  ' 
DO  1=1, ND 

WRITE(43,FMT='(I6,3F14.6)')  I,  Y(I),  Z(I),  (Y(I)-  Z(I)) 
END  DO 

IF  (CASE  .EQ.  1)THEN 

WRITE(43,FMT='(//,5X,A4,8X,A4,6X,A9,3X,A9,2(6X,A3,3X))') 
*'ACID','BASE','FREE  BASE'.'PROT  B ASE','AH I'.'AHO' 
END  IF 

IF  (CASE  .EQ.  2)  THEN 

WRITE(43,FMT='(//,4X,A4,6X,A4,3X,A9,2X,A9,3(4X,A3,3X))') 
*'ACID','BASE',TREE  BASE','PROT  BASE','AH2VAH I'/AHO' 
END  IF 

IF  (CASE  .EQ.  3)  THEN 

WRITE(43,FMT='(//,4X,A4,5X,A4,4X,A9,2X,A9,4(4X,A3,3X))') 
*'ACID','BASE','FREE  BASE','PROT  BASE';AH3','AH2VAH1VAH0' 
END  IF 
NPR=  I 

CALL  CHEAT(BN0N,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 

CL0SE(43) 

STOP 

END 


c============================================================, 

C  SUBROUTINE  CHEAT(BNON,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 
C 

C  PURPOSE 

C    THIS  SUBROUTINE  SETS  UP  THE  CALCULATION  FOR  EACH  TITRATION 

C     POINT  AND  CALCULATES  THE  HEAT  OF  REACTION. 

C 

C  VARIABLES 

C    BNON    (R*8,6)  =  KS  AND  HS,  DEPENDING  ON  CASE. 
C    Z      (R*8,1000)=  PREDICTED  HEATS 
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C  ND     (1*4)    =  NUMBER  OF  DATA  POINTS 

C  ACID   (R*8, 1 000)=  TOTAL  ACID  CONCENTRATIONS 

C  BASE    (R*8, 1 000)=  TOTAL  BASE  CONCENTRATIONS 

C  CASE    (1*4)    =  1  =  1  PROTON  ACID,  2=2  PROTONS,  3=3  PROTONS 

C  A      (R*8)    =  ACID  WITHOUT  ANY  PROTONS 

C  AH     (R*8)    =  ACID  WITH  ONE  PROTON 

C  AH2    (R*8)   =  ACID  WITH  TWO  PROTONS 

C  AH3     (R*8)    =  ACID  WITH  THREE  PROTONS 

C  B      (R*8)    =  FREE  BASE 

C  BH     (R*8)    =  PROTONATED  BASE 

C  DEFIND  (R*8,6)  =  DEFINED  VARIABLES  IN  THE  ORDER  Kl,  K2,  K3, 
C  WTOT,  BTOT,  H. 

C  RESULT  (R*8,6)  =  CONCENTRATIONS  IN  ORDER  B,  BH,  WH3,  WH2,  WHl,  W 
C 

C================================================================ 

SUBROUTINE  CHEAT(BNON,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 
IMPLICIT  NONE 

REAL*8  DEFIND(6),  RESULT(6),  WTOT,  WHO,  WHl,  WH2,  WHS 
REAL  EXPJ 

REAL*8  ATOT,  BTOT,  ONE,  TWO,  THREE,  FOUR 

REAL*8  BN0N(6),  Z(  1 000),  BASE(  1 000),  ACID(  1 000),  VOL(  1 000) 

REAL*8  Kl,  K2,  K3,  A,  AH,  AH2,  AH3,  B,  BH,  C,  D,  E 

REAL*8  DELI,  DEL2,  BETA,  BESAVE(2),  PLOST,  PONB 

INTEGER  CASE,  I,  J,  ND,  NPR 

DATA  ONE  /l.ODOO/ 

DATA  TWO  /2.0D00/ 

DATA  THREE  /  3.0D00/ 

DATA  FOUR  /  4.0D00/ 

IF  (CASE  .EQ.  1)THEN 

Kl  =BN0N(1) 

DO  1=1, ND 

BTOT  =  BASE(I) 

ATOT  =  ACID(I) 

C  THE  CASE  OF  ONE  PROTON  CAN  BE  SOLVED  BY  A  QUADRATIC  EQUATION, 

C     E*A**2  +  D*A  +  C  =  0 

C    A  =  (-D  +-  SQRT(D*D  -  4*E*C))  /  (2*E) 

E  =  0NE-K1 
D  =  (ATOT  +  BTOT)  *  Kl 
C  =  -Kl  *  ATOT  *  BTOT 
C  =  DSQRT(D*D  -  FOUR  *  E  *  C) 
A  =  (-D  +  C)  /  (TWO  *  E) 
IF  (A.LT.0.0  .OR.  A.GT.ATOT  .OR.  A.GT.BTOT) 

*  A  =  (-D  -  C)  /  (TWO  *  E) 

IF  (A.LT.0.0  .OR.  A.GT.ATOT  .OR.  A.GT.BTOT) 

*PRINT  *,'  ERROR  IN  CALCULATING  CONCENTRATIONS  ON  POINT 

*  I 
BH  =  A 

AH  =  ATOT  -  A 
B  =  BTOT  -  BH 
Z(I)  =  VOL(I)  *  A  *  BN0N(2) 
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IF  (NPR  .EQ.  1)  WRITE(43,FMT='(6D12.3)')  ACID(I),  BASE(I), 

*  B,  BH,  AH,  A 
END  DO 
RETURN 

END  IF 

C  IN  CASE  2,  THERE  ARE  FIVE  EQUATIONS  FOR  FIVE  UNKNOWNS,  AH2,  AH,  A, 

C    B,  AND  BH  (ASSUMING  Kl  AND  K2  ARE  KNOWN).  I  REDUCED  THE  FIVE 

C    EQUATIONS  IN  FIVE  UNKNOWNS  TO  TWO  EQUATIONS  IN  TWO  UNKNOWNS. 

C    HERE  IS  AN  ATTEMPT  TO  SOLVE  FOR  THOSE  TWO  UNKNOWNS.  THE  KEY  IS 

C    TO  ASSUME  BETA  (=  BH/B,  WHICH  ALLOWS  THE  CALCULATION  OF  B  AND  BH) 

C    CALCULATE  AH,  AND  FROM  AH  CALCULATE  AH2  AND  A.  THEN  CHECK 

C    WHETHER  THE  ASSUMED  BETA  IS  CORRECT  BY  DOING  A  MASS  BALANCE  ON 

C    THE  PROTONS  (THE  PROTONS  ON  BH  MUST  HAVE  BEEN  THE  ONES  LOST  FROM 

C    AH  AND  AH2).  WHEN  THE  PROTONS  BALANCE,  THE  EQUATIONS  HAVE  BEEN 

C  SOLVED. 

IF  (CASE  .EQ.  2)  THEN 
Kl  =BN0N(1) 
K2  =  BNON(2) 
DO  I=I,ND 
BTOT  =  BASE(I) 
ATOT  =  ACID(I) 
BESAVE(1)  =  0.0 
BESAVE(2)  =  0.0 
DEL2  =  0.00000  IDOO  *  ATOT 
DO  J=l,50 
IF  (J  .LT.  3)  THEN 
EXPJ  =  20. 

IF  (J  .EQ.  2)  EXPJ  =  -20. 
EXPJ  =  IO.O**EXPJ 
BETA  =  EXPJ 
END  IF 

IF  (J  .GT.  2)  BETA  =  DSQRT(BESAVE(1)  *  BESAVE(2)) 
BH  =  BTOT  /  (ONE  +  ONE/BETA) 
B  =  BH  /  BETA 

AH  =  ATOT  /  (ONE  +  BETA/Kl  +  K2/BETA) 
A  =  K2*AH/BETA 
AH2  =  AH*BETA/K1 
PLOST  =  TWO*A  +  AH 
PONB  =  BH 

DELI  =  DABS(PLOST  -  PONB) 

IF  (PLOST  .GT.  PONB)  BESAVE(l)  =  BETA 

IF  (PLOST  .LT.  PONB)  BESAVE(2)  =  BETA 

IF  (DELI  .LT.  DEL2)  GO  TO  10 

END  DO 

PRINT  *,'  WARNING,  POINT',  I, '  DID  NOT  CONVERGE.' 

10       Z(I)  =  VOL(I)  *  (AH  *  BNON(3)  +  A  *  BNON(4)) 

IF  (NPR  .EQ.  1)  WRITE(43,FMT='(7D10.3)')  ACID(I),  BASE(I), 

*  B,  BH,  AH2,  AH,  A 
END  DO 

RETURN 
END  IF 
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C  IN  CASE  3  THIS  PROGRAM  WILL  USE  A  PREVIOUS  SOBROUTINE  WHERE  THE  SPECIES 
C  CONCENTRATIONS  ARE  CALCULATED  SEPARETEDLY  BY  OTHER  SUBROUTINE 
"CONCEN" 

IF  (CASE  .EQ.  3)  THEN 
Kl  =BN0N(1) 
K2  =  BNON(2) 
K3  =  BN0N(3) 
DO  1=1, ND 

BTOT  =  BASE(I) 
WTOT  =  ACID(I) 

IF  (ABS(BTOT-3.0*WTOT)  .LT.  0.01* WTOT)    GO  TO  100 

DEFIND(1)  =  BN0N(1) 

DEFIND(2)  =  BN0N(2) 

DEFIND(3)  =  BN0N(3) 

DEFIND(4)  =  WTOT 

DEFIND(5)  =  BTOT 

DEFIND(6)  =  3.0D00*WTOT 

IF  (I  .EQ.  1)THEN 

RESULT(l)  =  0.0000000001 

RESULT(2)  =  BTOT 

RESULT(3)  =  0.00000000001 

RESULT(4)  =  0.00000000001 

RESULT(5)  =  0.000000001 

RESULT(6)  =  WTOT 

END  IF 

IF  (I  .GT.  1)  RESULT(2)  =  BTOT 

IF  (BTOT-3.0D0*WTOT  .GT.  0.0)  RESULT(2)  =  3.0D0*WTOT 

CALL  CONCEN  (RESULT,  DEFIND,  I) 

B  =RESULT(1) 
BH  =RESULT(2) 
WHO  =  RESULT(3) 
WHl  =RESULT(4) 
WH2  =  RESULT(5) 
WH3  =  RESULT(6) 

Z(I)  =  VOL(I)  *  (WH2*BNON(4)  +  WH1*BN0N(5)  +  WH0*BNON(6)) 
IF(NPR.EQ.  1) 

*  WRITE(43,FMT='(8D10.3)')  ACID(I),  BASE(I),  B,  BH,  WH3,  WH2, 

*  WHl,  WHO 


100  CONTINUE 
END  DO 
RETURN 
END  IF 

END 

C============================ 

C  SUBROUTINE  CONCEN(RESULT,  DEFIND,  NP) 
C 

C  PURPOSE 
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C    GIVEN  ALL  THE  DEFINED  VARIABLES,  THIS  ROUTINE  CALCULATES  THE 
C     CONCENTRATIONS  OF  B,  BH,  WH3,  WH2,  WH 1 ,  WHO.  THIS  IS  AN 
C    ITERATIVE  SEARCH  ROUTINE. 

C 

C  VARIABLES 

C    DEFIND  (R*8,6)  =  DEFINED  VARIABLES  IN  THE  ORDER  Kl,  K2,  K3, 
C  WTOT,  BTOT,  H. 

C    RESULT  (R*8,6)  =  CONCENTRATIONS  IN  ORDER  B,  BH,  WH3,  WH2,  WHl,  WHO 
C    DEFI  (R*4,4)=  INPUT  VARIABLES  K1/K2,  WTOT,  H/K2/K3 
C    DEFE  (R*4,4)  =  ESTIMATES  OF  THE  VARIABLES  IN  DEFI 
C 

C  TOM  MAIER  2/96 

SUBROUTINE  CONCEN(RESULT,  DEFIND,  NP) 
IMPLICIT  NONE 

REAL*8  DEFIND(6),  RESULT(6),  DELTA(4),  WF,  WFF 
REAL*8  B,  BH,  WH3,  WH2,  WHl,  WHO,  Kl,  K2,  K3,  K12,  K23 
REAL*8  DEFI(4),  DEFE(4) 

REAL*8  WTOT,  BTOT,  HTOT,  ERRAIM(4),  DELW 

REAL*8  WHON,  WHIN,  WH2N,  WH3N,  F012,  F123,  BHB,  DENOM 

INTEGER  I,  J,  KOUNT,  NL012,  NL123,  NP 

INTEGER  NF012(1 155),  NF123(1 155) 

KOUNT  =  0 

F012  =  l.ODOO 

F123  =  l.ODOO 

B  =RESULT(1) 

BH  =RESULT(2) 

WHO  =  RESULT(3) 

WHl  =RESULT(4) 

WH2  =  RESULT(5) 

WH3  =  RESULT(6) 

Kl  =DEFIND(1) 

K2  =  DEFIND(2) 

K3  =DEFrND(3) 

K12  =  K1/K2 

K23  =  K2/K3 

WTOT  =  DEFIND(4) 

BTOT  =  DEFIND(5) 

HTOT  =  DEFIND(6) 

DEFI(1)  =  K12 

DEFI(2)  =  WTOT 

DEFI(3)  =  HTOT  -  BTOT 

DEFI(4)  =  K23 

IF  (DEFI(3)  .LT.  0.0)  THEN 

DEFI(3)  =  HTOT 

BH  =  3.0*WTOT 

B  =  BTOT - BH 

BHB  =  BH/B 

DENOM  =  1  .ODOO  +  BHB/K3  +  BHB*BHB/(K3*K2)  + 
*  BHB*BHB*BHB/(K3*K2*K1) 
WHO  =  WTOT/DENOM 
WHl  =  WH0*BHB/K3 
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WH2  =  WH1*BHB/K2 
WH3  =  WH2*BHB/K1 
RESULT(1)  =  B 
RESULT(2)  =  BH 
RESULT(3)  =  WHO 
RESULT(4)  =  WHl 
RESULT(5)  =  WH2 
RESULT(6)  =  WH3 
RETURN 
END  IF 

ERRAIM(1)=  1.0D-05*DEFI(1) 
ERRAIM(2)  =  1.0D-04*DEFI(2) 
ERRAIM(3)=  I.0D-04*DEFI(3) 
ERRAIM(4)  =  1.0D-05*DEFI(4) 
IF  (ERRAIM(3)  .LT.  l.OD-08)  ERRAIM(3)  =  l.OD-08 


C  ESTIMATE  THE  KNOWN  VALUES,  DEFI  ARE  INPUT,  DEFE  ARE  ESTIMATES. 
C  DEFI(l)  =  K1/K2,  DEFI(2)  =  WTOT,  DEFI(3)  =  H-BH,  DEFI(4)  =  K2/K3 

10  KOUNT  =  KOUNT+  1 
WF  =  WH3  +  WH2  +  WHl  +  WHO 
DEFE(l)  =  WH2*WH2/(WH1*WH3) 
DEFE(2)  =  WH3  +  WH2  +  WHl  +  WHO 
DEFE(3)  =  3.0D00*WH3  +  2.0D00*WH2  +  WHl 
DEFE(4)  =  WH1*WH1/(WH2*WH0) 

C  ADJUST  TOTAL  W  BY  MULTIPLYING  CALCULATED  CONCENTRATIONS  BY  WF. 

WF  =  DEFI(2)  /  DEFE(2) 
WH3  =  WF  *  WH3 
WH2  =  WF  *  WH2 
WHl  =  WF*  WHl 
WHO  =  WF  *  WHO 

C  ADJUST  BASE  CONCENTRATIONS  BY  MULTIPLYING  BY  WF. 

IF  (KOUNT  .GT.  2)  THEN 
WF  =  BTOT  /  (B  +  BH) 
B  =  WF  *  B 
BH  =  WF  *  BH 
DEFI(3)  =  HTOT  -  BH 
END  IF 

C  ADJUST  TOTAL  H  CONCENTRATION  BY  ADJUSTING  WHi  CONCENTRATIONS. 

DEFE(3)  =  3.0D00*WH3  +  2.0D00*WH2  +  WHl 
WF  =  DEFI(3)  -DEFE(3) 
15  1  =  0 

IF  (DABS(WF)  .LT.  l.OD-10)  I  =  1 
IF  (WF  .GT.  0.0)  THEN 
IF  (WF  .LT.  WHO)  THEN 
1=1 

WHO  =  WHO  -  WF 
WHI  =  WHl  +  WF 
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GO  TO  20 
END  IF 

IF  (WF.LT.  WH1)THEN 
1=1 

WHl  =  WHl  -  WF 
WH2  =  WH2  +  WF 
GO  TO  20 
END  IF 

IF  (WF  .LT.  WH2)  THEN 
1=1 

WH2  =  WH2  -  WF 
WH3  =  WH3  +  WF 
GO  TO  20 
END  IF 
END  IF 

IF  (WF  .LT.  0.0)  THEN 
WFF  =  -  WF 

IF  (WFF.LT.  WH3)  THEN 
1=1 

WHS  =  WH3  -  WFF 
WH2  =  WH2  +  WFF 
GO  TO  20 
END  IF 

IF  (WFF.LT.  WH2)  THEN 
1=1 

WH2  =  WH2  -  WFF 
WHl  =  WHl  +  WFF 
GO  TO  20 
END  IF 

IF  (WFF.LT.  WHl)  THEN 
1=1 

WHl  =  WHl  -  WFF 
WHO  =  WHO  +  WFF 
GO  TO  20 
END  IF 
END  IF 
20  IF  (I  .EQ.  0)  THEN 
WF  =  0.8D0*WF 
GOTO  15 
END  IF 

C  CHECK  TO  SEE  HOW  WELL  WE  HAVE  DONE. 

IF  (WHS  .LE.0.0)  WHS  =  l.OD-10 
IF  (WH2  .LE.  0.0)  WH2  =  l.OD-10 
IF  (WHl  .LE.  0.0)  WHl  =  l.OD-10 
IF  (WHO  .LE.  0.0)  WHO  =  l.OD-10 
DEFE(l)  =  WH2*WH2/(WH1*WH3) 
DEFE(2)  =  WH3  +  WH2  +  WHl  +  WHO 
DEFE(3)  =  3.0D00*WH3  +  2.0D00*WH2  +  WHl 
DEFE(4)  =  WH1*WH1/(WH2*WH0) 
DO  1=1,4 

DELTA(I)  =  DEFI(I)  -  DEFE(I) 
END  DO 

IF  (DABS(DELTA(1))  .LT.  ERRAIM(l)  .AND. 
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*  DABS(DELTA(2))  .LT.  ERRAIM(2)  .AND. 

*  DABS(DELTA(3))  .LT.  ERRAIM(3)  .AND. 

*  DABS(DELTA(4))  .LT.  ERRAIM(4))  THEN 
RESULT(1)  =  B 

RESULT(2)  =  BH 
RESULT(3)  =  WHO 
RESULT(4)  =  WHl 
RESULT(5)  =  WH2 
RESULT(6)  =  WHS 
RETURN 
END  IF 

C  THE  ABOVE  CALCULATIONS  LINEARLY  BALANCED  CONCENTRATIONS.  NOW 

C  MAKE  ADJUSTMENTS  BASED  ON  THE  EQUILIBRIUM  CONSTANTS,  ACTUALLY  THE 

C  THE  RATIOS  OF  THE  Ks:  K1/K2  and  K2/K3  (THESE  RATIOS  ELIMINATE  B/BH) 

NL012  =  0 

IF  (WH2  .LT.  WHO)  NL012  =  2 
NLI23  =  1 

IF  (WHS  .LT.  WH I )  NL 1 23  =  3 

IF(NL123  .EQ.  1)THEN 

WHIN  =  WH2*WH2/(WH3*KI2) 

DELW  =  F12S  *  (WHI  -  WHIN) 

IF  (DELW  .GT.  0.6D0*WH1)  DELW  =  0.6D0*WH1 

IF  (DELW  .GT.  0.6D0*WH3)  DELW  =  0.6D0*WHS 

IF  (-2.0D0*DELW  .GT.  0.6D0*WH2)  DELW  =  -O.SDO*WH2 

WHl  =  WHl  -  DELW 

WH2  =  WH2  +  DELW  +  DELW 

WHS  =  WHS  -  DELW 

END  IF 

IF  (NL123  .EQ.  3)  THEN 

WHSN  =  WH2*WH2/(WH1*K12) 

DELW  =  F12S  *  (WHS  -  WHSN) 

IF  (DELW  .GT.  0.6D0*WH1)  DELW  =  0.6D0*WH1 

IF  (DELW  .GT.  0.6D0*WHS)  DELW  =  0.6D0*WHS 

IF  (-2.0D0*DELW  .GT.  0.6D0*WH2)  DELW  =  -O.SDO*WH2 

WHl  =  WHl  -  DELW 

WH2  =  WH2  +  DELW  +  DELW 

WHS  =  WHS  -  DELW 

END  IF 

IF  (NL012  .EQ.  0)  THEN 

WHON  =  WH1*WH1/(WH2*K2S) 

DELW  =  F012  *  (WHO  -  WHON) 

IF  (DELW  .GT.  0.6D0*WH0)  DELW  =  0.6D0*WH0 

IF  (DELW  .GT.  0.6D0*WH2)  DELW  =  0.6D0*WH2 

IF  (-2.0D0*DELW  .GT.  0.6D0*WH1)  DELW  =  -O.SDO*WHI 

WHO  =  WHO  -  DELW 

WHl  =  WHl  +  DELW  +  DELW 

WH2  =  WH2  -  DELW 

END  IF 

IF  (NL0I2  .EQ.  2)  THEN 

WH2N  =  WHl  *WH1/(WH0*K2S) 

DELW  =  F012  *  (WH2  -  WH2N) 

IF  (DELW  .GT.  0.6D0*WH0)  DELW  =  0.6D0*WH0 

IF  (DELW  .GT.  0.6D0*WH2)  DELW  =  0.6D0*WH2 

IF  (-2.0D0*DELW  .GT.  0.6D0*WHI)  DELW  =  -O.SDO*WHI 
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WHO  =  WHO  -  DELW 

WHl  =  WHl  +  DELW  +  DELW 

WH2  =  WH2  -  DELW 

END  IF 

IF  (WH3.GT.WH2)  B  =  WH2*BH/(WH3*K1) 
IF  (WH2.GT.WH3)  B  =  WH1*BH/(WH2*K2) 
IF  (WHO.GT.WHI)  B  =  WH0*BH/(WH1*K3) 

DEFE(l)  =  WH2*WH2/(WH1*WH3) 
DEFE(2)  =  WH3  +  WH2  +  WHl  +  WHO 
DEFE(3)  =  3.0DOO*WH3  +  2.0D00*WH2  +  WHl 
DEFE(4)  =  WH1*WH1/(WH0*WH2) 
DO  1=1,4 

DELTA(I)  =  DEFI(I)  -  DEFE(I) 
END  DO 

NF012(KOUNT)=  1 

IF  (DELTA(4)  .LT.  0.0)  NF012(KOUNT)  =  -1 
NF123(KOUNT)=  1 

IF  (DELTA(l)  .LT.  0.0)  NF123(KOUNT)  =  -1 

IF  (KOUNT  .GT.  2)  THEN 
J  =  NF012(KOUNT-1)*NF012(KOUNT) 
IF  (J.LT.O  .AND.  F0I2.GT.0.001)  F012  =  0.8D00*F012 
J  =  NF123(KOUNT-l)*NF123(KOUNT) 
IF  (J.LT.O  .AND.  FI23.GT.0.001)  FI23  =  0.8D00*F123 
END  IF 

C  CHECK  TO  SEE  HOW  WELL  WE  HAVE  DONE. 

IF  (DABS(DELTA(I))  .LT.  ERRAIM(l)  .AND. 

*  DABS(DELTA(2))  .LT.  ERRAIM(2)  .AND. 

*  DABS(DELTA(3))  .LT.  ERRAIM(3)  .AND. 

*  DABS(DELTA(4))  .LT.  ERRAIM(4))  THEN 

RESULT(1)  =  B 
RESULT(2)  =  BH 
RESULT(3)  =  WHO 
RESULT(4)  =  WHl 
RESULT(5)  =  WH2 
RESULT(6)  =  WH3 
RETURN 
END  IF 
DO  1=1,4 

DELTA(I)  =  DELTA(I)  /  ERRAIM(I) 
END  DO 

IF  (KOUNT.LT.  1 1 50)  GO  TO  1 0 

55  PRINT  *,'WARNING,  SLOW  CONVERGENCE  ON  POINT',  NP 
RESULT(1)  =  B 
RESULT(2)  =  BH 
RESULT(3)  =  WHO 
RESULT(4)  =  WHl 
RESULT(5)  =  WH2 
RESULT(6)  =  WH3 
RETURN 
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END 

C  SUBROUTINE  NLREG  (NPTOT,  ND,  X,  Y,  D,  BNON,  ICON,  Z,  YH,  DY,  NPF,  NPV, 

C  IN,  ITMAX,  NSIG,  VAR,  BASE,  ACID,  VOL,  CASE,  STDE) 

C 

C  ***  NON-LINEAR  REGRESSION  SUBROUTINE  *** 

SUBROUTINE  NLREG  (NPTOT,  ND,  X,  Y,  D,  BNON,  ICON,  Z,  YH,  DY,  NPF, 

♦  NPV,  IN,  ITMAX,  NSIG,  VAR,  BASE,  ACID,  VOL,  CASE,  STDE) 

C     IMPLICIT  NONE 

REAL*8  X(ND,NPTOT),  Y(ND),  BNON(l),  D(ND,NPTOT),  Z(ND),  YH(ND) 

REAL*8  BASE(ND),  DY(ND),  ACID(ND),  VOL(ND),  STDE(6) 

REAL*8  XPX(6,6),  XPY(6),  C0REC(6),  C0R(12),  C(6) 

REAL*8  DD,  DEL,  G,  GF,  GG,  QBEST,  QNEW,  QOLD,  SIG,  SUM,  TP,  VAR 

INTEGER  CASE,  I,  ICON,  ITER,  J,  NCOR,  ND,  NPF,  NPR,  NPV,  NSIG 

INTEGER  ICOR,  IN(NPTOT),  ITMAX,  K,  KI,  KJ,  NBEST,  NPTOT 

DATA  COR/1 .0D0,0.5D0,0.2D0,0.4D0,0.6D0,0.8D0, 

*  0.02D0,0.04D0,0.06D0,0.1D0,-0.1D0,-0.5D0/,NCOR/12/ 

c 

C  INPUT  ARGUMENTS  ••♦ 

c 

C  NPTOT  TOTAL  NUMBER  OF  NON-LINEAR  PARAMETERS. 

C  ND  TOTAL  NUMBER  OF  OBSERVED  DATA  POINTS. 

C  X  ND  X  NPTOT  MATRIX. 

C  Y  ND  VECTOR  OF  RESPONSES. 

C  D  MATRIX  OF  SIZE  AT  LEAST  ND  X  NPTOT. 

C  BNON  VECTOR  OF  LENGTH  NPTOT 

C  ON  INPUT,  CONTAINS  INITIAL  VALUES  FOR  NON-LINEAR 

C  PARAMETERS.  ON  OUTPUT,  CONTAINS  ESTIMATES  OF 

C  PARAMETERS. 

C  ICON  CONVERGENCE  INDICATOR  RETURNED  FROM  NLREG: 

C  ICON=NO.  OF  ITERATION  (  CONVERGENCE  ATTAINED) 

C  =-1  TOO  MANY  ITERATIONS  (  HAS  NOT  CONVERGED) 

C  =-2  XPX  MATRIX  SINGULAR 

C  =-3  NO  IMPROVEMENT. 

C  Z  WORK  VECTOR  OF  SIZE  AT  LEAST  ND. 

C  YH  PREDICTED  VALUES  RETURNED  FROM  NLREG. 

C  NPF  NUMBER  OF  NON-LINEAR  PARAMETERS  FIXED  (NOT  VARYING). 

C  NPV  NUMBER  OF  NON-LINEAR  PARAMETERS  VARYING. 

C  IN  NPTOT  VECTOR  INDICATING  INDICES  OF  NON-LINEAR  PARAMETERS 

C  VARYING. 

C  ITMAX  MAXIMUM  NUMBER  OF  ITERATIONS  ALLOWED.  ON  OUTPUT,  ITMAX 

C  CONTAINS  THE  ACTUAL  NUMBER  OF  ITERATIONS. 

C  NSIG  NUMBER  OF  SIGNIFICANT  DIGITS  DESIRED  IN  PARAMETER 

C  ESTIMATES. 

C  VAR  ESTIMATED  VARIANCE  OF  OBSERVATIONS  RETURNED  FROM  NLREG 
C 

C  STDE  STANDARD  DEVIATION  FOR  EACH  CALCULATED  PARAMETER 
C 

C     *****  NOTE:  ***** 
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C    *****    THIS  SUBROUTINE  IS  CURRECTLY  SET  UP  FOR  A  MAXIMUM  OF 
C    *****    6  PARAMETERS. 

Q  *♦*****************»*****♦♦♦**♦♦*****♦♦♦♦♦♦***♦*♦*************»»***** 

SIG=NSIG 
TP=10.0**(-SIG) 
ICON=0 
ITER=0 
DO  1=1,6 
COREC(I)=0.0D0 
END  DO 
100  NPR  =  0 

CALL  CHEAT(BNON,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 

QOLD=0.0 

DO  1=1, ND 

YH(I)  =  Z(I) 

DEL  =  Y(I)  -  YH(I) 

QOLD  =  QOLD  +  DEL*DEL 

END  DO 

IF  (M0D(ITER,4)  .EQ.  1)  THEN 
PRINT  *,'  • 

PRINT  * '  ***************************************************' 
WRITE(6,FMT='(3X,A21,14,3X,A23,D12.5)')  'AT  END  OF  ITERATION  ', 

*  ITER, '  SUM  OF  SQUARES  ERROR  =',  QOLD 
IF  (CASE  .EQ.  1)  PRINT  *,'  Kl  AND  HI  ARE:' 

IF  (CASE  .EQ.  2)  PRINT  *,'  Kl,  K2,  HI,  AND  H2  ARE:' 

IF  (CASE  .EQ.  3)  PRINT  *,'  Kl,  K2,  K3,  HI,  H2,  AND  H3  ARE:' 

WRITE(6,FMT='(2X,6D1 1.3)')  (BN0N(I),I=1,2*CASE) 

PRINT  *  '  ***************************************************' 

PRINT  *,' ' 
END  IF 

CALL  DERV  (X,Y,BNON,D,Z,YH,NPF,NPV,IN,NPTOT,ND,BASE,ACID, 

*  VOL,  CASE) 

CALL  XPRMX  (D,ND,ND,NPV,XPX) 

DO  J=1,ND 

DY(J)  =  Y(J)  -  YH(J) 

END  DO 

CALL  XPRMY  (D,DY,ND,NPV,XPY) 
CALL  INVERT  (XPX,NPV,DD) 
IF  (DD.NE.O.ODO)  GO  TO  120 

C  ***  XPX  MATRIX  IS  SINGULAR  -  IC0N=-2  -  ZERO  OUT  MATRICES  AND  RETURN  * 

DO  110I=1,NPTOT 
DO  110  J=1,NPT0T 
110D(J,I)=0.D0 
VAR=0.0 
DO  1=1, ND 

VAR=VAR+(Y(I)-YH(I))*(Y(I)-YH(I)) 
END  DO 

VAR=VAR/(ND-NPV) 

IC0N=-2 

RETURN 
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C  ***  END  XPX  MATRIX  SINGULAR  SECTION  ***** 

Q  +  *  + 

C  ***  FINISH  IF  ICON  IS  NOT  EQUAL  TO  ZERO  OR  ITMAX  IS  ZERO 

120  IF  (ICON  .NE.  0)  GO  TO  200 
IF  (ITMAX  .EQ.  0)  GO  TO  200 
C  ***  COREC  IS  THE  AMOUNT  BY  WHICH  TO  ADJUST  EACH  BNON  ELEMENT. 
DO  1=1, NPV 
KJ=IN(1) 
SUM=0.0D0 
DO  J=1,NPV 

SUM=SUM  +  XPX(I,J)*XPY(J) 
END  DO 

COREC(KJ)=SUM 
END  DO 
1TER=1TER+1 
GF=  1.8D0 
DO  1=1, NPV 

IF  (COREC(l)  .NE.  0.0)  THEN 

GG  =  DABS(BNON(I)/COREC(I)) 

IF  (GG  .LT.  GF)  GF  =  GG 

END  IF 

END  DO 

GF  =  0.5D0*GF 

C  ***  PARTIAL  CORRECTION  ROUTINE  FOR  PARAMETERS 

C  ***  TRY  THE  FRACTIONAL  CORRECTIONS  CONTAINED  IN  THE  VECTOR  COR. 

NBEST  =  0 

QBEST  =  QOLD 

DO  1C0R=1,NC0R 

G  =  GF*C0R(1C0R) 

DO  K=1,NPT0T 

C(K)=BNON(K)+G*COREC(K) 

END  DO 

NPR  =  0 

CALL  CHEAT(C,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 

QNEW=0.0 

DO  1=1, ND 

YH(1)=Z(I) 

QNEW=QNEW+(Y(1)-YH(I))*(Y(I)-YH(I)) 
END  DO 

IF  (QNEW  .LT.  QBEST)  THEN 
NBEST  =  ICOR 
QBEST  =  QNEW 
END  IF 
END  DO 

IF  (NBEST  .GT.  0)  THEN 

G  =  GF*COR(NBEST) 

DO  K=l,NPTOT 

C(K)=BNON(K)+G*COREC(K) 

END  DO 

GO  TO  130 

END  IF 
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C  ♦**  ALL  CORRECTIONS  TRIED;  NO  IMPROVEMENT 
C***  DEFINE  IC0N=-3 

ICON  =  -3 

GO  TO  200 

C  ***  CHECK  PARAMETERS  FOR  CONVERGENCE  *** 
130  DO  J=l,NPTOT 
IF  (BNON(J).EQ.O.ODO)  GO  TO  140 
IF  (DABS((C(J)-BNON(J))/C(J)).GT.TP)  GO  TO  150 
END  DO 

C  CONVERGED,  DEFINE  ICON  =  NUMBER  OF  ITERATIONS  USED 

140  ICON  =  ITER 
150  DO  J=l,NPTOT 

BNON(J)  =  C(J) 

END  DO 

IF  (ICON.GT.O)  GO  TO  200 

IF  (ITER  .LT.  ITMAX)  GO  TO  100 

C  **♦  TOO  MANY  ITERATIONS  -  SET  ICON  =  -1 

IC0N  =  -1 

C  CALCULATE  STD  ERRORS  FOR  ALL  PARAMETERS 
200  NPR  =  0 

CALL  CHEAT(C,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 
VAR  =  0.0 
DO  1=1, ND 
YH(I)=Z(I) 

VAR=VAR+(Y(I)-YH(I))*(Y(I)-YH(I)) 
END  DO 

VAR  =  VAR  /  (ND-NPV) 

DO210I=l,ND 

DO210  J=1,NPT0T 
210D(I,J)  =  0.0 

DO  215  1=1,6 

D0  215  J=l,6 
215  XPX(I,J)  =  0.0 

C  CALCULATE  XPX  MATRIX 

CALL  DERV(X,Y,BNON,D,Z,YH,NPF,NPV,IN,NPTOT,ND,BASE,ACID,VOL,CASE) 
CALL  XPRMX  (D,ND,ND,NPV,XPX) 
DO  220  l=l,NPTOT 
DO  220  J=1,I 
220  D(J,1)  =  O.ODO 

CALL  INVERT  (XPX,NPV,DD) 
IF  (DD  .EQ.  O.ODO)  THEN 
WRITE  (6,FMT='(/,A40,/)') 

*'  XPX  MATRIX  SINGULAR  ON  FINAL  ITERATION' 

RETURN 

END  IF 
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DO230I=l,NPV 
KI  =  IN(I) 
DO  230  J=1,I 
KJ  =  IN(J) 

D(KJ,KI)  =  XPX(I,J)*VAR 
IF  (KI  .EQ.  KJ)  THEN 
STDE(KI)  =  DSQRT(D(KJ,KI)) 

PRINT      KI,  BNON,  STDE  =',  KI,  BNON(KI),  STDE(KI) 
END  IF 
230  CONTINUE 


RETURN 
END 


SUBROUTINE  INVERT  (XPX,NROW,DD) 
IMPLICIT  NONE 
REAL*8  XPX(6,6),  AI(36),  DD 
INTEGER  I,  J,  K,  L(6),  M(6),  NROW 


C  SUBROUTINE  TO  INVERT  A  MATRIX 


K  =  0 

DO  100  1=1, NROW 
DO  100  J=1,NR0W 
K  =  K+  1 
100  AI(K)  =  XPX(I,J) 
CALL  DMINV  (AI,NROW,DD,L,M) 
K  =  0 

DO  200  1=1, NROW 
DO  200  J=1,NR0W 
K  =  K+  1 
200  XPX(I,J)  =  AI(K) 
RETURN 
END 


SUBROUTINE  XPRMY  (X,Y,NR,NCOL,XPY) 

IMPLICIT  NONE 
REAL*8  X(NR,NCOL),  Y(NR),  XPY(NCOL) 
INTEGER  NCOL,  NR,  I,  J 


C  SUBROUTINE  TO  CALCULATE  XPY  MATRIX 


DO  1=1, NCOL 
XPY(I)  =  O.ODO 
DO  J=1,NR 

XPY(I)  =  XPY(I)  +  Y(J)*X(J,I) 

END  DO 

END  DO 

RETURN 

END 

SUBROUTINE  XPRMX  (X,ND,NR,NCOL,XPX) 
C     IMPLICIT  NONE 

REAL*8  X(ND,NCOL),  XPX(6,6),  SUM 
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INTEGER  I,  J,  K,  NCOL,  ND,  NR 

C  CALCULATE  XPX  MATRIX 

DO  10  1=1, NCOL 
DO  10  J=1,I 
SUM  =  O.ODO 
DO  K=1,NR 

SUM  =  SUM  +  X(K,I)*X(K,J) 
END  DO 
XPX(I,J)  =  SUM 
10XPX(J,I)  =  SUM 
RETURN 
END 

C============================================================= 

SUBROUTINE  DERV  (X,Y,BNON,D,Z,YH,NPF,NPV,IN,NPTOT,ND,BASE, 
*ACID,VOL,CASE) 

C  THIS  SUBROUTINE  CALCULATES  DERIVATIVES  AS  DELTA  YH/DELTA  BNON 
C  FOR  EACH  DATA  POINT  (YH)  AND  EACH  BNON  VARIABLE. 

C     IMPLICIT  NONE 

REAL*  8  X(ND,NPTOT),  Y(l),  BNON(l),  D(ND,NPTOT),  YH(ND) 
REAL*8  Z(ND),  BASE(ND),  ACID(ND),  VOL(ND),  YHSAVE(IOOO) 
REAL*8  BHOLD,  DEL,  DELTA,  VAR 

INTEGER  CASE,  I,  IN(1),  J,  KI,  ND,  NPF,  NPR,  NPTOT,  NPV 

DEL=0.1D-3 

DO  J=1,ND 

YHSAVE(J)  =  YH(J) 

END  DO 

DO  100  1=1, NPV 

KI  =  IN(I) 

BHOLD  =  BNON(KI) 

DELTA  =  DEL  *  BNON(KI) 

IF  (ABS(DELTA)  .LT.  0.000001)  DELTA=0.000001 

BNON(KI)  =  BNON(KI)  +  DELTA 

NPR  =  0 

CALL  CHEAT(BNON,  Z,  ND,  BASE,  ACID,  VOL,  NPR,  CASE) 

VAR  =  0.0 

DO  J=1,ND 

YH(J)  =  Z(J) 

DEL  =  Y(J)  -  YH(J) 

VAR  =  VAR  +  DEL*DEL 

END  DO 

DO  J=1,ND 

D(J,I)  =  (YH(J)  -  YHSAVE(J))  /  DELTA 
END  DO 

BNON(KI)  =  BHOLD 
100  CONTINUE 

DO  1=1, ND 
YH(I)  =  YHSAVE(I) 
END  DO 
RETURN 
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END 

Q    **********************************************  t*****************^*;^* 

SUBROUTINE  DMINV(A,N,D,L,M) 
IMPLICIT  NONE 

REAL*8  DABS,  A(l),  D,  BIGA,  HOLD 

INTEGER  I,  U,  IZ,  J,  JI,  JP,  K,  KI,  KK,  L(l),  M(l),  N,  NK 

INTEGER  IK,  JK,  JQ,  JR,  KJ 

C       SEARCH  FOR  LARGEST  ELEMENT 

D=I.O 

NK=-N 

DO  80  K=I,N 

NK=NK+N 

L(K)=K 

M(K)=K 

KK=NK+K 

BIGA=A(KK) 

DO  20  J=K,N 

IZ=N*(J-I) 

DO  20  I=K,N 

IJ=IZ+I 

IF(DABS(BIG  A)-DABS(A(IJ)))  1 5,20,20 
15  BIGA=A(IJ) 

L(K)=I 

M(K)=J 
20  CONTINUE 

C       INTERCHANGE  ROWS 

J=L(K) 

IF(J-K)  35,35,25 
25  KI=K-N 

DO  30  1=1, N 

KI=KI+N 

HOLD=-A(KI) 

JI=KI-K+J 

A(KI)=A(JI) 
30  A(JI)  =HOLD 

C      INTERCHANGE  COLUMNS 

35  I=M(K) 

IF(I-K)  45,45,38 
38  JP=N*(I-I) 

DO  40  J=1,N 

JK=NK+J 

JI=JP+J 

HOLD=-A(JK) 

A(JK)=A(JI) 
40  A(JI)  =HOLD 


C  DIVIDE  COLUMN  BY  MINUS  PIVOT  (VALUE  OF  PIVOT  ELEMENT  IS 
C       CONTAINED  IN  BIGA) 
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45  IF(BIGA)  48,46,48 

46  D=0.0 
RETURN 

48  DO  55  1=1, N 

IF(I-K)  50,55,50 
50  IK=NK+I 

A(IK)=A(IK)/(-BIGA) 
55  CONTINUE 

C       REDUCE  MATRIX 

DO  65  1=1, N 
IK=NK+I 
HOLD=A(IK) 
IJ=I-N 

DO  65  J=1,N 
IJ=IJ+N 

IF(I-K)  60,65,60 
60  IF(J-K)  62,65,62 
62  KJ=IJ-I+K 

A(IJ)=HOLD*A(KJ)+A(IJ) 
65  CONTINUE 

C      DIVIDE  ROW  BY  PIVOT 

KJ=K-N 

DO  75  J=1,N 

KJ=KJ+N 

IF(J-K)  70,75,70 
70  A(KJ)=A(KJ)/BIGA 
75  CONTINUE 

C      REPLACE  PIVOT  BY  RECIPROCAL 

A(KK)=1.0/BIGA 
80  CONTINUE 

C      FINAL  ROW  AND  COLUMN  INTERCHANGE 

K=N 
100  K=(K-1) 

IF(K)  150,150,105 
105 I=L(K) 

IF(I-K)  120,120,108 
108  JQ=N*(K-1) 

JR=N*(I-1) 

DO  110  J=1,N 

JK=JQ+J 

HOLD=A(JK) 

JI=JR+J 

A(JK)=-A(JI) 
110  A(JI)  =HOLD 
120  J=M(K) 

IF(J-K)  100,100,125 


125  KI=K-N 
DO  130  I=1,N 
KI=KI+N 
HOLD=A(KI) 
JI=KI-K+J 
A(KI)=-A(JI) 

130  A(JI)  =HOLD 
GO  TO  100 

150  RETURN 
END 
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